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NOMENCLATURE 


A         Skempton's  pore  pressure  coefficient 

A         alternate  form  of  Skempton's  pore  pressure 
coefficient  given  by  A  =  AB 

A.p        Skempton's  pore  pressure  coefficient 
at  failure 

AA        change  in  area 

a         interface  adhesion  at  soil/construction 
material  interface  (psi  or  N/m2) 

B         Skempton's  pore  pressure  parameter  defined 
by  B  =  [(Au)/(Aa3)] 

C         capacitance  (Farads) 

C_        shape  factor  for  hydraulic  flow  through 
a        soil  medium 

c'        soil  cohesion  in  terms  of  effective  stress 
(psi  or  N/m2) 

D         diameter 

2 
E        modulus  of  elasticity  (psi  or  N/m  ) 

o 
E,        dynamic  elastic  modulus  (psi  or  N/m  ) 

e        void  ratio 

F(u,v)    Fourier  transform  representing  the  complex 
amplitude  as  a  function  of  position 
coordinates  (u,v)  in  the  transform  plane 
(back  focal  plane) 

f         focal  length  of  lens 

f(x,y)  input  function  as  a  function  of  position 
coordinates  (u,v)  in  the  transform  plane 
(back  focal  plane) 

f         interface  shear  resistance  factor  for 
cohesion  f  =  a/c  (psi  or  N/m2) 

f  /        interface  shear  resistance  factor  for  friction 
0  f /  =  5/t    (psi  or  N/m2) 

f         characteristic  frequency  (Hz) 
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NOMENCLATURE  (Cont'd.) 

2 
G,        dynamic  shear  modulus  (psi  or  N/m  ) 

H         height 

AH  fluid  head  loss 

I    (k)  a  Bessel   function  of  the   first  kind  of 

pure   imaginary  argument 

i  D.C.    current    flow    (amps) 

k  dielectric   constant 

K  active  earth  pressure  coefficient   that 


a 


defines  the  stress  ratio  at  failure 


ku        hydraulic  permeability  coefficient 

K        coefficient  of  earth  pressure  at  rest 
o 

K        passive  earth  pressure  coefficient  that 
P        defines  the  stress  ratio  at  failure 

n        porosity 

P_        internal  pressure  (psi  or  N/m  ) 

Pr        overconsolidation  ratio 

Ps        instantaneous  shut  in  pressure  for 
hydraulic  fracture  (psi  or  N/m^) 

Q        hydraulic  flow  rate  (volume  per  unit  time) 

R„        hydraulic  radius 

r        number  of  spikes  between  D.C.  (zero 

frequency)  and  the  first  order  spike  with 
height—  some  constant  specified  fraction 
of  the  height  of  the  first  order  spike, 
for  load  a. 

r,        same  as  for  r  ,  for  load  b 

r_        mean  value  of  r  for  both  sides  of  the 
profiles 

r,        mean  value  of  r,  for  both  sides  of  the 

D  c.   -i  D 

profile  s 
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NOMENCLATURE  (Cont'd.) 

o 
S         shear  strength  (psi  or  N/m  ) 

o 
S         vertical  shear  strength  (psi  or  N/m  ) 

2 
S,        horizontal  shear  strength  (psi  or  N/m  ) 

2 
S..       total  stress  tensor  (psi  or  N/m  ) 

T         torque  required  to  shear  vanes  of  various 
aspect  ratio  (lb-ft  or  N-m) 

2 
u         pore  pressure  in  the  sample  (psi  or  N/m  ) 

2 
u         initial  pore  pressure  at  depth  Z  (psi  or  N/m  ) 

V  measured  velocity  (ft/sec  or  m/sec) 

V  dilatational  wave  velocity  (ft/sec  or  m/sec) 

Vg        distortional  wave  velocity  (ft/sec  or  m/sec) 
AV        impressed  voltage  (volts) 
x3y,z     cartesian  coordinates 

a         angle  coordinate  of  center  of  gravity  of 
°  distribution 

a         angle  coordinate  of  center  of  gravity  of  sample 

taken  as  maximum  likelihood  estimate  of  « 

g 

3      /  3 

V  density  (lb/ft   or  kg/m  ) 


Y 


w 


ev»  eh 


unit  weight  of  water 


5         interface  friction  angle 

*         Kronecker's  delta 

°PY 

a  Kronecker's   delta 

ctY 


vertical  and  horizontal  strains 
€         permittivity  constant  (Farads /meter) 
t)         viscosity  of  fluid 
k        a  measure  of  concentration  about  the  mean 
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NOMENCLATURE  (Cont'd.) 

\  wavelength  of  incident  light  (meter) 

v         Poisson's  ratio 

vd        dynamic  Poisson's  ratio 

3        3 
p         mass  density  (slug/ ft   or  kg/m  ) 

2 
(jTiGojG'i     principal  stresses  (psi  or  N/m  ) 

g        confining  pressure 

q. .       tensor  representing  the  three  dimensional 
L-J       state  of  stress  at  a  point 

a  mean  stress 

°m 

a  normal  stress 

°n 


av>crh 


vertical  and  horizontal  normal  stress  (psi  or  N/iri) 


a  ,ct  normal  stress  in  x  and  z  directions 

Qnr  DC  conductivity  (mho/m) 

cr.p  AC  conductivity  (mho/m) 

0-',  q/  effective  vertical  and  horizontal  normal  stress 

tt  Interface  shear  resistance 

T   ,  i  q  shear  stress  in  xz  or  r0  plane 

0  angle  of  internal  friction 

0'  angle  of  friction  based  on  effective  stress 
notation 

i!r  critical  angle  to  vertical  for  crack  propagation 
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DETERMINATION  OF  IN-SITU  STRESSES 
IN  SOIL  MASSES 

1 .   INTRODUCTION 

The  behavior  of  a  soil  mass  under  loading  or  during 
excavation  depends  upon  the  stress  field  that  naturally 
exists  in-situ.  This  stress  field  defines  the  extent  to 
which  the  soil  shear  strength  may  be  already  mobilized  and 
the  origin  of  the  stress  paths  that  will  be  followed  on  sub- 
sequent additional  loading,  unloading  and  changes  of  the 
boundaries  of  the  soil  masses.   Usually,  the  in-situ  state  of 
stress  at  a  point  in  the  soil  mass  is  not  the  result  of  elastic 
or  plastic  equilibrium  alone  but  also  is  influenced  by  the 
stress  history  and  other  factors  (e.g.,  chemical,  temperature). 

At  present  computer  methods  allow  us  to  estimate 
ground  movements  and  changes  in  the  state  of  stress  during 
and  after  the  construction  of  structures.   The  reliability 
of  the  computation  and,  in  particular,  the  evaluation  of 
stability  conditions,  the  estimation  of  ground  movements  and 
stresses  on  the  structures  depend  equally  upon  a  realistic 
soil  model  and  the  initial  stress  state.   A  realistic  soil 
model  itself  requires  knowledge  of  all  the  stresses,  including 
the  intermediate  stresses.   The  significance  of  the  inter- 
mediate stress  on  soil  behavior  has  been  pointed  out  by 
Henkel  (1958),  Duncan  (1966)  and  others.   The  factor  of 
safety  depends  on  the  initial  stress  state,  the  degree  to 
which  the  soil  strength  is  mobilized  and  the  origin  of  the 
subsequent  loading  path.   It  is  necessary  that  the  full  in-situ 
stress  tensor  existing  in  the  field  be  known. 

Movements  predicted  by  modern  mathematical  techniques 
are  sensitive  to  the  initial  stresses  input  into  the  program. 
Current  practice  is  to  make  simplifying  assumptions  regarding 
the  existing  in-situ  stresses,  and  under  some  geotechnical 
conditions  these  assumptions  are  justified.   The  assumptions 


often  are  hidden  within  the  programming,  so  that  the 
geotechnical  engineer  may  not  be  aware  that  initial  in-situ 
stress  levels  are  based  on  estimated  values  of  elastic  moduli. 
In-situ  stresses  should  be  measured  by  geotechnical  engineers- - 
not  estimated  by  computers. 

This  study  addresses  the  problem  of  the  in-situ 
state  of  stress  in  soils.   Under  simple  conditions,  such  as 
a  uniform,  normally  consolidated  soil  under  a  horizontal 
ground  surface,  the  state  of  stress  can  be  determined  from 
theoretical  considerations.   In  other  cases,  it  should  be 
measured  in  the  field  or  in  laboratories . 

A  number  of  field  methods  for  measuring  in-situ 
stresses  exists:  hydraulic  fracturing,  borehole  pressuremeter , 
anisotropic  vane  shear  test,  acoustic  velocity  measurements 
and  others.   The  use  of  these  techniques  permits  us  to  make 
measurements  in-situ  of  the  least  principal  stress,  the  mean 
stress  in  the  horizontal  plane,  or  the  stress  ratio  cru/av- 
In-situ  stresses  may  be  also  estimated  by  laboratory  tests 
in  combination  with  some  knowledge  of  the  stress  history  of 
the  site.   A  typical  procedure  is  to  determine  the  over- 
consolidation  ratio  by  standard  means  and  to  determine  the  ratio 
K  for  the  particular  type  of  soil  under  consideration  through 
a  load-release  cycle  to  achieve  the  same  overconsolidation 
ratio  as  exists  in-situ. 

Existing  field  and  laboratory  methods  are  limited. 
In-situ  methods  measure  only  single  components  of  the  stress 
tensor  or  a  stress  ratio.   Laboratory  techniques  are  too 
dependent  on  information  that  may  not  be  available,  and  the 
assumption  that  the  geologic  history  of  a  site  may  be 
reproduced  on  a  laboratory  time  scale. 

New  methods  and  hardware  capable  of  measuring  the 
complete  in-situ  stress  tensor  are  required  to  take  full 
advantage  of  modern  analysis  tools.   The  studies  reported 
herein  were  designed  as  a  first  definitive  step  in  the 


fulfillment  of  these  needs.   The  scope  includes  evaluation 
of  existing  stress  measurement  methods  and  the  development 
of  improved  concepts  for  determining  the  full  stress  tensor 
The  specific  objectives  were  to: 

(A)  Summarize  current  methods  applicable  to 
the  determination  of  the  state  of  stress 
in  soils.   Included  are  in-situ  field 
determination  of  in-situ  normal  and  shear 
stresses,  and  the  principal  stress  directions, 
as  well  as  field  and  lab  methods  for  the 
estimation  of  stress  history.   This  is 
discussed  in  Section  2  of  this  report. 

(B)  Study  the  feasibility  of  applying  fabric 
analysis  in  soils  for  determining  the  in-situ 
stresses.   A  number  of  subordinate  objectives 
were  included  to  define  the  extent  to  which 
fabric  analysis  can  be  applied  to  soils  and 
the  results  which  may  be  expected.   This  is 
discussed  in  Section  3  of  this  report. 

(C)  Evaluation  of  the  feasibility  of  utilizing 
anisotropic  properties  of  undisturbed  soil 
for  determining  in-situ  stresses.   The 
following  technologies  were  reviewed  for 
this  purpose. 

o  Soil  mechanics 

o  Water  permeability 

o  Sound  conduction 

o  X-rays 

o  Gamma -rays 

o  Electric  conductivity 

o  Electron  beam 

o  Electromagnetic  resonance 

o  Infrared 

o  Electron  microscopy 

These  are  discussed  in  Section  2.5  and  Section  4. 

(D)  Bring  to  the  level  of  a  conceptual  design 
the  novel  technique  most  promising  for  the 
determination  of  in-situ  stresses.   Included 
are  theoretical  analysis,  development, 
evaluation  and  the  delivery  of  preliminary 
drawings  together  with  a  work  plan  to 
realize  the  projected  hardware.   This  is 
discussed  in  Section  5. 


Because  all  users  of  this  report  will  not  be 
geotechnical  engineers,  brief  reviews  of  soil  stress  notation 
and  a  glossary  of  fabric  analysis  terms  are  appended. 


2.   EXISTING  TECHNOLOGY 

2.1  STATE  OF  STRESS 

2.1.1  SOIL  STRESSES 

The  three  dimensional  state  of  stress  at 
a  point  may  be  expressed  as  a  second-order  symmetric 
tensor  by 


iJ 


Gll 

a12 

a  21 

a22 

a31 

a32 

°13 
a  23 

a33 


where  a.,  are  normal  stresses  on  three  mutually 


perpendicular  planes  and  aAT>  =  a^  are  the  corresponding 
shear  stresses  on  the  same  planes.   If  these  planes  are 
normal  to  the  principal  stress  directions,  the  tensor 
becomes 


a.  . 


in  which  the  a 
stresses  a 


Gll 

0 

0 

0 

a22 

0 

0 

0 

a33 

are  the  principal  stresses,  and  the  shear 


'AB 


AA 
are  zero. 


Several  options  are  available  in  the  determination 
of  the  stress  tensor.   The  shear  and  normal  stresses  on 
three  arbitrary  orthogonal  planes  may  be  determined  or 
the  principal  stresses  and  their  directions  may  be 
determined.   The  choice  in  any  given  case  depends  on  the 
technique  used  to  make  the  measurement.   Because  direct 
measurements  of  shear  stress  in-situ  are  usually  difficult, 
direct  measurement  techniques  are  generally  used  to  monitor 
normal  stress.   Certain  indirect  techniques,  e.g.,  fabric 
analysis,  deal  with  shear  stresses  more  effectively. 


Traditionally,  the  intermediate  principal  stress 
has  not  been  considered  important  in  geotechnical  appli- 
cations.  Its  role  is  now  the  subject  of  an  increasing 
amount  of  attention  (Shibata  and  Karube,  1965;  Broms  and 
Casbarian,  1965;  Saada  and  Zamani  1969)  and  it  is  important 
that  methods  of  measuring  the  full  stress  tensor  in-situ  be 
developed  to  make  use  of  these  and  other  efforts. 

The  geologic  events  that  lead  to  anisotropy  of 
in-situ  stresses  are  numerous.   Any  event  that  changes  the 
conditions  of  lateral  constraint,  such  as  river  valley  erosion, 
tectonic  movements,  or  nonlevel  ground  surface  may  be 
expected  to  result  in  modification  of  the  lateral  stresses. 
For  this  reason, a  knowledge  of  the  geologic  history  alone  is 
not  sufficient  to  estimate  in-situ  stresses. 

2.1.2   EARTH  PRESSURE  AT  REST 

The  coefficient  of  earth  pressure  at  rest,  K  ,  is 

o 

defined  as  the  ratio  of  horizontal  to  vertical  effective 

stresses,  for  the  case  where  the  principal  directions  are 

vertical  and  horizontal  and  the  lateral  strain  is  zero. 

For  example,  "at  rest"  conditions  may  be  obtained  in  a  tri- 

axial  test  with  the  appropriate  stress  ratio,  or  in  a 

natural  deposit  constrained  against  lateral  deflection.   To 

the  extent  that  sidewall  friction  exists,  K  conditions  are 

'   o 

not  present  in  a  consolidation  apparatus  or  a  silo.   An 
infinite  slope  or  inclined  strata  do  not,  of  course,  impose 
restraint  conditions  for  zero  lateral  strain. 

The  range  of  values  that  may  be  assumed  by  K  is 
limited  by  the  strength  of  the  soil  to 


a  —    o  —   p 

where  K  and  K  are  the  active  and  the  passive  earth  pressure 
a      p 

coefficients  that  define  the  stress  ratios  at  failure.   For 
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a  soil  having  cohesion  c ',  friction  0' ,  and  density  7, 
these  active  and  passive  limits  are 

Ka  =   tan   (45  -  j)      yT~   tan  ^45  ~I~) 

K   =   tan2  (45  +  |)   +  yf^-  tan  (45  +  |^-  ) 

at  depth  z.  The  term  yz  replaces  the  major  principal  stress 
a,  in  the  cohesion  terms  because  K  is  defined  only  when  the 
principal  directions  are  vertical  and  horizontal. 

Alternatively,  K  may  be  expressed  in  terms  of 
Poisson's  ratio  v  as 

K 


Although  Poisson's  ratio  is  frequently  assumed  to  lie 
between  zero  and  one-half,  soil  is  not  an  elastic  material 

and  K  values  greater  than  one  often  occur, 
o        ° 

The  history  of  a  site  is  of  major  importance  in 
the  magnitude  of  K  .   For  a  history  of  monotonically  increasing 
vertical  stress,  as  in  a  normally  consolidated  soil,  the 
value  of  K  must  be  less  than  one  and  is  frequently  near  the 

value  K  =  1/2.   If  the  vertical  stresses  are  decreased  by 

o  J 

the  removal  of  overburden,  the  horizontal  stresses  tend  to 
remain  locked  in  at  relative  high  levels.   For  heavily  over- 
consolidated  soils,  K  may  approach  the  passive  failure 

condition  K  .   Such  a  stress  history  is  indicated  in 

P  J 

Figure  2-1  .   Alternatively,  the  horizontal  stresses  may  be 

released,  as  by  river  down  cutting,  so  that  the  active 

failure  condition  is  approached  as  indicated  in  Figure  2-2. 

Tectonic  stresses,  changes  in  ground  water  level  and  other 

events  may  also  modify  the  in-situ  stresses. 

The  behavior  of  a  soil  mass  is  strongly  dependent 
upon  the  existing  in-situ  stresses.   Bjerrum  (1967)  suggested 
that  progressive  slope  failure  is  encouraged  by  conditions 
of  high  lateral  stress  and  low  residual  strength.   Henkel 
(1970)  used  a  simple  stress  space  to  indicate  the  influence 


Horizontal  Stress  a. 


Figure  2-1.   Stress  History  for  K  Overconsolidated  Site 
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Horizontal  Stress  a, 


Figure  2-2.   Stress  History  for  Normally  Consolidated  Site 
with  Lateral  Stress  Relief 


of  initial  stresses  and  stress  path  and  soil  behavior.   His 
curves  are  combined  in  Figures  2-3  and  2-4.   Typical  stress 
paths  for  normally  consolidated  and  overconsolidated  soils 
are  shown  for  both  active  and  passive  conditions.   The 
differences  in  behavior  that  result  from  the  different  stress 
conditions  and  stress  paths  are  obvious.   Duncan  and  Dunlop 
(1969)  present  finite  element  simulations  of  excavation  in 
a  normally  consolidated  soil  (K  =  0.81)  and  in  a  heavily 
overconsolidated  soil  (K  =  1.60).   Predicted  shear  stresses 
were  found  to  depend  strongly  on  the  value  of  K  .   This  may 
be  expressed  by  the  factor  of  safety  necessary  for  the  0  =0 
analysis  if  shear  stresses  are  to  be  everywhere  less  than 
the  shear  strength.   The  equilibrium  of  the  normally  consoli- 
dated soil  (K  =  0.81)  required  a  safety  factor  of  3;  the 
equilibrium  of  the  overconsolidated  soil  (K  =  1.60)  required 

a  safety  factor  of  5.   The  K  values  were  held  constant  with 

J  o 

depth  in  this  analysis,  as  indicated  in  Figure  2-5.   A  highly 
overconsolidated  soil  may  approach  passive  failure  conditions, 
as  was  the  case  of  Bradwell  (Skempton,  1961).   In  such  a  case, 
K  would  decrease  with  depth  to  the  limiting  value  given  by: 

KQ  =  Kp  =  tan2  (45  +  \  )  +  ff  tan(45  +  f) 

and  the  stress  distribution  would  be  as  shown  in  Figure  2-6. 
At  Bradwell,  K  behaved  in  a  more  complex  fashion 
increasing  with  depth  to  a  maximum  value  then  decreasing. 
In  such  an  extreme  case,  the  analysis  of  Duncan  and  Dunlop 
would  doubtless  have  indicated  even  greater  differences 
between  the  equilibrium  requirements  of  the  normally  and 
overconsolidated  cases. 

2.2   EVALUATION  OF  THE  COEFFICIENT  OF  THE  EARTH  PRESSURE 
AT  REST 

If  the  stress  history  of  a  soil  mass  is  known  to  be 
normally  consolidated,  or  overconsolidated  without  lateral 
strain,  and  the  vertical  is  a  principal  stress  direction,  an 
estimate  of  K  may  be  obtained  from  the  soil  properties. 
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Overconsolidated 


ov,  =  1.6  a 
n        v 


Stress 


Figure  2-5.   Stress  Distribution  for  K  Constant  with 


o 


Depth  in  Normally  and  Overconsolidated  Soils 


Overconsolidated 


a,  =  1.6 
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a+2 
v    c 


Stress 


Figure  2-6 


Stress  Distributions  for  Normally  Consolidated 
(K  constant)  and  for  Overconsolidated (K  =K  ) 
°  Soils  °   P 
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2.2.1  THEORETICAL  EVALUATION,  NORMALLY 

CONSOLIDATED  K 
o 

The  approximation  of  Jaky's  expression  K  =  1  -  sin  0 
is  frequently  used  for  normally  consolidated  soils.   The 
development  of  this  expression  (Kezdi)  is  based  on  the  equi- 
librium of  a  wedge-shaped  mass  (Figure  2-7)  with  side  slopes 
at  an  inclination  0.   Two  zones,  ABO  and  CDO ,  exist  in  which 
the  shear  strength  is  fully  mobilized.   These  zones  must  be 
in  plastic  equilibrium,  with  the  state  of  stress  known: 

2 

az  =  T  1   sin  0 


sin  0  cos  0 

ax  =  t  ctn  0 

•  2  w 
T   =  T  =  yz  sin  0  cos  0  -  x  sin  0 

Between  the  two  plastic  zones  there  is  a  transition  zone  BOC 
in  which  the  shear  strength  is  not  fully  mobilized.   The 
stresses  within  this  wedge  are  computed  on  the  basis  of 

1)  the  continuity  of  shear  stresses  across  the 
^oundarys  BO  and  CO . 

2)  assumption  that  the  shear  stress  t   is 
distributed  parabolically  in  the   '  x  direction. 

3)  equilibrium  and  compatibility  conditions. 
On  the  surface  BO,  with  x-^  =  z   tan  (45°  -  0/2) 

the  stresses  are  given  by: 


CTz   =     Yxl 


1  +  sin2  0 


cos  0 


Gx  =  Yxl  cos  ^ 


t      =     x_Y    sin  0 
xz  1 
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The  stresses  in  BOC  are  not  known,  but  t   =  0  at  x  =  0  and 

'      xz 

t   =  xy  sin  0  on  OB  and  OC .   A  parabolic  distribution  is 
assumed: 

2/  2 

T     =  Tn  X  /X, 
XZ      11 

The  equilibrium  condition 

ho  £t 

— ?  +  — 2«  =    Y 

Sz       dx 

gives    the  condition  a      =  Yz  at  x  =    0.      From 

5gx  5Txz 

dx  Bz 

the  condition  a  =  Yz  (1-sin  0  )  1  +  2/3  sin0  on  x  =  o. 
x  1  +  sin0 


Jaky  and  his  followers  assumed  that  the  x  =  0  plane 

has  zero  lateral  strain  because  of  symmetry  conditions.   He 

formed  the  ratio  a  /a   which,  if  this  assumption  is  valid, 

is   the  ratio  K   : 
o 


tfx  1  +   (2/3)    sin0 

Ko  =     -£     =      (l-sin0)  \      '  

u         ay  1  +  sin  V 

Various  approximations  of  this  expression  are 

K  =  1  -  sin0  and  K  =  0.9  (l-sin0)  . 
o  o 

These  are  plotted  in  Figure  2-8.  The  simple  K  =  l-sin0 
is  the  form  most  frequently  employed  to  estimate  K  .   The 
plane  x  =  0  may  not  be  truly  in  the  at-rest  condition,  the 
zero  deflection  condition  at  x  =  0  being  the  result  of 
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symmetry  rather  than  necessarily  representing  a  state  of 
zero  lateral  strain.   The  scatter  in  reported  K  values  is 
greater  than  the  difference  in  the  two  expressions  and  the 
simpler  form  is  commonly  used. 

2.2.2  EMPIRICAL  EVALUATION,  OVERCONSOLIDATED  K 

o 

Considerable  data  have  been  published  relating 

lateral  stresses  to  overconsolidation  ratio  and  soil  index 

properties.   Typical  relationships  include: 


KQ  =  \   +  a  (Pr  -  1) 


Sherif  and  Strazer  (1973) 


and 


K 


(C   P*) 


Schmidt  (1973) 


in  which 


nc 
K   =  normally  consolidated  K 

o  J  o 


P   =  overconsolidation  ratio 
r 


X,   a,  h   =  empirical  constants. 

Sherif  and  Strazer  generate  the  empirical  parameters 
a  and  A  by  plotting  K  as  a  function  of  overconsolidation 
ratio  (Figure  2-9).  Values  of  a  and  \   for  various  soils  are 
plotted  in  Figure  2-10.   If  Figure  2-10  is  used  to  predict 
a  and  A,  errors  are  in  the  range  of  +  0.1  and  A  +  0.02  for 
a.   Then 

A  KQ  =  A  X  +  (Pr-1) A  a 

A  KQ  =  ±  0.1  +  .02  (Pr-1) 

and  for  very  high  overconsolidation  ratios,  the  variation 
in  a  may  become  important. 
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The  parameter  h  of  the  representation  suggested 
by  Schmidt  is  found  as  the  slope  of  a  log-log  plot  of  K 
against  overconsolidation  ratio.   Its  value  ranges  between 
0.3  and  0.6. 

The  K  vers.  P   data  of  several  investigators  are  plotted 
or  o  r 

in  Figure  2-11.   The  influence  of  soil  type  is  most  apparent 

at  large  overconsolidation  ratios.   Both  Skempton  (1961)  and 

Yudhbir  and  Prasad  (1973)  indicate  maximum  K  values  when 

v    '  o 

P  is  in  the  range  20  to  40;  in  Skempton 's  case  this  was 

attributed  to  surface  weathering  and  cracking  at  the  Bradwell 

site.   The  maximum  value  of  K  is  limited  by  the  passive 

earth  pressure  K  ,  indicating  that  the  continuously  upward 

trends  of  the  data  of  Sherif  and  Strazer  and  of  Brooker  and 

Ireland  must  at  least  level  off  at  high  P  values.   If  K 

°   r  o 

approaches  K  ,  vertical  swelling  of  the  soil  deposit  acts 
to  cause  unloading  along  the  K  yield  surface,  as  indicated 
in  the  typical  stress  trajectory  inset  in  Figure  2-11.   The 
Brooker  and  Ireland  (1965)  data  were  originally  shown  as 
functions  of  plasticity  index  (Figures  2-12  and  2-13) . 

It  can  be  seen  that  an  approximate  value  of  K  may 
be  determined  if  the  soil  mass  is  known  not  to  have  experi- 
enced lateral  deflection.   In  some  areas,  local  experience 
may  provide  guidance  if  knowledgeably  applied.   However, 
in-situ  lateral  stresses  can  be  modified  radically  by 
localized  conditions  and  more  direct  methods  of  estimating 
their  actual  value  are  needed. 

2.3   LABORATORY  DETERMINATION  OF  STRESSES 

A  number  of  techniques  are  available  to  estimate 
in-situ  stresses.   Some  measure  the  soil  behavior  under 
K  conditions  and  require  the  following  assumptions: 

o  Principal  stresses  are  vertical  and 

horizontal 
o  The  soil  mass  has  a  history  of  zero 

lateral  strain 
o  The  vertical  stress  history  is  known 
o  The  soil  response  over  geological 

periods  are  reflected  in  short-term 

laboratory  tests. 
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Others  measure  some  property  of  the  sample  that  reflects  the 
stresses  experienced  in-situ. 

2.3.1   SPECIAL  TRIAXIAL  TESTS 

The  value  of  K  may  be  estimated  by  laboratory 
tests  under  conditions  of  zero  lateral  deflection.   At  the 
1958  Brussels  Earth  Pressure  Conference,  Bishop  (1958) 
described  triaxial  tests  to  maintain  zero  lateral  deflection. 
The  basis  of  one  method  rests  on  the  equality  between 
vertical  and  volumetric  strains  when  lateral  strains  are 
zero.   This  is  applied  in  the  K  -test  cell  shown  schematically 
in  Figure  2-14.   As  long  as  the  chamber  and  confining  fluid 
may  be  considered  as  allowing  no  lateral  strain  in  the  soil, 
the  equality  of  vertical  and  volumetric  strains  in  the  soil 
is  enforced  by  the  fact  that  the  loading  piston  and  soil 
specimen  diameters  are  identical.   The  test  is  invalid  if 
pore  pressure  gradients  are  established  in  the  soil  or  if 
the  soil  is  not  homogeneous  since  either  condition  may  induce 
localized  nonzero  lateral  strains  in  the  sample.   Bishop 
stated  that  previous  investigators  were  unsuccessful  in 
applying  this  technique.   Bishop's  second  method  includes 
monitoring  the  lateral  strain  across  the  specimen  diameter 
and  manipulating  the  stress  ratio  to  maintain  zero  lateral 
strain  and  uses  a  number  of  devices  (Bishop  and  Henkel,  1961) . 
Other  workers  have  developed  sensitive  lateral  strain  sensors 
in  standard  triaxial  cells  to  maintain  K  loading  conditions. 
Another  method  mentioned  by  Bishop  was  the  measurement  of 
overconsolidation  ratios  in  both  the  horizontal  and  vertical 
directions  on  undisturbed  samples. 

Altschaeffl  and  Mishu(1970)  used  a  steel  ring  sur- 
rounding the  soil  specimen  as  a  noncontacting  radial  deflection 
capacitive  transducer.  The  capacitance,  C,  between  coaxial 
cylinders  of  radius  a  and  b  is  given  by 


C  = 


2tt£  Lk 

Ml" 
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where 

e   =   permittivity  constant  = 

8.85  x  10    Farads /meter 
k   =   dielectric  constant 
L   =   length  of  the  shorter  cylinder 

For  a  given  set-up,  the  capacitance  depends  on  the  radius, 
a,  of  the  soil  sample.   The  gage  was  driven  by  a  standard 
capacitance  bridge  as  shown  in  Figure  2-15.   An  alternative 
electronics  package  for  this  type  of  noncontacting  probe 
exists  in  the  form  of  proximity  gage  electronics.   Such  a 
system  provides  a  DC  voltage  analog  to  the  change  in 
capacitance,  without  requiring  continuous  manual  nulling 
and  phase  balancing.   The  capacitance  bridge  package  should 
perform  well  in  K  testing  since  the  intent  is  to  avoid 
any  changes  in  the  radius  of  the  soil  specimen.  With  proper 
control  of  the  stress  ratios,  the  capacity  of  the  ring-soil 
capacitor  should  not  change  except  with  changes  in  the  k  of 
the  fluid  with  temperature  or  pressure.   Geometries  other 
than  the  sample  in  a  ring  used  by  Altschaeffl  and  Mishu  are 
also  useful.   Proximity  gage  transducers  frequently  use 
the  capacitance  between  a  sharp  point  and  a  conductive  plane 
to  make  local  measurements.   In  these  cases,  a  mechanical 
micrometer  is  included  to  calibrate  the  system  directly. 

In  general,  capacitive  transducers  have  the 
following  advantages: 

o  a  large  sensing  area  to  average  minor 
point-to-point  variations 

o  sensitivity  and  range  that  is  easily 

scaled  by  the  probe  geometry.   (Extremely 
sensitive  transducers  may  be  designed  if 
a  small  operating  range  is  permitted.) 

o  zero  contact  force. 

Their  disadvantages  are  that  : 

o  careful  shielding  is  required  to  avoid 
extraneous  signals  resulting  from  the 
capacity  of  the  probe  with  respect  to 
its  environment,  and 

o   calibration  needs  may  be  extensive. 
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Non-contacting  Lateral  Strain  Sensor 
(Altschaeffl  and  Mishu,  1970) 
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The  optical  methods  suggested  by  Dudley  (1971) 
include  photogrammetry  and  optical  tooling  techniques 
to  monitor  the  sample  diameter  in  a  triaxial  cell.   The 
optical  tooling  technique  provides  real  time  measurement  of 
specimen  diameter,  as  needed  for  a  K  test.  Two  10-second 
theodolites  were  set  up  47  in.  (1.2  m)  from  the  specimen  center, 
to  monitor  the  angles  subtended  by  two  specimen  diameters. 
This  system  is  said  to  determine  the  specimen  dimension  to 
an  accuracy  of  0.005  in.  (13x10  m)  .   This  is  equivalent  to 
an  uncertainty  of  about  0.27o  strain  for  the  2.42-in. 
(61.5x10  m)  diameter  specimens  used.   Dudley  suggests  that  the 
accuracy  might  be  increased  to  0.001  in.  (3x10  m)  by 
improved  methods.   The  triaxial  cell  wall  acts  as  an  optical 
lens,  increasing  sensitivity. 

C.A.  Moore  (1971)  used  foil  strain  gages  mounted 

on  a  foil  band  to  monitor  radial  strain  in  K  tests  on  soils 

o 

containing  mica.   The  sensor  (Figure  2-16)  is  simply  a  gaged 
foil  band  wrapped  around  the  specimen,  tightened  by  an 
attached  tubing  clamp.   The  pin  and  hole  arrangement  consti- 
tutes a  remotely  operated  release  to  permit  removal  of  the 
transducer  from  the  specimen  before  soil  failure.  Moore 
used  his  lateral  strain  sensor  to  provide  a  feedback  signal 
in  a  closed  loop  stress  control  system  for  his  triaxial 
apparatus.   The  confining  pressure  was  increased  at  a  constant 
rate  while  the  axial  stress  was  automatically  controlled  for 
zero  lateral  strain.   The  closed  loop  system  was  added  to 
avoid  the  cummulative  effect  of  incremental  loading  steps 
and  re -adjustments  that  may  be  expected  in  a  less  sophisticated 
system. 

The  problem  of  accumulating  errors  as  a  zero  deflec- 
tion control  system  hunts  for  the  correct  stress  ratio  was 
approached  differently  by  Andrawes  and  El-Sohby  (1972, 1973) . 
These  authors  conduct  several  constant  stress  ratio  tests  at 
slightly  different  stress  ratios.  The  stress  ratio  at  which 
axial  and  volumetric  strains  are  equal  is  interpolated  as 
shown  in  Figures  2-17.   Figure  2- 17 (a)  indicates  that  the 
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value  of  K  for  the  dense  sand  is  between  1/2.3  and  1/3.0. 
Interpolation  to  zero-strain  ratio  in  Figure  2-17 (b)  yields 
the  value  of  K  =  1/2.9  =  .34.  Accuracy  as  great  as 
justified  by  the  apparatus  and  instrumentation  being  used 
may  be  obtained  by  additional  triaxial  tests  at  stress  ratios 
closer  to  the  interpolated  value  of  1/K .  Andrawes  and 
El-Sohby  (1973)  measured  axial  strain  by  dial  indicator  and 
volumetric  strain  by  burette. 

Utilization  of  K  triaxial  tests  to  estimate  in-situ 
stresses  requires  that  the  soil  mass  be  consolidated  under 
K  conditions,  that  the  major  principal  stress  be  known,  and 
that  the  triaxial  test  adequately  simulate  the  stress-strain 
history  of  the  soil  mass. 

2.3.2  KQ  DEVICES 

Special  test  devices  have  been  designed  to  enforce 

conditions  of  zero  lateral  strain  (Brooker  and  Ireland, 

1965;  Obrcian,  1969;  and  Sherif  and  Strazer,  1973).   These 

devices  are  characterized  by  rigid  sidewalls,  instrumented 

for  lateral  stress  but  permitting  little  or  no  lateral 

strain.  Unlike  K  triaxial  tests  in  which  variations  in 

o 

lateral  strain  may  occur  which  average  to  zero  over  the 
measured  portion  of  the  specimen  ,  the  rigid  sidewalls 
force  the  lateral  strain  to  remain  zero.  This  advantage 
is  countered  by  the  disadvantage  that  the  rigid  sidewalls 
cannot  be  made  perfectly  frictionless,  thus  may  induce 
shear  stress  on  what  should  be  a  principal  stress  plane. 
For  this  reason,  the  specimen  height  is  made  small  relative 
to  the  diameter,  as  is  done  in  one -dimensional  consolido- 
meters  for  the  same  reason. 

The  apparatus  described  by  Brooker  and  Ireland  (1965) 
employs  a  steel  membrane  to  restrict  lateral  strain  as  shown 
in  Figure  2-18.   Foil  strain  gages  on  the  rear  of  the 
membrane  monitor  circumferential  strain,  which  is  maintained 
within  about  10  microstrains  by  oil  pressure  on  the  back. 
Typical  data  obtained  are  indicated  in  Figures  2-19  and  2-20. 
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The  stress-time  curve  in  Figure  2-19  indicates  the  pore 
pressure  decay  at  each  loading  step.   The  stress  path  in 
Figure  2-20  was  constructed  using  the  stresses  at  100% 
consolidation.   This  figure  also  shows  that  initial  loading 
K  is  constant  for  this  material  to  at  least  2200  psi 

(15  x  10  Newton/meter  ). 

Obrcian  (1969)  employed  a  K  device  (Figure  2-21) 
having  a  heavy  wall  split  confining  ring  instrumented  to 
monitor  lateral  stress  and  side  wall  friction.   These  ^ata 
indicate  that  mean  shear  stresses  induced  on  the  confining 
ring  were  on  the  order  of  10  percent  of  the  major  principal 
stress,  that  is,  roughly  1/4  to  1/3  the  magnitude  of  the 
minor  principal  stress  at  the  K  values  reported.   Typical 
data  shown  in  Figure  2-22  appear  to  indicate  a  decrease  in 
K  values  with  increasing  stress.   Obrcian  points  out  that 
this  behavior  is  consistent  with  increasing  density  at 
higher  stress  levels,  yet  this  appears  to  conflict  with 
investigators  who  find  that  K  is  constant  for  virgin 
loadings.   If  the  data  in  Figure  2-22  are  used  to  back 
calculate  the  lateral  stresses,  the  relationships  in 
Figure  2-23  emerge,  indicating  that  the  hyberbolic  relation- 
ship between  K  and  a-,  may  be  the  result  of  a  seating  error 

°  4    2 

of  approximately  1.5  psi  (10  N/m  )  in  the  lateral  stress. 

Thus  there  appears  to  be  no  significant  discrepancy  between 

the  Obrcian' s  data  and  the  findings  of  other  investigators 

that  on  initial  or  virgin  loading  KQ  is  constant. 

Sherif  and  Strazer  (1973)  employed  the  University 
of  Washington  Stress-Meter  (Figure  2-24)  in  their  work. 
Typical  data  obtained  with  this  K  device  are  shown  in 
Figure  2-25. 

2.3.3  LABORATORY  TESTS  FOR  IN -SITU  STRESSES 

At  least  two  techniques  that  use  laboratory  tests 
to  evaluate  the  stresses  to  which  the  soil  sample  is 
subjected  in  situ  are  available. 
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The  lateral  consolidation  test  mentioned  by  Bishop 
(1958)  appears  to  have  been  developed  by  Poulos  and  Davis 
(1970,  1972).   The  procedure  is  to  consolidate  an  undis- 
turbed sample  in  the  horizontal  direction  while  holding  the 
vertical  stress  at  its  assumed  in-situ  value. 

o  A  triaxial  specimen  is  loaded  to  the 
estimated  in-situ  vertical  stress 
value  using  a  low  stress  ratio.   The 
confining  pressure  should  be  sufficiently 
large  to  just  prevent  failure,  and  must 
be  less  than  the  anticipated  value  of 

V 

o  The  confining  pressure  is  gradually 
increased  while  the  axial  stress  is 
held  constant  at  the  in-situ  vertical 
stress  level. 

o  The  observed  void  ratio  is  plotted 
against  confining  pressure.  The  plot 
should  resemble  an  e-log  P  consoli- 
dation curve  and  display  a  break  near 
the  in-situ  lateral  stress. 

o  The  sequence  may  be  repeated  using 
specimens  of  varying  disturbance.   If 
the  estimated  lateral  stress  is  plotted 
against  disturbance  ratio,  the  results 
may  be  extrapolated  back  to  zero 
disturbance. 

Degree  of  disturbance   is   defined  as  Ae/    (Ae) 
to  v      'max 

where  Ae  is  the  difference  between  the  field  void  ratio 

and  that  of  the  laboratory  specimen  at  the  preconsolidation 

pressure  and  (Ae)    is  the  difference  between  the  field 
r  x  'max 

void  ratio  and  that  of  a  remolded  specimen  at  the  preconsolida- 
tion pressure.   Typical  data  on  laboratory-consolidated 
Kaolin  are  shown  in  Figures  2-26  through  2-28.   Even  samples 
remolded  at  constant  moisture  content  (disturbance  ratio  =  1.0) 
retain  a  "memory"  of  the  in-situ  stress.   Note  that  the  test 
series  consolidatd  at  K0=  1.0  displays  a  less  distinct  break 

relative  to  the  K  -consolidated  test  series.   The  method  cannot 

o 

be  used  for  heavily  overconsolidated  clays. 
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A  rather  indirect  method  for  dealing  with  over- 
consolidated  clays  was  employed  by  Skempton  (1961)  at  the 

Bradwell  site.   In  this  case,  the  value  of  K  was  determined 

o 

by  an  indirect  stress -relief  method  employing  the  transfer 
of  lateral  effective  stress  to  neutral  stress  that  occurs 
during  the  sampling  procedure  itself.   Briefly,  when  a 
sample  is  removed  from  the  soil  mass,  the  total  stress  in 
the  sample  is  reduced  to  zero.  The  relative  compressibilities 
of  the  soil  matrix  and  the  pore  water  result  in  a  transfer  of 
effective  stress  from  the  soil  to  tension  in  the  pore  water. 

That  is,  when  the  sample  is  removed  from  the  ground, 
the  total  stress  on  the  sample  is  reduced  from  the  in-situ 
value  to  zero.  The  soil  matrix  expands  slightly,  inducing 
tension  in  the  pore  water.  Thus  we  have  a  change  in  stress 
(caused  by  removing  the  sample  from  the  ground)  that  is 
related  to  a  change  in  pore  pressure,  and  can  make  use  of 
Skempton' s  pore  pressure  coefficient  A.  The  total  stresses 
in-situ  are: 

a  -  a'  +  u 
v     v    o 

a,  -  (a-u)K  +  u  =  a1  K  +u 
h    v  v  o'   o    o    vo    o 

in  which  u  is  the  pore  pressure  at  depth  Z.   After  removal 

from  the  ground,  the  effective  stress  in  the  sample  is  -u 
where  u  is  the  pore  pressure  in  the  sample.   The  sample  pore 
pressure  is  given  by 


u  ■  u  +  Au 
o 


u  -  u  +  1/3  (Aa1  +  2Aa3)  +  (A  -  1/3) 

The  change  in  stress  due  to  the  release  of  the 
in-situ  stress  is: 

A  1       v       o    v 


Act-  -  Aa~ 


A°3  =  _0h  =  -°v  Ko  "  uo 
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leading  to 

<u/ctv  >  +  A    for  K  >  1 

K.   -   O 

°      A  -  1 

-(u/a')+(A-2/3) 

v for  K  <  1 

Ko  '     A  +  1/3        ° 

in  which  A  is  the  appropriate  pore  pressure  coefficient  and 
the  other  symbols  take  their  accepted  meanings.   An  extensive 
series  of  tests  was  conducted  at  the  Bradwell  site  to 
determine  the  neutral  stress  both  in-situ  and  in  samples.   The 

resulting  K  values  are  shown  in  Figure  2-29.   The  decrease 
in  K  at  depths  less  than  20  ft  (6m)  is  attributed  to  the 
reduction  of  lateral  stress  resulting  from  weathering  and 
softening  of  the  London  clay  prior  to  the  disposition  of 
the  post-glacial  Marsh  clay.   At  depths  near  30  ft  (9m)  , 
K  is  near  the  approximate  value  of  K  ,  indicating  a 
possible  failure  zone.   The  method  employed  by  Skempton 
at  the  Bradwell  site  was  in  part  a  solution  in  response 
to  a  specific  problem.   Although  it  has  since  been  applied 
by  other  investigators,  e.g.,  Blight  (1969),  it  may  be 
somewhat  elaborate  for  routine  site  investigation  programs. 
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2.4   IN-SITU  DETERMINATION  OF  STRESSES 

Several  techniques  that  may  be  used  in-situ  for 
stress  measurement  are  discussed  below.   These  devices 
measure,  in-situ,  some  soil  property  or  effect  that  may  be 
related  to  one  or  more  components  of  the  in-situ  stress 
field.   Pre-knowledge  of  the  stress  history  or  over- 
consolidation  ratio  is  not  generally  required  in  a  direct 
in-situ  technique. 

2.4.1  ACOUSTIC  WAVE  PROPAGATION 

The  words  acoustic  and  sonic  are  taken  as  descrip- 
tive of  low  amplitude  waves  which  may  be  of  any  convenient 
frequency.   The  propagation  velocity  of  such  waves  in  a 
continuum  is  given  by: 

2 


-i  (1  + 


l-vj  -2v 


d   -» 


1/2 


for  the  dilatational  wave  velocity  and  by 

-.1/2 


V. 


r"   [2  a+o] 


d/JJ 


for  the  distortional  wave  velocity,  in  which 

Ei  =  dynamic  elastic  modulus 

vd  =  dynamic  Poisson's  ratio 

Gi  =  dynamic  shear  modulus 

p  =  mass  density 

Acoustic  velocity  measurements  are  frequently  used 
to  determine  dynamic  elastic  moduli  (e.g.,  Afifi  and 
Richart,  1973;  Bamert ,  et  al. ,  1965).   Because  the  elastic 
moduli  of  rock  and  soil  are  typically  stress-dependent, 
the  wave  propagation  velocities  may  be  correlated  directly 
with  in-situ  stress.   Johnson  and  Nelson  (1970)  report  a 
seismic  field  study  that  included  determination  of  elastic 
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moduli  in  several  directions  at  a  rock  site.   The  velocity- 
data  and  a  corresponding  in-situ  velocity  ellipse  are 
shown  in  Figure  2-30.   It  is  assumed  that  the  velocity 
ellipse  corresponds  to  the  in-situ  stress  ellipse. 

Hardin  and  Richart  (1963)  conducted  resonant  column 
tests  on  Ottawa  sand,  crushed  quartz  sand  and  crushed  quartz 
silt  to  determine  the  parameters  influencing  wave  velocities. 
They  proposed  correlations  having  the  form 

V   =   (A  -  Be)   a  C 
s     v      '        o 

in  which    V   =   shear  wave  velocity, 

s  J 

e   =  void  ratio 

a   =  confining  pressure 

A,B  &  C   =  empirical  constants. 

For  dry  sands,  gradation  and  relative  density  had  no  effect 
except  in  so  far  as  the  void  ratio  changed.   Crushed  quartz 
silt  exhibited  a  dependency  on  time  when  measurements  were 
made  one  hour  after  loading.   Substituting  the  empirical 
constants,   the  shear  wave  velocities,  V  ,  were  found  to  be: 

For  dry  Ottawa  sand: 

V  =   (170-78.2  e)  aj;/4         for      (<*>  2000  psf) 
so  o 

V  =   (119-56.0  e)  a  3/1°       for      (a«*  2000  psf) 

S  O  u 

For  dry  crushed  quartz  sand: 

Vs  =   (159-53.5  e)  a*/4 

For  dry  crushed  quartz  silt  (not  overconsolidated)  : 
Vc  =   (167-87.8  e)  Ql/4 
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Typical  data  for  Ottawa  sand  are  shown  in  Figure  2-31. 
The  velocity  data  are  predicted  to  within  about  +  10%  by  the 
correlations  noted.   If  the  correlations  are  reversed, 
i.e.,  measured  velocity  data  used  to  predict  stresses,  the 
associated  errors  would  be  on  the  order  of  +  40%  since 
"&c /a   =  4  BV/Vo   For  overconsolidated  quartz  silt,  there  is 
little  variation  of  wave  velocity  with  confining  pressure, 
so  that  the  correlation  given  holds  only  for  normally 
consolidated  quartz  silt. 

The  relationships  of  Hardin  and  Richart  were  applied 
by  Stokoe  and  Woods  (1972)  in  the  cross-borehole  method 
shown  in  Figure  2-32.   These  investigators  measured  shear 
wave  velocity  between  boreholes.   As  pointed  out  by  Haupt 
(1973)  in  a  discussion  of  this  paper,  refraction  of  the 
wave  path  as  indicated  in  Figure  2-33  makes  necessary  correc- 
tions when  measurements  at  shallow  depths  are  attempted. 
The  required  corrections  to  the  Stokoe  and  Woods  data  are 
shown  in  Figure  2-33 (b) . 

The  one -fourth- root  dependence  of  acoustic  velocity 
is  an  unattractive  feature  of  the  concept  because  the 
experimental  error  in  determination  of  velocity  is  reflected 
in  an  error  four  times  as  large  in  the  inferred  stresses, 
i.e. ,  5  percent  error  in  the  velocity  measurement  results 
in  a  20  percent  error  in  the  reported  stress;  A<j/q  =  4(AV  /V  ) 
Compounding  this  is  the  number  of  other  parameters  that 
affect  the  acoustic  velocity.   Circumvention  of  this 
problem  would  require  development  of  an  in-situ  calibration 
capability.   This  could  be  accomplished  by  direct  measurement 
of  in-situ  stress  in  the  formation  being  mapped,  or  by 
pressurizing  a  section  of  the  borehole  wall  in  which  a  local- 
ized velocity  versus,  stress  calibration  may  be  computed.   The 
major  advantages  of  the  cross  borehole  concept  are  its 
rapidity,  and  the  fact  that  measurements  are  conducted  on  the 
undisturbed  soil  far  from  the  boreholes.   In  all,  it  appears 
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Variation  of  Shear  Wave  Velocity  with  Void 
Ratio  for  Various  Confining  Pressures,  Grain 
Sizes,  and  Gradations  in  Dry  Ottawa  Sand 
(Series  I)  (Richart  &  Hardin,  1963) 
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to  be  a  good  technique  for  obtaining  a  generalized  view  of 
the  gross  stress  field,  but  will  not  give  results  of  great 
precision  in  any  localized  volume.   As  such,  it  would  be 
especially  valuable  when  used  in  conjunction  with  another 
method  capable  of  detailed  and  precise  measurements  in 
small  volumes. 
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2.4.2  BOREHOLE  DEFORMATION  TECHNIQUES 

A  number  of  techniques  based  on  deformation 
measurements  in  boreholes  or  shafts  have  been  suggested 
or  used.   In  rock  mechanics,  these  include  overcoring, 
(Obert  and  Duvall,  1967,  and  Jaeger  and  Cook,  1964)  and 
flat  or  curved  jacks ,  (Glossop,  1968).   It  is  difficult 
to  visualize  overcoring  or  stress  relief  methods  being 
applied  to  soil  although  it  should  not  be  immediately 
discarded.   The  major  difficulty  results  from  the  changes 
in  the  stress  field  near  a  borehole,  particularly  in  the 
case  of  yielding.   If  rock  yields  during  coring  (or  over- 
coring)  the  phenomenon  of  core  discing  occurs,  but  at 
least  an  estimate  of  the  in-situ  stresses  may  be  recovered 
(Obert  and  Duvall,  1967).   Applications  of  borehole 

deformation  techniques  in  soils  include  the  jacking  of  plates 
or  semi-cylindrical  sleeves   (Glossop,  1968;  Sherif  and 
Strazer,  1973) .   Many  investigators  have  reported  on  the 
development  or  application  of  various  types  of  pressure- 
meters  . 

The  reported  jacking  tests  (Glossop;  Sherif  and 
Strazer)  were  directed  toward  the  measurement  of  in-situ 
deformation  moduli.   However,  since  the  opening  of  a  bore- 
hole or  a  shaft  implies  stress  relief  of  the  radial  component 
at  the  wall,  the  initial  phase  of  a  jacking  test,  if 
properly  conducted,  results  in  an  estimate  for  the 
in-situ  lateral  stress,  much  as  the  odometer  test  provides 
a  value  for  the  vertical  consolidation  stress.   This  would 
be  subject  to  the  same  constraints  as  a  consolidation  test, 
which  measures  the  maximum  historical  applied  stress. 

Such  tests  have  been  described  by  Obert  and  Duvall  (1967) , 
Glossop  (1968),  Wroth  and  Hughes  (1972,  1973),  Menard  (1967), 
Shepherd  and  Charleson  (1971),  Baguelin,  et  al.,  (1972), 
Higgins  (1969)  ,  Meigh  and  Greenland  (1965)  ,  Smith  and  Smith 
(1968) ,  Gibson  and  Anderson  (1961) ,  and  Campbell  and  Hudson 
(1969) .   Although  earlier  work  was  done  in  Germany,  the 
development  of  the  pressuremeter  appears  to  rest  with  Menard. 
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Recently  major  improvements  were  made  at  Cambridge  by  Wroth 
and  Hughes  and  in  France  by  Baguelin  and  Jezequel.   The 
classical  Menard  pressuremeter  (Figure  2-34)  is  simply  an 
inflatable  probe  which  is  lowered  into  an  open  borehole 
and  expanded  against  the  soil. 

The  data  are  typically  a  plot  of  probe  volume  or 
diameter  as  a  function  of  pressure  and,  usually,  time.   The 
probe  pressure  is  increased  in  increments  and  the  volume 
measured  at  specified  intervals  after  the  application  of 
each  increment.   Idealized  results  are  indicated  in 
Figure  2-35.   The  "creep11  curve  is  the  volume  change 
recorded  between  readings  taken  30  and  60  seconds  after 
each  pressure  increment.   The  P-V  curve  represents  the  60- 
second  volume  readings  and  the  "inertia"  curve  represents 
the  pressure -volume  relationship  for  the  pressuremeter 
itself.   The  in-situ  stress  may  either  be  estimated  directly 
from  the  curve  or  it  may  be  evaluated  from  theoretical 
cavity  expansion  analysis.   The  pressure -volume  curve  may 
provide  an  estimate,  since  most  authors  make  a  statement 
such  as  "the  initial  sloping  portion  of  the  curve  results 
from  the  restoration  of  the  earth  pressure  to  the  undisturbed 
condition  prior  to  the  removal  of  the  soil  from  the  hole," 
(Campbell  and  Hudson,  1969).   Also,  the  first  knee  on  the 
P-V  curve  appears  to  correspond  roughly  to  an  abrupt  de- 
crease in  the  creep  values.   Such  an  abrupt  increase  in 
modulus  would  be  expected  at  the  level  of  the  lateral  in- 
situ  stresses.   Seldom,  however,  do  data  as  clean  as  these 
idealized  curves  result  from  actual  testing.   This  is 
partially  due  to  the  fact  that  most  investigators  are  inter- 
ested in  modulus  and  ultimate  strength  values,  so  that  the 
initial  portions  of  the  P-V  curve  may  be  completely  missed. 

The  classical  approach  to  the  stress  problem  is  to 
deal  with  the  ultimate  pressure  at  which  indefinitely  large 
deflections  occur.   This  is  dealt  with  by  theoretical 
analysis  of  the  expanding  cavity  problem  (Menard,  1957; 
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Gibson  and  Anderson,  1961;  Wroth  and  Hughes,  1972;  and 
Vesic,  1972),  which  relates  the  ultimate  pressure  to  the 
in-situ  stress „   Assumptions  regarding  the  soil  behavior 
are  made,  and  mathematical  models  are  developed  to 
predict  the  pressure.   The  result  is  a  relationship  between 
the  limit  pressure  and  the  soil  properties  and  volumetric 
behavior,  and  other  parameters  considered  appropriate  by 
the  particular  investigator.   The  problem  is  made  more 
difficult  by  the  post-yield  behavior  of  the  soil  near  the 
pressuremeter.   A  plastic  zone  is  formed  early  in  the  test 
which  modifies  the  stress  distribution  around  the  pressure- 
meter.   The  radius  of  the  plastic  zone  is  considered  by 
Vesic  (1972) .   Such  analysis  requires  detailed  knowledge 
of  the  post-yield  soil  behavior,  which  may  be  lacking. 
This  is  a  problem  that  requires  a  certain  amount  of  experience 
and  judgment  if  confidence  is  to  be  placed  in  the  results. 

Many  of  the  tests  conducted  with  the  pressuremeter 
in  this  country  were  directed  toward  the  evaluation  of 
specific  hardware.   Thus,  many  of  the  disadvantages  cited 
by  some  authors  are  not  inherent  in  the  concept,  but 
reflect  hardware  developmental  problems.   This  may  be  a  factor 
in  the  lack  of  general  acceptance  of  this  device.   Hardware 
development  is,  in  fact,  continuing,  as  demonstrated  by  the 
efforts  of  the  Cambridge  group  under  Wroth  and  Hughes,  and 
the  French  group  under  Baguelin  and  Jezequel. 

Both  devices  avoid  the  problem  of  stress  relief  and 
disturbance  by  drilling  into  the  soil  without  an  open  bore- 
hole.  The  Cambridge  probe,  shown  schematically  in  Figure  2-36, 
resembles  a  vertical  mole.   The  auger  level  should  be  main- 
tained such  that  there  is  no  tendency  for  the  soil  to 
yield,  either  into  or  ahead  of  the  advancing  probe.   The 
thin  membrane  is  held  flat  against  the  probe  body  by  a 
slight  vacuum  during  insertion  of  the  probe  into  the  soil. 
Note  that  the  membrane  is  approximately  the  same  thickness  as 
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a  standard  triaxial  membrane,  about  0.012  inch  (3x10  m) 
This  is  much  thinner  than  the  membrane  used  for  the  Menard 
borehole  pressuremeter  and  can  be  so  because  the  membrane 
is  everywhere  supported  by  the  soil  since  it  is  not  in  an  over- 
size borehole.   An  additional  difference  is  that  the  membrane 
is  inflated  by  gas  rather  than  water.  This  allows  use  at 
greater  depths  since  the  pressure  due  to  a  high  head  of  water 
in  the  drill  string  is  avoided.   Deflection  is  measured  by 
feeler  gages  acting  on  the  interior  of  the  bladder.   The 
Cambridge  device  may  be  used  in  conjunction  with  a  companion 
direct  lateral  stress  measurement  device  which  is  said  to 
be  more  accurate  than  the  pressuremeter  for  stress  measure- 
ments.  However,  when  using  the  pressuremeter  alone,  the 
in-situ  lateral  stress  should  be  the  level  at  which  expansion 
is  initiated. 

The  French  device,  shown  in  Figure  2-37,  differs 
from  the  Cambridge  device  in  the  following  respects.   The 
membrane  is  supported  by  a  thin  water  layer  during  insertion. 
The  parameter  monitored  during  a  test  is  the  volumetric 
expansion  of  the  probe  rather  than  the  measurement  of  radial 
deflection.   The  pressure  behind  the  probe  is  monitored 
during  insertion  and  can  provide  certain  indications  if 
the  slurry  pressure  is  so  high  as  to  force  mud  up  around 
the  outside  of  the  probe.   Should  this  occur,  the  soil  is 
already  disturbed  and  a  test  at  the  same  location  would  be 
useless.   Because  the  membrane  is  fluid-backed,  any  inequal- 
ity between  the  lateral  stresses  tends  to  redistribute 
during  insertion. 

During  insertion,  disturbance  to  the  in-situ  stresses 
may  be  caused  by  interface  shear  stresses  induced  at  the 
surface  of  the  probe.   This  shear  stress  is  reduced  by  the 
use  of  polished  steel  instrument  surfaces.   Potyondy  (1961) 
reported  shear  resistance  for  soil-polished  steel  in  terms 
of  frictional  and  cohesive  parameters  fw  =  5/0  and  f  =  a/c 
in  which  6  is  the  interface  friction  angle:  a  is  the  interface 
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adhesion.   Thus,  the  Mohr  failure  curve  may  be  expressed 
in  terms  of  the  soil  shear  strength  parameters  0  and  c  by 
t  =  c  f  +  a   tan  (0  fw) .   Potyondy  recommends  the  following 
values  for  use  with  polished  steel: 

dry  dense  sand  0.54 

saturated  dense  sand  0.64 

dry  dense  cohesionless  silt  0.79 

saturated  dense  cohesionless  silt   0.68 

saturated  loose  cohesionless  silt   0.40 

cohesive  granular  soil  0.40 

clay  0.50      0.50 

Leonards  (1965)  reports  work  that  included  the  effect  of 
lubricants  on  soil-polished  steel  interface  shear  resistance 
in  terms  of  a  coefficient  of  friction  p,  =  t/g.   Since 
|i  =  tan  6 ,  these  data  can  be  cast  into  the  same  form  as 
that  of  Potyondy,  giving  the  following  results: 


20-30  sand 

60 

-80   sand 

f0 

f0 

0.47 

0.55 

0.46 

0.76 

0.34 

0.45 

polished  steel 
teflon  coated  steel 
graphite  coated  steel 

This  indicates  that  a  teflon  coating  may  actually  increase 
the  shear  resistance  of  a  steel-sand  interface,  whereas 
graphite  performs  better  as  a  lubricant.  Potyondy 's  data 
show  that  polished  steel  surfaces  introduce  shear  stresses 
into  the  soil  that  range  from  40  to  80  percent  of  the  soil 
shear  strength.   For  dense  silt  a  polished  steel  surface  is 
little  better  than  a  wholly  rough  surface. 

Pressuremeters  have  been  used  in  rock  and  granular 
soils,  but  the  Cambridge  and  the  French  devices  are  considered 
most  reliable  in  the  softer  clays. 
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2.4.3   STRESS  PROBE 

A   direct  approach  to  the  in-situ  stress  measurement 
problem  is  to  somehow  insert  a  stress  gage  into  the  ground 
with  so  little  disturbance  that  accurate  measurements  can 
be  made.   It  should  be  noted  that  the  problem  of  making 
stress  measurements,  even  in  artificial  specimens  where  the 
gage  may  be  placed  with  practically  zero  disturbance,  is  not 
itself  trivial.   Nonetheless  quite  effective  soil  stress 
gages  are  commercially  available. 

The  first  attempt  at  a  stress  probe  may  have  been 
undertaken  by  Kenney  (1967)  while  at  the  Norwegian 
Geotechnical  Institute.   The  device  was  an  open  bottom  tube 
built  up  from  sheet  piling  sections,  in  all  approximately 

10  meters  in  length.   Total  stress  and  pore  pressure  gages 
were  placed  flush  with  the  outer  surface  at  several  locations 
along  the  tube.   The  device  was  allowed  to  sink  into  the 
quick  clay  under  its  own  weight.   Although  the  clay  was 
observed  to  displace  into  the  tube,  certain  disturbance  along 

the  outside  was  inevitable  as  the  tube  entered  the  soil. 

Stress  measurements  were  made  over  a  long  period  as  the  clay 

relaxed  against  the  tube.   Kenney  pointed  this  out  and 

did  not  use  the  term  K  in  his  report.   This  device  was 

o 

intended  to  resolve  a  discrepancy  between  reported 
laboratory  values  of  K  and  reported  shear  strength  para- 
meters for  Norwegian  normally  consolidated  quick  clays. 
In  fact,  in-situ  lateral  stresses  were  shown  to  be  greater 
than  laboratory  K  tests  had  indicated,  so  that  the  stress 
probe  successfully  dealt  with  the  specific  question  it  was 
designed  to  answer. 

The  current  extension  of  the  stress  probe  is  embodied 
in  the  device,  Figure  2-  3  8,  designed  at  Cambridge  (Wroth 
and  Hughes;  1972,  1973).   This  is  a  heavy  wall  open  tube 
carrying  flush  mounted  lateral  stress  gages  and  an  internal 
auger  to  reduce  insertion  disturbance.   The  time  required 
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for  the  soil  stresses  to  relax  is  on  the  order  of  10  hours 
as  opposed  to  up  to  10  days  for  the  NGI  stress  probe. 
Typical  data  are  shown  in  Figures  2-39  and  2-40.   The 
device  is  similar  to  the  Cambridge  pressuremeter  in  its 
manner  of  insertion.   At  one  point,  both  the  pressure- 
meter  and  stress  probe  units  were  mounted  end  to  end. 
Because  of  the  difference  in  test  duration,  the  pressure- 
meter  is  probably  more  efficiently  used  in  a  separate  probe 
even  if  both  units  are  used  at  the  same  site. 

The  stress  probe  configurations  reported  in  the 
literature  surveyed  monitor  lateral  stress  in  only  one 
direction.   If  the  lateral  stresses  are  not  equal  and  the 
principal  directions  are  not  known,  the  devices  have  to  be 
rotated  in  the  ground  to  make  three  separate  sequential 
measurements . 

2.4.4  HYDRAULIC  FRACTURING  METHODS 

A  common  stress  measurement  tool  in  rock  mechanics 
is  hydraulic  fracturing  (Haimson,  1968,  1971,  1972;  Kehle, 
1964;  Obert  and  Duvall,  1967).   Initially  employed  as  a 
technique  to  improve  well  production  by  fracturing  the 
near-by-strata,  the  technique  may  be  extended  to  measure 
the  least  principal  stress  and  rock  tensile  strength. 
The  technique  as  mathematically  idealized  by  Kehle  is  shown 
in  Figure  2-41.   A  section  of  the  well  borehole  is  isolated 
by  packers  and  pressurized  at  sufficiently  high  levels  to 
open  a  tensile  fracture  normal  to  the  least  principal  stress. 
The  shut-in  pressure  Pg  in  Figure  2-42  is  shown  to 
approximate  the  least  principal  stress.   If  the  least 
principal  stress  happens  to  be  vertical,  little  information 
is  gained  since  this  stress  component  is  generally  known, 
given  the  density  of  overburden  material.   Haimson  describes 
the  stress  information  that  can  be  retrieved  from  a  know- 
ledge of  the  fracture  direction  and  pressures: 
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For  the  most  common  case  of  vertical  fracture  initiation  and 
extension,  the  in-situ  stress  field  is  well  defined. 

The  use  of  hydraulic  fracturing  in  soils  apparently 
results  from  investigations  of  anomalous  values  of  field 
permeability  values  reported  by  Bjerrum  and  Anderson  (1972) 
and  Bjerrum,  et  al. ,  (1972) .   Bozozuk  of  the  National  Research 
Council  of  Canada  (NRC)  and  Schmertmann,  recently  at  NRC, 
are  extending  the  method.   This  work  employs  standard 
piezometers  that  are  backpressured  to  fracture  the  soil. 
Bjerrum' s  field  apparatus  and  data  are  shown  in  Figures  2-43 
and  2-44.   If  paraffin  is  used  as  a  fracturing  fluid  there 
is  no  flow  until  the  fracture  actually  occurs.   If  water 
is  used,  minor  flow  into  the  soil  occurs  before  fracture. 
This  is  the  pressure  range  that  would  be  employed  in  a 
permeability  test.  When  fracture  does  occur,  the  flow 
rate  is  drastically  increased. 

Bozozuk 's  considerably  refined  apparatus  is  said 
to  be  capable  of  measuring  lateral  stress  and  tensile 
strength  for  depths  down  to  125  ft  in  fine  grained  soils 
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having  K  as  high  as  1.25.   In  general,  lateral  stresses 
in  excess  of  the  overburden  pressure  can  be  measured  only 
if  the  soil  has  some  tensile  strength.   One  of  the  major 
advantages  of  the  technique  is  the  rapidity  with  which 
measurements  can  be  made  once  the  piezometers  are  established 
in  place. 

Hydraulic  fracturing  would  be  highly  recommended  for 
fine  grained  soils  if  the  in-situ  stresses  are  to  be  moni- 
tored at  a  point  during  some  period  of  time.   Once  the  probe 
is  in  place  and  the  soil  stresses  have  relaxed,  measurements 
can  be  made  at  the  rate  of  20  a  day.   It  would  appear  that 
major  earthwork  operations  particularly  lend  themselves 
to  this  technique  for  monitoring  soil  stress  build-up,  if 
the  degree  of  compaction  is  not  so  great  as  to  result  in  the 
minor  principal  stress  being  in  the  vertical  direction. 

Successful  use  of  the  technique  in  coarse- grained 
soils  is  limited  in  that  the  piezometer  tip  should  have 
high  permeability  as  compared  with  the  soil.   Certain 
approaches  are  obvious,  such  as  the  use  of  a  relative 
viscous  fracturing  fluid  in  an  open  ported  piezometer  tip. 
It  is  not  known  if  these  developments  have  been  attempted. 

2.4.5  VANE  SHEAR  TESTS 

The  strength  anistropy  of  a  soil  mass  may  be 
measured  by  means  of  multiple  vane  shear  tests  (Aas ,  1965, 
1967;  and  Blight,  1971)  using  vanes  of  different  aspect 
ratio.   The  concept  is  shown  in  Figure  2-45.   The  torque 
M  required  to  rotate  a  vane  of  diameter  D  and  height  H  is 
written  in  terms  of  the  horizontal  and  vertical  shear 

strengths  S,  and  S  : 

°     h      v 


M  = 


M 


7TDH 

-  S     +2  £±L-     _ 
2     v              4        3 

sh 

(after  Aas) 

2 

=     S     +  S, 
v          h 

"d" 

3H 

ttD2H 
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Then  if  vane  shear  tests  having  different  aspect  ratios 
D/H  are  plotted  on  2M/ttD  H— D/3H  coordinates,  the  result 
should  be  a  straight  line  having  intercepts  of  S  and 
S,/S  on  the  two  axes.   Data  from  four  Norwegian  sites  as 
reported  by  Aas  are  shown  in  Figure  2-46.   It  can  be  seen 
that  the  straight  line  relationship  holds  except  for  high 
aspect  ratios  (D/H  -  1/6).   The  values  for  Sv/Sh  =  1/(S,/S  ) 
range  from  about  1/2  for  normally  consolidated  sites  to 
1  for  slightly  overconsolidated  sites.   These  are  roughly 
the  expected  values  of  K  ,  so  it  appears  that  the  measured 
strength  anisotropy  may  reflect  in-situ  stress  ratios.   Aas 
states  that  this  is  supported  by  tests  in  a  vane-triaxial 
apparatus;  i.e.,  a  triaxial  cell  into  which  a  minature  vane 
is  inserted. 

A  somewhat  different  analysis  is  made  by  Blight 
(1971) .   The  torque  required  to  shear  vanes  of  various 
aspect  ratios  is  written  in  terms  of  the  cohesion,  friction 
angle  and  vane  dimensions.   The  following  notation  is  used: 

T  (D)  =  torque  for  vane  aspect  ratio  of  1 

T(D/3)  =  torque  for  vane  aspect  ratio  of  1/3 

T(2D)  =  torque  for  vane  aspect  ratio  of  2 

R  =  strength  ratio  S  /S, 

m     T(D)  -  T(D/3)       m     R  =  T(2D)  -  T(D/3) 
3T(D/3)  -  T(D)  6T(D/3)  -  T(2D) 

If  0  is  assumed  isotropic  and  c  is  assumed  negligible  in 
comparison  with  a  tan  0,  then 

c  +  o,     tan  0   cj-u, 

R  =  D =  -±  . 

c  +  a     tan  0   crv 

Blight  cautions  that  his  assumptions  may  only  be  valid  for 
intrinsically  isotropic  soils  in  which  cohesion  may  be 
ignored. 
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2.4.6   COMPARISON  OF  FIELD  DETERMINATIONS 

Five  in-situ  techniques  currently  available  for 
determining  in-situ  field  stresses  are: 

acoustic  velocity 
borehole  pressuremeter 
borehole  stress  probe 
hydraulic  fracturing 
anisotropic  vane  shear 

The  choice  of  a  particular  technique  for  use  at  a  specific 
site  depends  largely  on  the  circumstances  and  must  be  made 
after  a  number  of  factors  have  been  considered.   The  geo- 
technical  applicability  of  each  field  technique  is  indicated 
in  Table  2.1.   Matrix  locations  marked  X  indicate  geotechnical 
conditions  for  which  a  technique  is  particularly  suited, 
conditions  noted  N/A  indicate  geotechnical  conditions  under 
which  a  technique  is  not  appropriate.   Selection  of  a  tech- 
nique for  use  at  a  particular  site  requires  judgements 
that  cannot  be  included  in  a  simple  table  --  engineering 
judgement  requires  familarity  with  the  method,  the  site  and 
the  apparatus  available. 

Some  of  the  characteristics  of  the  existing  in-situ 
methods  are  summarized  in  Table  2.2.   The  accuracy  of  each 
method  obviously  depends  on  the  circumstances  of  the  test, 
and  the  ratings  applied  may  change  if  the  geotechnical 
conditions  are  not  appropriate  to  the  test  method  being 
used.   The  range  of  each  test  is  more  easily  identified. 
The  borehole  and  anisotropic  vane  shear  devices  respond 
to  local  effects  near  the  test  apparatus,  and  will  require 
great  care  to  avoid  insertion  disturbance.   Acoustic  velocity 
and  hydraulic  fracturing  are  long  range  tests  that  tend  to 
average  over  a  large  volume  of  soil  and  be  less  influenced 
by  local  disturbance.   The  uniformity  of  the  soil  mass  must 
also  be  considered  in  terms  of  the  sampling  range,  since 
a  long  range  technique  will  not  resolve  the  detailed  changes 
in  stress  that  may  be  expected  if  more  than  one  soil  type 
is  present  in  the  volume  sampled. 
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Table  2.1 

APPLICABILITY  OF  FIELD  METHODS  TO 
VARIOUS  GEOTECHNICAL  CONDITIONS 
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Table  2.2 


CHARACTERISTICS  OF  IN- SITU  STRESS  DEVICES 
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The  costs  associated  with  the  various  methods 
include  the  capital  expenditures  and  the  time  required  to 
conduct  the  actual  testing.   The  anisotropic  vane  shear 
test  requires  only  several  special  vanes  in  addition  to 
standard  equipment.  The  capital  cost  associated  with  this 
method  is  noted  as  negligible  in  Table  2.2,  indicating 
costs  less  than  $100.  The  capital  costs  associated  with 
the  cross -borehole  shooting  and  hydraulic  fracturing 
methods  range  from  minor  cabling  to  more  significant 
recording  and  transducer  purchases.  The  following  major 
items  are  suggested  for  the  two  methods: 

cross -borehole  shooting  hardware 

acoustic  receiver  (geophone) 

acoustic  transmitter 

oscilloscope  (preferably  dual  beam) 

adjustable  delay  line 

cabling 

hydraulic  fracturing  hardware 

standard  piezometer  tip 

electrical  pressure  transducer  (mounted  down 

hole) 
conditioning  electronics  for  transducer 

(bridge  balance) 
recorder  (X-Y  or  chart  recorder) 
variable-flow  fluid  source  (e.g.  ,  pressurized 

gas  bottle,  precision  regulator,  fluid 

bottle  and  metering  valve) 
flowmeter 

hydraulic  fracturing  hardware  (shallow  depths  only) 

standard  piezometer  tip 
mercury  manometer 
screw  volume  control 

Obviously  the  capital  expenditures  depend  on  the  instrumenta- 
tion available  and  the  sophistication  desired.   Total  system 
costs  range  from  several  hundred  to  several  thousand  dollars. 

The  techniques  requiring  the  most  expensive  hardware 
are  the  cross -borehole  pressuremeter  and  the  stress  probe. 
In  addition  to  instrumentation  requirements,  a  sophisticated 
down  hole  device  is  required.   Both  techniques  employ 
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self-boring  probes  and  presently  only  one  model  of  each  is 
known  to  be  on  the  market. 

The  operating  costs  associated  with  each  technique 
depend  on  the  structure  and  work  load  of  the  field  crew. 
Cross -borehole  shooting  probably  requires  only  several 
minutes  per  test  plus  drilling  time  to  reach  the  desired 
level.   The  anisotropic  vane  shear  technique  is  considered 
equivalent  to  three  standard  vane  shear  tests.   Hydraulic 
fracturing  requires  less  than  1/2  hour  per  test  after  the 
piezometers  are  established  and  excess  pore  pressures  due 
to  driving  are  allowed  to  dissipate.   The  self -boring 
pressuremeter  requires  1/2  to  several  hours  per  test  whereas 
the  stress  probe  requires  one  to  four  days.   The  latter  two 
instruments  also  require  drill  rig  and  augering  support  for 
insertion  at  the  desired  horizon. 

2.5  EVALUATION  OF  ANISOTROPIC  PROPERTIES 

2.5.1   INTRODUCTION 

Because  soil  fabric  tends  to  become  more  oriented 
with  increasing  consolidation  or  deformation,  the  stress 
history  is  inferred,  at  least  in  part,  from  the  degree  and 
orientation  of  anisotropic  behavior  in  a  soil  mass.   The 
fabric  elements  that  display  stress  dependency  thus  are  of 
interest  vary  considerable  in  size.   For  example,  Blight 
(1971)  applied  the  Griffith  criterion  to  determine  the 
fracture  patterns  expected  in  shrinking  or  swelling  clays. 
By  this  criterion,  the  angle  f   to  the  vertical  at  which 
cracking  occurs  is  given  by: 

av-Qh 
cos  2if   =  — - — - —  shrinking 


i     i 

cr,  -  a 

cos2(90°-^)   =  — 2- 1  swelling 

2K+CTv> 
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The  case  of  shrinking  is  given  by  a  mass  of  initially 
intact  clay,  saturated  throughout,  but  with  the  water 
table  a  distance  h  below  the  surface.   The  pore  water 
tension  in  the  saturated  clay  at  the  surface  is  given  by 
u  =  ay  gh«   Tne  factor  a  tends  toward  unity  if  evaporation 
ceases,  and  increases  above  unity  if  the  evaporation  rate 
increases.   The  resulting  shrinkage  stress,  which  depends 
on  the  stiffness  of  the  soil,  may  be  defined  by  a  -   PY,gh, 

£3  W 

in  which  p  will  probably  be  equal  to  or  greater  than  a. 
The  effective  and  total  stresses  may  now  be  computed, 
giving  for  the  Griffith  criterion 


cos  2ijr  =  2 —  at  the  surface 

2(2a-p) 


Yz(Ko-l)-PYw(h-z) 

cos   2tjr   =      -»g ■  '  '  '    -i       at   depth  z 


2[Yz(Ko+l)-(2a-P)Yw(h-z)] 


or  if  it  is  assumed  that  K  =  %  and  v  =  2y 

o        z  w 

-  [z+e(h-z)l 

cos   2i|r   =     — : — - ,  at  depth  z 

2  [3z-(2a-B)(h-z)| 

Predicted  crack  profiles  for  two  values  of  a   and  8  are 
shown  in  Figure  2-47,  a  corresponding  actual  crack  pattern 
in  Figure  2-48. 

If  shrinkage  cracks  become  filled  by  material  washed 
down  from  the  surface  and  the  soil  subsequently  swells, 
horizontal  swelling  will  be  restrained  and  horizontal  swelling 
pressures  may  uecome  large.   Two  sets  of  fissures  may  be 
expected  to  form  in  soil  that  alternately  shrinks  and  swells. 
The  swelling  cracks  will  form  at  orientations  of  50  to  60 
degrees  if  the  cracks  form  at  a  total  stress  ratio  of 
3<  ah/Qv  ^  ^'      "tn  F*-8ure  2-49,  a  frequency  distribution 
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of  fissure  angles  shows  a  concentration  of  shrinkage  cracks 
at  f   =  30  degrees  and  swelling  cracks  at  f   =  57  degrees. 
It  can  be  seen  that,  although  anisotropy  is  generally 
considered  to  be  the  result  of  the  orientation  of  small 
fabric  elements  such  as  soil  grains,   fabric  elements  many 
times  larger  may  be  useful  in  deducing  the  history  of  a 
soil  mass. 

The  following  tools  and  disciplines  were  reviewed 
to  evaluate  their  potential  use  with  respect  to  fabric 
analysis  and  the  broader  goal  of  in-situ  stress  determination, 

o  hydraulic  permeability 

o  electrical  conductivity 

o  acoustic  propagation 

o  electron  microscopy 

o  infrared 

o  gamma  rays 

o  electron  beam 

o  electromagnetic  resonance 

2.5.2  HYDRAULIC  PERMEABILITY 

Hydraulic  permeability  is  a  parameter  routinely 
measured  both  in-situ  and  in  the  laboratory.   The  anisotropy 
of  permeability  has  been  shown  to  significantly  influence 
soil  properties.   Ingles  and  Lee  (1971)  argue  that  strength 
anisotropy  in  fine-grained  soils  results  from  the  influence 
of  macropore  anisotropy  on  minimum  energy  fracture  paths. 
Duncan  and  Seed  (1966)  show  that  anisotropy  in  Skempton's 
pore  pressure  coefficient  at  failure  A^, 

Au  =  B  (Ao^)  +  Aikv^   -  Aa3) 

is  associated  with  anisotropy  of  undrained  strength. 
Experimental  data  given  show  that  Af  for  horizontal  samples 
is  consistently  greater  than  Af  for  vertical  samples  as 
indicated  in  Figure  2-50.  Since  Af  depends  on  the  soil 
matrix  stiffness  and  pore  geometry,  the  different  ~K~ 
values  in  the  two  directions  indicate  the  corresponding 
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differences  in  fabric  and  stresses.   It  is  anticipated 
that  the  soil  pore  fabric  may  prove  to  be  a  more  sensitive 
stress  indicator  than  the  soil  grain  fabric. 

Mercury  intrusion  techniques  were  employed  by 
Diamond  (1970)  and  Sridharan,  Altschaeffl  and  Diamond  (1971) 
to  determine  pore  size  distribution  curves.   As  static  compac- 
tive  effort  increases,  the  larger  voids  are  sequentially 
eliminated,  the  smaller  voids  being  largely  unaffected. 
The  pore  size  distribution  for  sedimented  clay  is  signifi- 
cantly different  from  that  in  compacted  clay,  with 
few  voids  larger  than  the  grain  size  (Figures  2-51 
and  2-52) .   Pore  size  distributions  in  cohesionless  material 
were  determined  by  Windisch  and  Soulie  (1970)  by  manual 
measurements  on  thin  sections  of  epoxy  treated  specimens. 
Few  data  except  a  pore  size  distribution  are  given  in  this 
paper  which  is  basically  the  description  of  a  fabric  analysis 
technique . 

Possible  methods  by  which  pore  fabric  may  be 
evaluated  include  miscroscopy,  hydraulic  permeability  and 
electrical  conductivity.   Hydraulic  permeability  is 
determined  in  the  laboratory  by  various  techniques  or  by 
piezometric  measurements  in  the  field.   The  anisotropy  of 
permeability  is  commonly  treated  by  assigning  separate 
values  to  the  permeability  coefficient  in  the  horizontal 
and  vertical  directions.   Taylor  (1967)  gives  a  model  for 
the  hydraulic  permeability  coefficient  kH  as 

*B.     =     Cs  ^w  n  4  I   * 

The  electrical  analog  for  laminar  flow  is  employed  by  Long 
and  Zimmie  (1973)  to  evaluate  the  shape  factor  C_  and 
estimate  the  hydraulic  radius  R„.   The  DC  electrical  flow 
in  a  porous  plug  can  be  expressed  by 
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C  n  = 


s        AvXA 


This  may  be  compared  with  the  fluid  flow  as  follows 

Cn=   JL  .L-  _L  =  Jl  JL  _JL__ 

S        Av   A    x      AH   A   Y  R2 


w  H 


or 


R^    =   2.  Av  r±\ 


i  AH  y 
w 


in  which  IL,  =  hydraulic  radius 
L   =  sample  length 

AH 
Q   =  hydraulic  flow  rate  =  k„  ~  A 

AH  =  fluid  head  loss 

r\        -     viscosity  of  fluid 

Yr.  =  unit  weight  of  fluid 
w  ° 

i   =  DC  current  flow 

Av  =   impressed  voltage 

X   =  specific  conductivity  of  electrolyte 

n   -  porosity 

A   =  area 
Sequential  measurements  of  fluid  flow  and  DC  current  were 
made.   Tests  on  fretted  glass  disks  were  conducted  and 
found  to  compare  reasonable  well  with  mercury  intrusion 
data  and  the  advertised  properties  of  the  glass  disks.   No 
measurements  in  soil  samples  were  reported.   Seemingly  soil 
anisotropy  would  be  strongly  displayed  by  the  use  of  this 
technique  in  samples  of  oriented  soil  taken  in  different 
directions . 

Application  of  pore  anisotropy  as  a  fabric  parameter 
may  be  implemented  by  two  approaches.   The  permeability  k„ 
may  be  evaluated  and  empirically  correlated  to  grain 
orientation  functions  by  assumption  of  the  shape  factors 
C  for  the  soil  in  question.   This  approach  may  be  useful 
after  a  significant  backlog  of  data  is  established.   No 
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published  relations  between  permeability  ratio  and  grain 
orientation  ratio  were  discovered  during  the  literature 
survey.   Alternatively,  combined  permeability  and  conductivity 
measurements  may  be  used  to  determine  the  effective  hydraulic 
radius  It,  in  various  directions,  and  a  more  direct  correla- 
tion made  to  fabric  parameters.   It  is  possible  that  sounding 
apparatus  may  be  devised  to  permit  in-situ  measurements. 

2.5.3   ELECTRICAL  CONDUCTIVITY 

Specific  values  of  the  resistivity  of  soil  cannot 
be  tabulated  to  identify  soil  types.  According  to  Moore 
(1961)  ,  most  clays  will  display  low  resistivities  ranging 
from  several  to  several  tens  of  ohm-meters.   Clean  sand 
and  gravel  may  range  up  to  10,000  ohm-m.   Dense  sand  or 
rock  such  as  granite  may  have  resistivities  of  10,000  to 
40,000  ohm-m,  whereas  hard  volcanic  rock  formations  may 
be  as  low  as  50  ohm-m. 

Electrical  conductivity  has  been  employed  as  a  tool 
for  the  evaluation  of  soil  micro structure  by  several  investi- 
gators. Mitchell  and  Arulanandan  (1968)  evaluated  the  soil 
properties  that  influence  electrical  dispersion.   The  soil 
was  modeled  as  a  parallel  R-C  circuit,  and  the  parameters  e1 
(real  portion  of  dielectric  constant) ,  e"  (imaginary  portion 
of  dielectric  constant)  and  a  (conductivity  in  mhos/m) 
were  plotted  as  a  function  of  frequency  (Figure  2-53) . 
Anisotropy  was  found  to  influence  Ae  ,  f  ,  and  a. 

The  value  Ae'  reflects  the  magnitude  of  polarization 
developed  at  low- frequencies ,  and  thus  is  a  measure  of 
mechanisms  requiring  relatively  long  times,  such  as  the 
displacement  of  double-layer  charges  adjacent  to  particle 
surfaces . 

The  characteristic  frequency,  f  ,  is  a  measure  of 
the  average  time  for  a  particular  polarization  process „   This 
reflects  the  distance  and  ease  with  which  charges  may  move, 
and  charge  size.   The  more  easily  charges  can  follow  the 
electric  field,  the  higher  will  be  f  . 
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Figure  2.53. 


Conductivity  &  Dielectric  Dispersion 
Characteristics  of  Saturated  Kaolinite 
Hydrite  R 
(Mitchell  and  Arulanandan,  1968) 


ELECTRICAL  PROPERTIES  OF  SAN  FRANCISCO  BAY  MUD 


Ae  *    f Q     a 
(x  10"4)  (Hz)   (deg) 


'DC 


UDC 
mhos/m 


Bay  Mud  parallel        100    300     18    0.96 
Bay  Mud  remolded        116    200     14    0.97 
Bay  Mud  perpendicular    82    500      5    0.93 
(Moisture  Content  ■  90%) 
(Pore  Fluid  Electrolyte  Contentration  =  0.3  N) 

(Mitchell  and  Arulanandan,  1968) 
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The  DC  conductivity,  o^,  measures  the  amount  of 
free  charge  and  the  ease  with  which  it  passes  through  the 
soil.   Thus  o~c   may  reflect  void  ratio,  flow  path  tortuosity 
and  pore  fluid  conductivity.   An  increase  in  conductivity 
with  frequency  is  represented  by  o >„   -  onr.   The  major 
portion  of  a*„   is  ascribed  to  electrokinetic  coupling 
between  current  and  fluid  flow. 

If  several  polarization  mechanisms  are  operative, 
polarization  will  occur  over  a  larger  frequency  range;  if  a 
single  uniform  mechanism  exists,  polarization  will  be  confined 
to  a  narrow  band  of  frequencies.   This  distribution  of 
relaxation  times  is  measured  by  the  parameter  a,  generated 
by  a  Cole-Cole  plot  of  e"  versus  e.'  .   The  wider  is  the 
distribution  of  relaxation  times,  the  larger  will  be  a. 

The  parameters  a  and  oJ)C   showed  consistent  trends,  whereas 
remolding  may  have  disrupted  cation  and  absorbed  water 
layers,  significantly  affecting  Ae1  and  f  .   The  apparatus 
used  to  make  resistance  and  capacitance  measurements  required 
manual  null -balancing  of  a  physical  R-C  circuit,  and  it  is 
likely  that  development  of  the  method  will  result  in 
specialized  electronics  that  are  much  easier  to  employ. 

The  difficulty  of  employing  null-balance  methods 
was  cited  by  Ramiah  and  Purushothamaraj  (1971)  as  the 
reason  for  measuring  only  o-qc   in  their  application  of  the 
method  to  direct  shear  testing.   Direct-current  dispersion 
characteristics  were  employed  as  a  -fabric  parameter  in  a 
series  of  direct  shear-reversal  tests  directed  at  the 
evaluation  of  initial  structure  on  residual  strength.   In 
this  study  the  DC  conductivity  was  measured  in  a  direct  shear 
box  along  a  path  crossing  the  shear  plane  at  an  angle  of 
nearly  zero.   The  conductivity  normal  to  the  shear  plane 
was  not  measured.   These  data  indicate  that  conductivity 
decreases  with  increasing  alignment  of  soil  particles  to 
the  current  flow,  counter  to  the  results  in  Bay  Mud.  A 
rough  correlation  between  stress  ratio  and  conductivity  ratio 
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Figure  2-54.  Typical  Stress-Deformation  Curve 
(Ramiah  &  Purushothamaraj  1971) 
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Figure  2-55. 


DC  Conductivity  vs  Deformation 
(Ramiah  &  Purushothamara,  1971) 
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was  also  found.   Portions  of  these  data  are  shown  in 
Figures  2-54  and  2-55.   The  decrease  in  conductivity  is 
attributed  to  electrostatic  shielding  resulting  from  the 
parallel  alignment  of  clay  particles  in  the  shear  zone,  so 
that  charge  carriers  housed  outside  that  zone  require  greater 
transit  times.   Many  factors  influence  the  response  of  a 
clay-water-electrolyte  system,  e.g.,  the  soil  fabric,  amount 
and  type  of  electrolyte,  moisture  content,  the  nature  of  the 
clay  surface  and  temperature.   The  change  in  conductivity  as 
soil  fabric  changes  depend  upon  the  relative  importance  of 
mechanisms  that  increase  or  decrease  charge  carrier  flow  in 
the  particular  case. 

Electromagnetic  techniques  are  not  limited  to  DC 
jr  low  frequency  methods.   The  literature  survey  included 
the  entire  electromagnetic  spectrum  from  DC  to  electron 
frequencies.   Later  discussion  will  include  the  use  of 
higher  frequencies. 

2.5.4   ELECTRON  MICROSCOPY 

The  electron  microscope  must  be  considered  to  be 
a  tool  for  use  in  fabric  analysis.   Many  excellent  micro- 
graphs of  kaolin  are  available  in  the  literature, 
but  montmorillonite  and  illite  micrographs  are  rare. 
The  scanning  electron  microscope  produces  excellent  repre- 
sentations of  kaolinite.   These  may  be  viewed  directly,  or 
grain  orientation  may  be  quantified  by  measuring  grain  angles 
or  spacing.   This  is  frequently  done  manually,  but  computer 
scanning  techniques  could  be  applied.   Various  optical 
conditioning  methods  are  also  available,  e.g.,  the  optical 
diffraction  and  convolution  square  techniques  of  Tovey  (1971). 

We  will  consider  microscopy  only  as  a  fabric  analysis 
tool.   Geotechnical  microscopy  may  be  classed  as  dealing  with 
structure  or  as  corroborating  the  results  of  other  techniques, 
as  the  use  of  x-ray  diffraction  in  quantitative  analysis. 
Discussions  of  micros true ture  given  by  Mor gens tern  and 
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Tchalenko  (1967);  Lafeber  and  Kurbanovic  (1965);  Barden, 
Sides  and  Karunaratne  (1970)  ,  Barden  (1962) ;  and  Kirkpatrick 
and  Rennie  (1972)  deal  more  with  the  factors  that  result  in 
anisotropy  of  soil  fabric.   Microscopy  is  used  as  a  tool  to 
corroborate  x-ray  diffraction  studies  by  Comer  (1960,  1962), 
and  Ormsby  and  Shartsis  (1960) . 

Full  utilization  of  the  scanning  electron  microscope 
requires  that  some  quantitative  parameter  representing  the 
soil  fabric  be  commonly  employed.   The  literature  indicates 
that  this  is  not  presently  the  case.   It  is  clear  that 
changes  in  soil  fabric  occur  under  stress,  but  different 
parameters  are  measured  by  various  investigators. 

Investigating  the  deformation  fabric  of  sand,  Oda 
(1972)  measured  grain  orientation  distribution  and  the 
ratio  of  grain-grain  contact  areas  projected  on  the  vertical 
and  horizontal  planes.   The  angle  9  between  the  long  axis 
of  the  grains  for  two  degrees  of  initial  orientation  are 
shown  in  Figure  2-  56  .   In  Figure  2 -56(a),  the  sand  is  initially 
oriented  horizontally,  and  the  preference  for  this  direction 
becomes  stronger  with  increasing  stress  ratio.   This  prefer- 
ence decreases  slightly  at  large  axial  strains.   In  Figure 
2-56(b),  the  initially  vertically  oriented  sand  also  tends 
toward  a  horizontal  orientation  with  increasing  axial  strain, 
with  the  result  that  the  initially  oriented  structure  is 
destroyed. 

The  projected  contact  area  ratio  Sz/Sx,  in 
(Figure  2-57 (a))  was  found  to  change  with  stress  ratio  as 
shown  in  Figure  2-57 (b) .   An  important  fact  is  that  the 
parameter  changes  continuously  with  stress  ratio  and  does 
not  display  changes  in  behavior  associated  with  peak  stress. 
Murayama  and  Matsouka  (1973)  and  Matsouka  (1973)  argue  that 
the  important  parameter  is  not  the  projected  area  ratio  S  /S  , 
but  the  distribution  of  the  contact  normal  directions.   In 

fact,  Oda  assumes  a  mean  contact  area  and  his  ratio  S  /S 

'  z'  x 

appears  to  be  a  function  of  the  contact  directions  only. 
Murayama  and  Matsouka  conducted  direct  shear  tests  on  rod 
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Frequency  Distribution  of  9  Alunimun  Rod  Mass 
(a)  Before  Shear  (b)  At  the  Peak  Strength 
(c)  At  the  Residual  Strength 
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masses  and  found  a  nonsymmetric  distribution  of  contact 
normal  directions  as  shown  in  Figure  2-58.   The  different 
behavior  reported  in  these  two  papers  may  be  a  result  of 
the  different  forms  of  symmetry  associated  with  triaxial 
tests  as  opposed  to  direct  shear  tests.   It  is  clear,  however, 
that  fabric  changes  associated  with  changes  in  stress  may 
be  observed  microscopically,  but  that  additional  study  will 
be  necessary  before  the  scanning  electron  microscope  can 
be  confidently  used  to  measure  stress-dependent  parameters. 
2.5.5   INFRARED 

The  use  of  infrared  on  soils  has  been  limited 
primarily  to  quantitative  analysis  to  determine  the  amount 
of  various  minerals  present  under  various  conditions,  e.g., 
MacKenzie  (1969) .   It  would  appear,  however,  that  IR  has 
some  application  in  fabric  analysis.   The  use  of  infrared 
techniques  is  attractive  because  wavelengths  in  this  portion 
of  the  spectrum  are  of  the  same  length  as  clay  size  particles. 
Thus,  optical  techniques  employing  IR  radiation  may  be  quite 
effective  in  soils.   A  wide  variety  of  laser  sources  is  now 
available-  so  that  illumination  with  wavelengths  between 
0.3  to  10.6  microns  may  be  employed.   Two  major  alternatives 
exist:   direct  detection  techniques,  such  as  probing  an  area 
with  beam  of  polarized  light,  and  pattern  analysis  based  on 
IR  images. 

Detection  techniques  may  be  based  either  on  the 
principle  of  diffraction  or  polarization.   If  the  soil  matrix 
is  considered  to  be  a  diffraction  grating,  the  diffraction 
depends  on  the  ratio  of  grain  spacing  to  wavelength.   In  a 
highly  oriented  soil,  the  grain  spacing  is  much  smaller 
when  the  plate-like  grains  are  viewed  edge-on.   Thus, 
qualitative  evaluation  of  the  degree  of  grain  orientation 
and  its  direction  should  be  immediately  available  from  the 
intensity  and  spacing  of  the  diffraction  pattern.   A  similar 
direct  approach  may  use  a  reflection  polariscope  or  ellipso- 
meter  to  measure  the  net  refractive  index  on  specimen 
surfaces  in  various  directions. 
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It  is  likely  that  extraction  of  quantitative  data 
with  direct  detection  techniques  will  be  difficult  and  that 
optical  or  electronic  image  processing  may  be  necessary „ 
This  is  the  second  major  alternative  of  image  analysis. 
Much  previous  work  in  the  application  of  fabric  analysis 
has  suffered  from  the  lack  of  rapid  and  exact  techniques 
for  measuring  the  soil  fabric.  When  actual  measurements 
are  made  on  soil  fabric  photographs,  they  are  frequently 
made  manually  at  what  must  be  a  relatively  high  cost  for 
the  amount  of  information  returned.   If  optical  or  computer 
image  conditioning  techniques  are  to  be  employed,  the 
advantages  of  IR  techniques  relative  to  visible  light  are 
decreased.   The  additional  difficulties  associated  with 
IR  technology  would  outweigh  any  advantages  of  increased 
information  (i.e.,  over  that  obtained  by  optical  diffrac-. 
tion  using  visible  light) . 

A  potential  use  of  IR  is  in  the  evaluation  of 

the  extent  of  diagenetic  bonding  in  particular  soil  specimens 
The  use  of  transmission  infrared  in  the  determination  of 

the  clay  content  of  soils  would  involve  the  use  of  the 
Si-0  absorption  band  in  the  9  to  llu  region.   This  band  is 
of  strong  intensity  and  could  possibly  serve  in  both 
qualitative  and  quantitative  roles.  However,  this  portion 
of  the  spectrum  is  a  general  group  frequency  range  where 
many  oxygenated  species  other  than  silicates  also  absorb  and 
consequently  would  interfere.  The  presence  of  such  anions 
as  borates,  phosphates,  sulfites,  sulfates  and  chlorates 
would  limit  the  use  of  this  band.   Accordingly,  the  use  of 
transmission  infrared  spectroscopy  would  be  subject  to  the 
presence  or  absence  of  a  large  number  of  inorganic  constitu- 
ents which  may  be  present  in  a  particular  soil  sample. 

2.5.6  GAMMA  RADIATION 

Gamma  rays  are  commonly  used  in  the  determination 
of  density  and  moisture  content  in  soils  and  rock.   A  number 
of  available  nuclear  gage  configurations  use  transmission 
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and/or  backscatter  techniques  for  cased  or  uncased  boreholes , 
ground  level  or  special  applications.   These  include  perme- 
ability measurements  using  tagged  fluid,  determination  of 
the  average  atomic  number  of  the  soil  media,  and  its 
elemental  composition.   No  techniques  by  which  the  soil 
stresses  can  be  easily  determined  are  available,  however. 

2.5.7  X-RAYS 

X-rays  are  used  to  monitor  the  motion  of  x-opaque 
markers  placed  in  a  soil  sample  under  test  and  to  determine 
the  orientation  of  crystal  planes  by  diffraction.   The  use 
of  x-rays  to  monitor  internal  deflections  does  not  appear 
to  be  of  great  usefulness  in  an  in-situ  stress  measurement 
system.   This  is,  however,  a  useful  technique  in  laboratory 
studies  of  anisotropic  soil  properties  and  is  mentioned  by 
Wroth  and  Hughes  (1972,1973),  Gerrard  and  Morgan  (1972), 
Arthur  and  Shamash  (1970)  ,  Arthur  and  Menzies  (1972) ,  and 
others . 

X-ray  diffraction  may  be  employed  to  evaluate  the 
mineralogy  of  soil  specimens  and  its  use  is  especially 
common  in  the  ceramics  literature  'e.g.,  Friedman,  1973; 
Bishop,  1958;  Kirkpatrick  and  Rennie,  1972;  Comer,  1960; 
Clark,  1970;  Dudley,  1970;  Meade,  1964;  Ormsby,  et  al. , 
1962;  and  Ormsby  and  Shartsis,  1960  . 

The  most  useful  form  of  x-ray  diffraction  is  the 
evaluation  of  preferred  orientation  ratios.   There  are 
several  forms  of  preferred  orientation  ratio  (POR) ,  each 
a  quantitative  measurement  of  clay  particle  orientation  which 
may  be  a  most  useful  fabric  analysis  parameter.   Because 
x-ray  diffraction  is  a  fabric  analysis  tool,  its  use  in 
the  context  depends  on  the  development  of  fabric  analysis. 
The  characteristics  of  x-ray  diffraction  include  the  ability 
to  make  quantitative  POR  determinations.   Sample  preparation 
is  somewhat  simpler  than  is  required  for  optical  methods 
because  the  x-rays  respond  to  a  layer  of  soil  up  to  50|jjn. 
Sophisticated  apparatus  and  techniques  are  required. 
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Meade  (1961)  suggests  that  the  ratio  of  the 
diffraction  pattern  at  the  001  and  020  peaks  be  compared 
on  vertical  and  horizontal  sections  to  obtain  a  quantitative 
orientation  ratio.   That  is,  the  orientation  ratio  is  given 
by 

ooiz 


020 
OR  =   — 


001  +  001 
,  x v_ 

020  +  020 
x     y 


in  which  001z  is  the  001  diffraction  peak  on  a  section 
cut  normal  to  the  z-axis  and 

020  is  the  020  diffraction  peak  on  the 
same  section. 

The  procedure  is  indicated  in  Figure  2-59. 

Martin  (1966)  employed  a  diffractometer  equipped 
with  a  pole  figure  device  in  his  study  of  kaolinite.   He 
defined  a  peak  ratio,  PR,  as  the  ratio  of  the  002  and  020 
reflections  as  shown  in  Figure  2-  60.  Clark  (1969)  modified 
Martin's  PR  ratio  by  taking  the  ratio  of  PR  on  orthogonal 
faces,  obtaining  a  preferred  orientation  ratio  POR.   This 
was  correlated  to  strain  as  indicated  for  various  test 
conditions  (Figure  2-61) . 

Clearly  x-ray  diffraction  can  be  used  to 
quantitatively  describe  the  degree  of  anisotropy  in  clay 
minerals.   The  direction  of  the  principal  stresses  may  be 
identified  within  several  degrees,  but  additional  effort 
will  be  required  before  quantitative  statements  regarding 
the  magnitude  of  the  principal  stresses  can  be  made. 

2.5.8  ELECTRON  BEAM  TECHNIQUES 

Electron  beam  techniques  are  summarized  in  a 
review  article  by  Tufts  (1969).   All  of  these  techniques 
are  influenced  by  the  short  range  of  the  electrons  (10  to 
15  mg/cm  for  electrons  with  100  keV  energy).   As  a  result, 
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Figure  2-59 


Measurement  of  Horizontal  Preferred  Orientation 
By  Comparison  of  X-ray  Diffraction  Peak  Heights 
(Meade,  1961) 
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Simple  Shear-50%  porosity 
Electrolyte  Solution 
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Preferred  Orientations  from  Distortion 
Experiments,  Plotted  as  a  Function  of  Total 
Strain  1  4-  e,/l  +  £0.  tn  -  approximate  slope 

of  straight  line  fitted  by  hand  to  data 
points.   A.  Simple  shear  experiments,  saturated 
samples,  50  percent  porosity.   B.  Flattening 
experiments,  saturated,  50  percent  porosity. 
C.  Simple  shear  experiments,  saturated,  58 
percent  porosity.   D.  Simple  shear  experiments, 
saturated  with  1  m  NaCl  solution,  50  percent 
porosity. 
(Clark,  1969) 


104 


the  beam  is  either  interacting  with  a  surface  layer  or 
passing  through  a  very  thin  (1.500  [±m)    sample.   Because 
tbin-wafer  samples  of  this  size  are  not  practical  for  soil 
core  samples,  surface  techniques  are  the  only  reasonable 
possibility.   Surface  preparation  is,  of  course,  critical 
if  bulk  characteristics  are  to  be  represented. 

The  closest  analogous  situation  in  conventional 
testing  of  metals  is  as  follows:   a  test  bar  is  fractured, 
producing  a  very  rough  irregular  surface.  What  is  the 
preferred  orientation  and  crystal  size  in  the  exposed 
surface?  Electron  and  x-ray  diffraction  can  provide  such 
information.  The  primary  difference  between  electron  and 
x-ray  diffraction  in  such  an  application  is  the  depth  of 
the  surface  layer  contributing  to  the  response.   X-ray 
penetration  is  generally  greater  than  electron  penetration. 
For  metal  fracture  it  is  important  to  restrict  the  measure- 
ment to  the  top  most  surface  layer,  i.e.,  use  electrons, 
in  order  to  correlate  the  fracture  with  the  local  structure. 
For  a  soil  core  sample,  this  would  not  be  necessary  since 
one  is  interested  in  bulk  properties.   The  ability  of 
electron  beam  techniques  to  emphasize  surface  characteris- 
tics is,  in  fact,  a  disadvantage  and,  as  a  result,  primary 
emphasis  should  be  placed  on  x-ray  diffraction  procedures 
(such  as  texture  analysis,  pole  figure  analysis)  which  are 
generally  easier  to  implement. 

2.5.9   SUMMARY  OF  ANISOTROPIC  PROPERTIES 
MEASUREMENTS 

Laboratory  tools  that  can  be  used  to  measure  the 
orientation  parameters  of  soil  fabric  include  x-ray 
diffraction  and  scanning  electron  microscopy.   Careful  sample 
preparation  is  required,  especially  for  the  SEM.   Field 
techniques  include  permeability  and  electrical  conductivity. 
Some  capability  is  available  in  IR,  gamma -rays  and  the 
electron  beam  but  in  general  these  techniques  are  more 
difficult  to  apply  and  return  less  information.   Quantitative 
measurements  of  soil  fabric  may  be  made  by  these  techniques 
as  a  connecting  link  between  measurements  of  quantities 
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(e.go,  pore  size  distribution  or  preferred  orientation 
ratio) ,  and  the  state  of  stress  that  induced  the  anisotropy 
is  established. 
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3„   FABRIC  ANALYSIS 

3.1   INTRODUCTION 

Fabric  analysis  is  the  analysis,  usually  but  not 
necessarily  quantitative,  of  the  spatial  data  for  the 
interior  of  the  material  (e.g.,  the  soil  sample) 
under  consideration.   Results  depend  on  the  methodology 
used  and  the  materials  to  which  the  fabric  analysis  is 
applied.,   Fabric  analysis  has  so  far  been  largely  a  strain 
analysis  tool  in  rock  mechanics.   Although  activity  in 
stress  analysis  has  increased  in  recent  years,  fabric 
analysis  in  soil  mechanics  has  not  advanced  as  far  as  in 
rock  mechanics. 

The  treatment  of  fabric  analysis  in  this  section 
is  more  applied  than  theoretical.  We  present  an  overview 
of  relevant  aspects  of  fabric  analysis,  and  proceed  to  a 
development  of  the  problems  and  possibilities  in  the  use 
of  fabric  analysis  to  study  the  state  of  stress  in  soils. 
These  include: 

o  Expectations  of  fabric  analysis 

o  Statistical  methods 

o  Anticipated  behavior  of  soil  types 

o  Orientation  of  measuring  devices  with 

respect  to  soil  fabric 
o  Experimental  work 
o  Plans  for  mapping  and  analyzing  soil 

fabric 
o  Field  sampling 

3.2  EXPECTATIONS  OF  FABRIC  ANALYSIS 

Fabric  analysis  applied  to  undeformed  materials 
provides  a  datum  to  which  analyses  of  deformed  materials 
can  be  compared.   From  such  analyses  of  undeformed  materials 
one  can  obtain  primary  or  initial  fabrics  (most  commonly 
exemplified  by  depositional  and  diagenetic  fabrics)  which 
provide  insights  into  the  mode  of  formation  of  the  materials. 
Primary  fabrics  may  also  permit  prediction  of  the  orientation 

107 


of  mechanical  anisotropy  under  subsequent  loading  and,  under 
some  circumstances  the  relative  degree  of  mechanical  aniso- 
tropy0   Figures  3-1  and  3-2  indicate  relations  between 
fabric  and  deformation  that  were  developed  in  rock  mechanics 
but  presumably  are  applicable  to  soil  mechanics „ 

Kinematic  fabric  analysis  is  essentially  a  geometric 
or  strain  analysis  approach  to  deformation,  based  on  the 
premise  that  the  symmetry  of  the  fabric  reflects  the  symmetry 
of  movement  responsible  for  the  fabric.   This  does  not  lead 
directly  to  the  determination  of  the  state  of  stress  in 
the  material,  although  it  may  lead  to  useful  inferences 
regarding  the  orientation  of  one  or  more  principal  stress 
axes.   In  rocks,  kinematic  analysis  has  been  practical  only 
where  deformation  has  been  intense;  it  is  likely  that  at 
least  some  soils  develop  useful  deformation  fabrics  with 
considerably  less  intense  deformation0   Figure  3-3  illustrates 
four  basic  types  of  fabric  symmetry  determined  from  geometric 
considerations.   Kinematic  analysis  has  also  been  used  to 
restore  fabrics  to  an  earlier  configuration  by  unrolling 
and  leveling  operations  and  may  be  useful  in  delineating 
an  inferred  primary  fabric  and  in  determining  if  more  than 
one  stage  of  deformation  has  taken  place. 

Dynamic  fabric  analysis  is  essentially  a  stress 
analysis  approach,  via  the  study  of  fabrics,  based  on  the 
premise  that  specimens  deformed  in  the  laboratory  can  be 
used  to  determine  the  orientation  and  relative  magnitudes 
of  the  principal  stresses  at  the  time  of  deformation  in  their 
presumably  comparable,  natural  counterparts.   Absolute 
magnitudes  of  the  principal  stresses  in  nature  cannot  be 
estimated  with  confidence  from  the  types  of  laboratory 
experiments  carried  out  so  far.   Relative  magnitudes  or 
ranges  of  magnitudes  may  sometimes  be  estimated  using  physical- 
chemical  principles.   Intense  deformation  is  not  a  necessary 
condition  for  application  of  dynamic  fabric  analysis  in  rocks, 
and  this  might  also  apply  to  some  soils.   The  resolution  of 
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Figure  3-1.  Kinematic  and  dynamic  petrofabrics  applied  to 
geologic  deformation.   Both  descriptive  and 
interpretive  phases  of  fabric  analysis  are 
indicated  (Friedman,  1964). 
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deformation  of  a  body.   Relation  of  initial 
and  observed  fabric  to  deformation  is  indi- 
cated  (Paterson  and  Weiss,  1961). 
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Figure  3-3.  Types  of  fabric  symmetry  (Turner  and  Weiss,  1963). 

(a)  and  (b)   axial  fabrics  (symmetry  of  an  ellipsoid  of 
rotation) 

(c)  orthorhombic  fabric  (symmetry  of  the  triaxial 
ellipsoid) 

(d)  monoclinic  fabric  (single  plane  of  symmetry) 

(e)  triclinic  fabric  (no  planes  of  symmetry) 

(Spherical  fabric,  viz.,  no  preferred  orientation,  not 
illustrated.   The  above  fabric  terms  are  not  equivalent 
to  crystallographic  terminology.,) 
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two  or  more  stress  systems  by  dynamic  fabric  analysis  is 
at  best  formidable  and  usually  not  feasible.   The  deforma- 
tion mechanisms  of  importance  here  are  fracture,  faulting 
and  flow  (cataclastic,  gliding,  recrystallization) . 

Under  very  favorable  circumstances,  it  should  be 
feasible  to  estimate  orientations  of  principal  stress  axes 
from  soil  fabric  analysis,  to  within  say  10  degrees.   Under 
favorable  circumstances,  estimates  to  within  30  degrees 
might  reasonable  be  expectedc   With  additional  study  and 
refinement,  the  30  degree  figure  might  be  reduced  signifi- 
cantly. The  biggest  advances  in  terms  of  orientation 
analysis  should  come  in  terms  of  getting  useful  results  in 
situations  that  are  currently  not  tractable  to  such  analysis, 
e.g. ,  1)  multiple  deformations  in  which  the  principal  axes 
of  stress  have  changed  orientation  by  45  degrees  or  more 
from  one  episode  of  deformation  to  the  next;  2)  where 
principal  axes  of  stress  have  acted  obliquely  with  respect 
to  an  initial  anisotropy,  the  orientation  of  which  cannot 
be  recovered  unambiguously;  3)  where  heterogeneous,  aniso- 
tropic material  has  been  acted  on  nonuniformly  over  the 
region  of  study.   At  this  time,  there  does  not  appear  to  be 
justification  for  estimating  stress  magnitudes  from  fabric 
analysis  of  soils  to  any  better  than  one  order  of  magnitude. 
With  additional  work,  it  should  be  feasible  to  improve  this 
to  a  half  or  quarter  order  of  magnitude  for  relatively 
uncomplicated  situations. 

To  work  out  the  stress  and  strain  history  of  soil 
from  fabric  analysis  requires  that  both  the  directional  and 
magnitude  estimates  described  be  feasible  and  that 
relative  ages  of  multiple  episodes  of  deformation  be 
established.   That  such  histories  can  be  derived  with 
confidence  in  some  cases  is  very  likely;  it  is  just  as  likely 
that  it  will  not  be  feasible  to  derive  such  histories  in 
other  cases,  at  least  in  the  foreseeable  future.   Studies 
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following  the  initial  inquiries  presented  in  this  report 
should  be  concerned  with  arriving  at  precise  statements  of 
the  conditions  and  combinations  of  conditions  that  will 
permit  the  derivation  of  reliable  stress  and  strain  histories, 
and  those  conditions  and  combinations  of  conditions  that 
necessarily  preclude  the  derivation  of  such  histories. 
Follow-up  studies  should  also  estimate  the  proportions  of 
tractable,  doubtful,  and  intractable  situations  , 
in  addition  to  advancing  the  state-of-the-art. 

Relative  scales  of  heterogeneities  in  the  earth 
materials,  of  fabric  elements  being  studied   and  of  the 
sample  size  are  of  considerable  importance  in  arriving  at 
valid  interpretations  of  fabric  (Figures  3-4  and  3-5) . 
Absolute  scales  are  of  considerably  less  significance  in 
fabric  analysis,  per  se,  although  they  may  be  crucial  in  a 
practical  design  problem  that  requires  fabric  analysis. 

Statistical  analysis  is  used  in  both  kinematic  and 
dynamic  approaches .   The  kind  of  outcome  to  be  expected 
depends  in  part  on  the  statistical  methods  used.   Descriptive 
statistics  yield  descriptive  results;  probability  statements 
result  from  procedures  such  as  hypothesis -testing.   Proba- 
bility statements  may  be  very  important  in  interpreting 
poorly  defined  fabrics,  fabrics  resulting  from  the  super- 
position of  one  fabric  on  another,  and  suites  of  apparently 
genetically  related  fabrics.   The  kind  of  interpretation  that 
can  be  made  will  also  vary  depending  on  whether  one  is  using 
parametric  or  nonparametric  methods,  Gaussian  or  von  Mises 
probability,  continuous  or  discrete  probability  density 
functions,  small-sample  or  large-sample  theory.   Although 
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the  statistical  armament arium  is  reassuringly  formidable, 
simple  descriptive  statements  that  can  be  rendered 
quantitatively  may  be  quite  adequate  for  determining  what 
one  needs  to  know  about  the  stress  field. 

The  type  of  earth  materials  being  studied  can  have 
considerable  impact  upon  the  type  of  fabric  analysis  per- 
formed and  its  relative  usefulness.   Figures  3-6  and  3-7 
illustrate  some  of  the  variety  of  planar  (surface)  and 
linear  fabrics  that  may  be  encountered. 

A  highly  oriented  depositional  fabric  may  be  quite 
insensitive  to  a  stress  field  oriented  with  maximum  com- 
pression normal  to  depositional  surfaces;  a  loosely  compacted 
sand  consisting  of  well-rounded,  spheroidal  particles  may 
yield  very  few  clues,  through  fabric  analysis,  of  its 
deformation  history;  homogeneous  plasmic  material  may  reveal 
preferred  orientation  patterns  only  through  x-ray  methods, 
with  little  or  no  noncrystallographic  fabric  data  in 
support.   Further,  a  material  in  which  two  or  more  sub  fabrics 
do  not  agree  in  symmetry  (heterotactic  fabric)  may  present 
insurmountable  interpretative  problems. 

Either  kinematic  or  dynamic  fabric  analysis  may 
yield  results  consistent  with  an  independently  evaluated 
in- situ  stress  field;  then  again  there  might  be  no  clear 
relationship.   Should  such  consistency  be  demonstrable  in 
a  particular  situation,  then  additional  fabric  studies  might 
be  useful  in  interpolating  or  perhaps  cautiously  extrapo- 
lating from  sites  of  verified  results.   In-situ  studies, 
including  such  techniques  as  x-ray,  acoustic  and  photoelastic 
coatings,  might  be  adapted  from  studies  of  rocks  in  the 
laboratory  and  in  the  field.   The  study  of  crack  generation 
and  propagation  in  stress-relieved  soils  in-situ  might  yield 
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very  useful  results  if  suitable  precautions  are  taken 

to  deal  with  anomalies  introduced  by  drying  and  excavation. 

Some  of  the  uncertainties  expressed  or  implied  in 
the  foregoing  statements  can  very  likely  be  dispelled  through 
experimentation  with  soils.   That  is,  it  should  be  feasible 
to  determine  whether  or  not  some  of  the  capabilities  described 
as  possibilities  can  in  fact  be  achieved  in  a  practical  way. 

3.3   STATISTICAL  METHODS 

As  in  any  investigation,  the  appropriateness  and 
efficiency  of  statistical  analysis  depend  in  part  on 
identification  of  the  variables  to  be  treated  and  the 
understanding  of  their  properties.   The  particular  tools 
chosen  will  depend  on  whether  a  variable  is,  for  example, 
continuous  or  discrete,  linear  or  angular,  infinite  or 
finite.   The  design  of  the  data-collection  program  and  the 
selection  of  analytical  methods  cannot  be  independent  of 
these  variables  and  their  statistical  properties.  Both 
kinematic  and  dynamic  fabric  analyses  are  served  by  such 
statistical  analyses.   The  variables  of  interest  in  relating 
soil  fabrics  to  strain  and  stress  history  and  state  of 
stress  can  be  divided  into  categories  such  as  direction, 
spacing,  and  pattern. 

Direction  applies  to  entities  such  as  orientation 
of  long  axes  of  peds,  grains  and  voids,  of  twin  lamellae 
and  fractures,  and  of  laminations  and  other  planar  or  linear 
expressions  of  composites.   Such  variables  are  angular  or 
cyclic.   Degree  of  preferred  orientation  (concentration  of 
direction)  and  direction  of  preferred  orientation  (mean 
direction)  may  be  determined.   Tests  to  ascertain  whether 
there  is  a  preferred  orientation  may  be  applied. 
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Spacing  applies  to  entities  such  as  size  of  grains, 
peds,  clods,  glaebules,  fragments,  and  voids,  and  width  of 
fracture  openings  and  tabular  fabric  elements.   Spatial 
frequency  (density  of  fractures  and  other  discontinuities), 
size  distributions,  and  elongation  ratios  (long: short 
dimension)  may  be  determined  from  spacing  data.   Tests  to 
ascertain  if  there  is  a  concentration  of  spacings  (or  their 
inverse,  spatial  frequencies)  may  be  applied. 

Pattern  applies  to  the  distribution  of  fabric 
elements  and  combinations  of  fabric  elements  into  features 
such  as  laminations,  bands,  or  layering,  smears,  clusters, 
and  tertiary  structure  of  secondary  peds.   In  turn,  spatial 
frequency,  size,  and  orientation  of  such  composite  features 
may  be  the  subject  of  study.   Spacing  and  direction  of 
constituent  fabric  elements  (subfabrics)  may  be  determined 
in  relation  to  each  other  within  and  between  composite 
features. 

The  kinds  of  directional  data  described  are 
examples  of  vectorial  fabric.   The  kinds  of  spacing  data 
described  are  examples  of  scalar  fabric;  under  the  heading 
of  scalar  fabric  come  also  such  diverse  attributes  as  size 
of  fabric  domains,  porosity,  and  composition.   Pattern  may 
include  both  vectorial  and  scalar  fabric. 

In  treating  vectorial  data,  the  von  Mises  probability 
density  functions  are  preferable  to  the  Gaussian  where  the 
range  of  the  data  occupy  a  substantial  portion  of  the  period 
of  the  variable.   Many  orientation  variables  in  the  earth 
sciences  are  characterized  by  such  large  dispersions. 

The  density  function  of  the  circular  normal  (von 
Mises)  distribution  is  given  by 
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k  cos(a-a  ) 

Y(a;a  .k)  -  S & 

8        2jtIo  (k) 


where 


a  =  angle  variate 

a  =  angle  coordinate  of  center  of  gravity 

o 

of  distribution. 
a  =  angle  coordinate  of  center  of  gravity 
of  sample  taken  as  maximum  likelihood 

estimate  of  a  •   a  is  calculated  from 

n  &        n 
tan  a  =   (Z   sin  a  )/(Z   cos  a  ) 

i=l     i   i=l     i 

Io(k)  =  a  Bessel  function  of  the  first  kind  of 
pure  imaginary  argument . 
k   =   a  measure  of  concentration  about  the 
mean;  determined  uniquely  from  a,  the 
vector  strength.   a  is  calculated  from 

a  =  R/n    for  n  measurements,  where 

n  n        2 

R2  -   (Z   sin  a   )2     +      (l   , cos  a> 
i-1     i        1=1    ± 

The  density  function  of  the  semicircular  normal  distribution 
is  given  by 

cos  2 (a  -  a J 

7  (a;  a  .k)   =  S S_ 

6         4ttIo  ( k) 

A  physical  model  of  the  circular  normal  distribution 
would  be  a  tilted  roulette  wheel.  The  circular  normal 
function  applies  to  variables  ranging  through  360  degrees, 
such  as  azimuth;  the  semicircular  normal  function  applies 
to  variables  ranging  through  180  degrees,  such  as  strike 
or  any  compass  orientation  with  direction  but  without  sense. 
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The  density  function  of  the  spherical  normal 
distribution  is  given  by 


7  (9,  a,  9g,  ag,  k) 


K 


4  irsinK    L 


exp  k[cos9cos9  +  sin9sin9  cos(a-a  )1 


where  k  satisfies  a   =  -1/&     +  coth  £. 

The  density  function  of  the  hemispherical  normal 
distribution  is  given  by 


y  (9,a;9g,ag,K)      = 


K 


K 


exp 


K[cos9cos9     +  sin9sin9_  c 


Ke 


e  -1 


g 
-1 


os(a-a  )] 


4  ir(e  -1) 
where  < satisfies  a  =   1/k 

(9  and  a  are  orientation  angles) 

8  =  angle  between  unit  vector  and  z-axis 

a  =  angle  between  projection  of  unit  vector  in 
xy-plane  and  x  axis . 

6  ,a   =  angle  coordinates  of  center  of  gravity  of 
distribution 

a  -  vector  strength 

Maximum  likelihood  estimates  of  parameters  are  denoted  by 
carets . 

The  orientation  of  planes  and  lines  in  the  earth  sciences 
is  commonly  depicted  with  reference  to  a  hemisphere  (dip 
and  strike  of  planes;  plunge  and  bearing  of  lines); 
accordingly,  the  hemispherical  function  above  is 
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of  particular  interest.   A  review  of  the  four  functions 

and  simplified  probability  tables  for  their  use  are  presented 

in  Pincus  (1953) . 

A   variety  of  tests  for  preferred  orientation  (or  non- 
uniformity)  is  available,  for  example,  including  those  in 
which  the  alternative  to  uniformity  is  characterized  by  a 
specified  density  function  and  those  in  which  the  alternative's 
characteristics  are  not  specified.   Both  two-and  three- 
dimensional  treatments  are  available;  two-dimensional  treat- 
ments in  mutually  perpendicular  planes  may  be  used  to 
approximate  three-dimensional  treatments  under  some  conditions 
and  where  inherent,  directionally-dependent  sampling  biases 
are  manageable.   In  studying  vectorial  soil  fabrics, 
particularly  in-situ,  two-dimensional  treatments  (both 
singly  and  multiply)  are  likely  to  be  considerably  more 
practical  than  three-dimensional  treatments. 

The  statistical  treatment  of  spacings  is 
essentially  the  same  as  that  applying  to  sizes.   Logarithmic 
scales  or  transformations  often  help  to  normalize  the  data 
so  that  Gaussian  and  other  commonly  used  probability  opera- 
tions can  be  undertaken  quite  easily.   In  dealing  with 
elongations,  the  statistics  of  ratios  comes  into  play,  and 
here  extreme  caution  is  necessary.  The  need  for  this 
caution  is  compounded  when  ratios  of  elongations  are  con- 
sidered, as  for  example  in  studying  relative  changes  in 
elongation  of  fabric  elements,  such  as  grains,  peds ,  or 
voids,  between  materials  having  been  subjected  to  different 
loads . 

The  density  function  of  the  ratio  z  of  two 
independent  and  normally  distributed  variates  x  and  y 
(z  =  y/x)  is  given  by 

f(z)   =   (1/  tt)  [l/  (1  +   z2)  ] 
which  has  neither  a  mean  nor  a  variance. 
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The  statistical  treatment  of  patterns  is  usually  a 
composite  of  vector  and  scalar  treatments,  the  exact  nature 
of  the  treatment  depending  on  the  features  involved.   A  soil 
consisting  of  micaceous  layers  alternating  with  granular 
quartzose  material  could  involve  vectorial  analysis  of  the 
orientation  of  the  mica  flakes,  scalar  analysis  of  the  size 
distributions  of  the  quartz  grains  and  lamination  thicknesses, 
and  vectorial  analysis  of  intra-,  inter-  and  trans-lamina 
fracture  orientations. 

Standardized  graphical  procedures  for  analysis  of 
fabric  data  have  been  available  for  many  years.   Figure  3-8 
shows  the  generation  of  the  most  commonly  used  projections 
for  displaying  and  graphically  processing  fabric  data. 
Figures  3-9  and^  3-10  show  graphical  treatment  of  such  data. 
Figures  3-11  and  3-12  display  specimen  plots  of  fabric  data 
illustrating  the  four  classes  of  symmetry  presented  earlier 
in  Figure  3-3. 

The  fact  that  many  earth-science  data  are  both 
multimodal  and  "noisy"  (erratically  dispersed)  places 
great  obstacles  in  the  way  of  using  many  conventional 
efficient  statistical  methods.   [Compare  Figure  3-12(c) 
which  is  multimodal  and  noisy  to  Figure  3-11  (b)J.  Sometimes 
one  must  settle  for  the  less  persuasive  results  of  non- 
parametric  methods.   Attempts  to  mitigate  these  undesirable 
qualities  by  increasing  sample  size  might  not  improve  the 
data  and,  in  fact,  might  introduce  even  more  complications  by 
including  part  of  yet  another  population.   Generally  speaking, 
whenever  the  data  have  more  than  one  mode  and  are  noisy,  the 
simpler  the  Statistical  method  chosen  (such  as  picking  modes 
instead  of  using  significance  tests  to  select  "concentrations") , 
the  more  easily  can  one  use  the  results,  particularly  in 
communicating  with  the  non-statistically-minded  user. 
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Figure  3-12.  Contoured  equal 
area  projections  illustrating 
monoclinic  sub fabrics  (a  andb) 
and  triclinic  fabrics  (c  andd); 
m  =  planes  of  symmetry 
(Paterson  and  Weiss,  1961). 
(Compare  to  Figure  3-3,  d  and  e) 


Figure  3-13.  Idealized  shear 
and  extension  fractures  above 
and  below  neutral  surface  at 
the  crest  of  a  fold,  and  above 
and  below  their  fabric  diagram; 
m  =  planes  of  symmetry 
(Friedman  and  Sowers,  1970). 
Fabric  diagrams  show  orienta- 
tions of  normals  to  fractures. 
Diagrams  show  parallel  ortho- 
rhombic  symmetry  above  and 
below  neutral  surface,  but 
state  of  stress  is  different 
in  each  domain. 
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Of  course,  what  appears  to  be  multimodal  might, 
on  closer  examination,  turn  out  to  be  two,  noncoincident 
unimodal  subfabrics,  e.g.,  two  nonparallel  sets  of  extension 
fractures  (say  one  from  loading  and  another  from  unloading) , 
or  a  combination  of  generally  related  modes  (Figure  3-13) . 
A  multimodal  fabric  may  be  replaceable  by  a  unimodal 
fabric  by  virtue  of  a  presumably  valid  physical  relation, 
e.g.,  substitution  of  the  acute  bisectrix  (axis  of  maximum 
compression)  for  conjugate  shear  pairs.   The  point  here  is 
that  the  data  may  be  simplified  via  nonstatistical 
considerations  and,  in  fact,  they  should  be  simplified  as 
much  as  possible  on  physical  premises  prior  to  statistical 
treatment . 

A  straightforward  analog  method  for  handling 
vectorial  and  scalar  data  simultaneously  is  through 
spatial  Fourier  analysis  by  optical  diffraction.   Figures 
3-14  to  3-24  illustrate  the  basic  procedures  and  lay-outs  in 
optical  diffraction  analysis  that  have  been  applied 
successfully  to  fabric  studies  in  rock  mechanics.   Both 
orientation  and  spacing  data  are  transformed  directly  into 
a  polar  coordinate  diagram  (diffraction  pattern  or  two- 
dimensional  Fourier  transform)  from  the  input,  which  is 
a  transparent  replica  (usually  reduced)  of  the  surface  of 
the  earth  material . 

The  two-dimensional  Fourier  transform  is  given  by 

■ 
F   (u,v)   =   //f(x,y)  exp   -2iri(ux  +  vy)  Af  dxdy 

where  f (x,y)   =   input  function  as  a  function  of 
position  coordinates  x,y  in  the  object 
plane  (front  focal  plane) 

F(u,v)   =   complex  amplitude  as  a  function 
of  position  coordinates  u,v  in  the 
transform  plane  (back  focal  plane) 

f  =   focal  length  of  lens 

A   =  wavelength  of  incident  light 
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Figure  3-15.  Relation  of  diffraction  pattern  to  distribution 
parameters  (Dobrin,  et  all,  1964). 

r  is  proportional  to  spatial  frequency 

9  is  direction  normal  to  the  lineation  giving  rise  to 
diffraction  pattern  AA' 
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A)  NARROWER 
SLIT   SPACING 

($,) 


B)WIDER     SLIT 
SPACING     (SJ 


r 

FOCAL 
PLANE 

(TRANSFORM 
OR  SPECTRUM 
PLANE) 


ZERO  ORDER 
FIRST  ORDER    (UP) 

FIRST   ORDER    (DOWN) 


(i 


A- FIRST  ORDER    POINT    (UP) 

B- ZERO  ORDER    POINT 

C- FIRST  ORDER     POINT    (DOWN) 


f 


=    V2 
=    2Y, 


Figure  3-18.  Spatial  frequency  analysis  by  diffraction 
via  lens   (After  Shulman,  1970). 
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FOURIER  TRANSFORM  PAIRS 
FOR  TYPICAL  SIGNALS 

Figure  3-19  Fourier  Transforms  of  simple  functions 
(Peterson  and  Dobrin,  1966) 

The  functions  in  each  column  are  the  transforms  of 
the  corresponding  functions  in  the  adjoining  column. 

A  white  line  on  a  black  background  would  be  represented 
by  b),  the  square  wave. 

The  diffraction  pattern  that  is  seen,  photographed,  or 
mapped  is  actually  the  square  of  the  transform  function, 
so  all  negative  values  in  the  transform  are  rendered 
positive. 
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Figure  3-20.   Fourier  Transforms  of  truncated  functions 
(Peterson  and  Dobrin,  1966). 

Functions  such  as  those  in  3-19  are  truncated  with  a 
window  function.   The  convolution  process,  denoted  by 
*,  explains  the  derivation  of  the  transform  of  the 
truncated  function  from  the  infinite  function  and  the 
window  function. 
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Figure  3-24  Filtering  of  input  of  polar  coordinate  paper 
Inputs  and  outputs  are  negative  photographs. 
(Pincus,  1969a). 

A  -  Unfiltered  input  (0°  reference  at  6  o'clock) 

B  -  Transform  of  A 

C  -  Directionally  filtered  image  (10°  wedge) 

D  -  Filtering  for  C  (filter  centered  at  90°  and  270°) 

E  -  Directionally  filtered  image  (40°  wedge) 

F  -  Filtering  for  E  (filter  centered  at  135°  and  315°) 

G  -  High  cut  filtered  image 

H  -  Filtering  for  G  (circular  filter) 
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This  equation  expresses  the  relation  between  the 
spatial  distribution  of  light  passing  through  the  input  image 
in  the  focal  plane  and  the  distribution  in  the  back  focal 
plane  where  the  diffraction  pattern  (Fourier  transform)  is 
formed. 

Spatial  filtering  is  rendered  practical  by  the  fact 
that  the  Fourier  transform  of  the  Fourier  transform  is  the 
input  function,  i.e.,  the  output  function  is  the  reconstructed 
image  of  the  input  image.   Such  filtering  is  accomplished,  in 
its  simplest  form,  by  placing  filters  (masks)  in  the  transform 
plane  and  an  additional  transform  lens  (reconstruction  lens) 
one  focal  length  beyond  the  transform  plane. 

Through  directional  and  frequency  filtering,  some 
features  can  be  suppressed  or  eliminated  in  the  reconstructed 
image  (spatially  filtered  image)  so  that  less  prominent 
details  of  the  fabric  can  be  studied  more  easily. 

The  optical  diffraction  method  has  been  used  to 
study  changes  in  fabric  of  rocks  under  load,  and  has  equally 
great  promise  for  study  of  soils.   Both  the  overall  style 
of  fabric  and  changes  or  contrasts  in  fabric  can  be 
quantified  by  this  method.   Features  of  fabric  at  widely 
different  scales  can  be  studied  making  it  feasible,  for 
example,  to  compare  soil  fabric  viewed  by  the  scanning 
electron  microscope  with  that  viewed  by  the  light  micro- 
scope, in  hand  specimens  and  in-situ  (in  the  excavation). 

Optical  diffraction  analysis  is  further  illustrated 
and  applied  later  in  the  section  titled  Application  of 
Optical  Diffraction  Analysis  to  Soil  Fabrics.   Here,  trans- 
form (diffraction  pattern)  analysis  is  illustrated;  spatial 
filtering  of  soil  fabric  data  is  not  illustrated  in  the 
current  report. 

Results  equivalent  to  those  of  optical  diffraction 

analysis  can  also  be  obtained  digitally,  but  less  conveniently 

and  with  less  versatility  in  varying  data-processing  para- 
meters.  Digital  processing  can  be  very  helpful,  however,  in 
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preparing  and  analyzing  equal-area  diagrams  obtained  from 
universal-stage  work  and  field  measurements. 

The  analysis  of  x-ray  data  obtained  from  soils  can 
yield  useful  information  on  crys tallographic  orientations. 
In  a  sense,  such  data  are  already  processed  statistically 
by  the  very  nature  of  the  procedures  used  to  obtain  the  x- 
ray  data. 

3.4  ANTICIPATED  BEHAVIOR  OF  VARIOUS  SOIL  TYPES 

The  change  in  fabric  associated  with  a  particular 
loading  history  depends  on  factors  such  as  the  direction, 
intensity,  and  rate  of  loading  and  the  constitution  of  the 
material  that  has  been  loaded.   Included  in  these  constitu- 
tional factors  are  scalar  and  vector  fabrics,  as  well  as 
(or  overlapping  with)  the  usual  compositional  expressions 
of  soil  mechanics. 

Here  it  is  important  that  the  primary  fabric 
of  the  soil,  i.e.,  the  fabric  developed  during  formation 
of  the  soil  and  prior  to  subsequent  deformation,  be 
known  or  accurately  inferred  in  order  to  characterize  the 
deformation  in  some  useful  way.   Whether  the  soil  initially 
has,  for  example,  single-grained,  honeycombed,  or  flocculent 
structure  will  determine,  in  part,  what  kinds  of  inter- 
pretations can  be  made  following  deformation. 

In  soils  with  single-grained  structure,  i.e., 
cohesionless  soils  such  as  sands  and  gravels,  individual 
grains  may  be  displaced  or  rotated  to  save  space  in  the 
direction  of  maximum  compression,  grains  may  be  crushed 
thereby  increasing  the  dispersion  of  particle  sizes  and 
providing  a  means  for  packing  intergranular  voids  with  fine 
material,  and  individual  grains  may  be  deformed  elastically, 
plastically,  or  by  fracture.   Elastic  deformation  of 
individual  grains  or  aggregates  of  grains  might  be  deter- 
mined in-situ  by  photoelastic  or  x-ray  methods  adapted 
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for  this  purpose;  plastic  deformation  might  be  detected 
by  studies  of  grain  elongation  or  twin-gliding.   The 
deformation  of  intergranular  voids  in  soils  with  single- 
grained  structure  requires  essentially  the  same  treatment 
as  in  the  study  of  changes  in  grain  shape  and  size. 

In  soils  with  honeycombed  structure,  viz.  fine 
enough  for  cohesion  to  be  significant  such  as  in  fine  silts 
and  clays,  the  undisturbed  honeycomb  structure  itself  could 
provide  a  datum  against  which  to  measure  deformation.   Of 
course,  a  pervasive  destruction  of  the  honeycombs  might 
leave  no  evidence  of  the  initial  structure.   Individual 
grains  might  undergo  the  same  changes  as  considered 
for  grains  in  single-grained  structure.   With  cohesion, 
structural  units  of  higher  organization  such  as  peds  and 
clods  can  provide  yet  another  scale  at  which  deformation 
may  be  characterized. 

In  soils  with  flocculent  structure,  viz.,  those 
consisting  of  colloidal-size  particles,  measures  of  deforma- 
tion may  be  based  on  changes  in  grains,  floes  (honeycombs 
of  grains)  and  honeycombs  of  floes  (second-order  honeycombs) . 
In  such  fine-grained  material,  the  practical  approach  may  be 
to  study  changes  in  larger-scale  features  such  as  peds, 
fractures,  and  bands. 

Soils  of  mixed  structures  provide  a  larger  variety 
of  features  with  which  to  attempt  kinematic  and  dynamic 
fabric  analysis,  but  this  also  makes  interpretation  much 
more  difficult.   The  presence  of  more  than  one  type  of 
structure,  unless  the  structures  are  homogeneously  distributed 
through  the  material,  creates  interpretation  problems  of  as 
great  difficulty  as  those  arising  from  a  history  of  more  than 
one  episode  of  deformation. 

The  presence  of  water  (or  its  absence)  certainly 
affects  the  course  of  deformation  and  consequently  the  fabric 
associated  with  deformation.   For  example,  the  total 
absence  of  water  in  granular  noncohesive  soils  is  likely 
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to  yield  obviously  cataclastic  features  with  deformation. 
The  role  of  water  is  increasingly  dramatic  with  decreasing 
particle  size;  in  fact  it  can  be  asserted  that  the  same 
fine-grained  soil,  wet  and  dry,  is  in  fact  not  the  same 
material  under  these  two  sets  of  conditions,  even  with  no 
changes  in  soil  mineralogy.   Despite  all  these  contrasts, 
it  is  not  likely  that  the  presence  or  absence  of  water  will 
enhance  or  reduce  the  usefulness  of  fabric  analysis.   And, 
experiments  to  determine  unambiguously  the  role  of  water  in 
fabric  change  should  be  limited,  at  least  initially,  to  those 
ranges  of  saturation  most  likely  to  be  found  in  the  major 
soil  types. 

Table  3.1  shows  possible  deformation  indicators 
that  could  be  used  for  soil  types  recognized  in  the  Unified 
Soil  Classification  Systems.   Table  3.2  shows  a  comparable 
breakdown  by  type  of  fabric  component  in  soils  as  recognized 
in  pedology  (Brewer,  1964).   Both  the  Unified  Soil  Classi- 
fication System  and  Brewer's  breakdown  of  soil  components 
appear  to  lend  themselves  to  useful  associations  with 
deformation  indicators. 

An  important  step  in  follow-up  studies  will  be  to 
relate,  where  feasible,  features  of  soil  fabric  to  basic 
soil  mechanics  properties  such  as  sensitivity  and  Atterberg 
limits,  and  to  contrasts  in  properties  such  as  those  between 
remolded  and  unremolded  soil.   Microfabric  features,  such  as 
average  space  between  particles  (pore  fabric) ,  might  explain 
transitional  points  such  as  liquid  limits  in  rational,  use- 
ful terms.   Figure  3-25  indicates  some  interesting 
possibilities . 

3.5   ORIENTATION  OF  MEASURING  DEVICES  WITH  RESPECT 
TO  SOIL  FABRIC 

The  design  of  ins trumentational  configurations  for 
measuring  stress  or  strain  in  earth  materials  is  usually 
rendered  difficult  by  the  materials '  anisotropy  and 
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(a) 


(b) 


(c) 


Figure  3-25.  Alignment  of  grains  in  experimentally 

deformed  clays  (Turner  and  Weiss,  1963). 

a)  Rolling  in  direction  and  sense  R  about  axis  A 

b)  Flattening  of  a  sphere  to  a  disk  normal  to  axis  F 

c)  Extrusion  of  a  cylinder  parallel  to  axis  E  . 
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heterogeneity,  the  presence  of  discontinuities  (both  natural 
and  man-made) ,  the  perturbations  in  behavior  of  the  earth 
materials  induced  by  the  instrumentation  itself,  the  frequent 
nonparallelism  between  principal  stress  and  principal 
directions  of  mechanical  anisotropy,  and  departures  from 
ideal  elastic  and  plastic  behavior  of  many  earth  materials 
even  where  intact,  homogeneous,  and  approximately  isotropic 
materials  are  involved.   Interpretations  of  the  results  of 
such  measurements  are  rendered  difficult  by  the  same 
factors . 

Generally  these  difficulties  are  encountered 
whether  one  is  measuring  (a)  stress  or  strain  resulting 
from  a  known  external  load  or  (b)  stress  or  strain 
associated  with  residual  stress.   The  latter  measurements 
are  more  subject  to  error  than  are  the  former.   Strain 
measurements  are  typically  more  reliable  than  stress  measure- 
ments (at  least  some  of  which  are  not  really  measurements  of 
stress)  and  the  strain  field  can  probably  be  more  easily 
related  to  fabric  anisotropy. 

As  mechanical  anisotropy  increases,  these 
difficulties  in  design  and  interpretation  may  be 
accentuated.   Where  fabric  anisotropy  is  related  to 
mechanical  anisotropy  in  some  direct,  systematic  way, 
fabric  analysis  can  be  of  help  in  anticipating  design  and 
interpretation  problems.   In  addition,  fabric  analysis 
might  help  directly  in  specifying  orientation  of  some 
components  of  the  measuring  device  and  in  calculating 
meaningful  results.   For  example,  if  fabric  analysis  reveals 
a  single,  strong,  preferred  orientation  with  axial  symmetry 
(Figures  3-3a  and  3-lla)  ,  (e.g.,  parallel  channel-voids), 
which  is  probably  associated  with  a  strong,  similarly  oriented 
mechanical  anisotropy,  then  one  could  orient  one  of  the 
measuring  axes  parallel  to  this  direction  in  order  to 
increase  the  sensitivity  of  the  entire  measuring  system. 
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On  the  other  hand,  should  fabric  analysis  reveal 
that  the  fabric  is  heterotactic ,  i.e.,  the  symmetries  of  the 
constituent  sub fabrics  do  not  coincide,  then  the  most  practical 
strategy  might  be  to  ignore  the  fabric  in  orienting  the 
instrumentation  but  to  recall  the  details  of  the  fabric  in 
interpreting  the  results. 

All  of  the  foregoing  remarks  apply  as  well  to  time- 
dependent  and  location-dependent  measurements. 

Considering  the  complexities  outlined  in  the  first 
paragraph  of  this  section,  it  seems  wise  to  investigate 
possibilities  for  making  measurements  initially  in  many  or 
all  directions  in  a  given  plane,  to  do  this  in  three  mutually 
perpendicular  directions,  and  to  use  fabric  analysis  to 
help  explain  the  results.   From  this  might  come  sufficient 
understanding  to  use  fabric  data  to  specify  orientations  of 
measuring  devices.   Photoelastic  borehole  inclusions, 
photoelastic  coatings,  brittle  lacquer  coatings,  and  Moire 
coatings  should  all  be  considered  for  this  purpose.  The 
coating  and  its  adhesive  medium  should  be  selected  and  applied 
so  as  not  to  result  in  a  significant  reinforcing  effect. 
Coatings  with  thicknesses  less  than  several  times  that  of 
fabric  elements  such  as  grains  may  respond  to  grain-by-grain 
strain  gradients;  accordingly  the  thickness  must  be  considered 
as  part  of  the  instrument  design.   Photoelastic  borehole 
inclusions  can  be  fabricated  from  a  variety  of  materials  that 
can  be  used  to  achieve  desired  matches  or  mismatches  between 
mechanical  properties  of  the  inclusion  and  the  surrounding 
earth  materials. 

Despite  the  importance  of  fabric  anisotropy  in  all 
of  these  conisderations ,  it  seems  likely  that  the  local 
vertical  will  usually  be  one  of  the  directions  of  great 
interest,  and  at  least  one  direction  in  the  horizontal  will 
be  of  at  least  secondary  interest.   (In  fact,  both  local 
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vertical  and  some  direction  or  directions  in  the  horizontal 
are  likely  to  be  important  fabric  directions  in  some  and 
perhaps  many  cases).   Where  the  fabric  is  unknown  or  diverse 
(e.g.,  heterotactic) ,  vertical  and  horizontal  orientations 
may  be  included  in  the  measuring  program.   Where  there  are 
also  conspicuous  discontinuities  of  moderate  dip,  such  as 
joint  sets  (Figure  3-13)  or  layering,  a  third  (or  fourth) 
orientation  could  be  taken  perpendicular  to  these  surfaces. 

The  design  of  the  instrumental  configuration  should 
include  provision  for  measuring  as  many  variables  as  would  be 
required  in  the  worst  case  possible  in  the  situation  being 
studied.   As  additional  information  about  anisotropy  and 
loading  at  the  locality  are  obtained,  the  number  of  measure- 
ments required  may  be  decreased  from  the  original  theoretical 
maximum.   With  increasing  experience,  additional  reductions 
may  be  effected. 

3.6   EXPERIMENTAL  WORK  TO  INVESTIGATE  RELATIONS  AMONG 
CHANGES  IN  FABRIC,  STRESS  HISTORY,  AND  INITIAL 
FABRIC  AND  COMPOSITION 

A  principal  objective  of  experimental  work 
should  be  to  determine  under  what  conditions' and  with 
which  soil  types  dynamic  fabric  analysis  is  feasible,  the 
combinations  of  conditions  and  materials  permitting  useful 
kinematic  fabric  analysis,  and  the  combinations  of  conditions 
and  materials  that  will  not  permit  effective  use  of  any  kind 
of  fabric  analysis. 

Materials  selected  should  represent,  at  the  very 
least,  single-grained,  honeycombed  and  flocculent  initial 
structures.   Preferably,  the  selection  should  give 
representation  to  perhaps  half  to  two-thirds  of  the  15 
categories  in  the  Unified  Soil  Classification.   Beyond  this, 
specimens  should  also  be  selected  to  show  effects  on  larger- 
scale  features  such  as  peds  and  clusters,  and  to  oriented 
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features  such  as  elongated  grains,  layering  and  smears.   As 
noted  elsewhere  (Section  3.4  ANTICIPATED  BEHAVIOR  OF  VARIOUS 
SOIL  TYPES) ,  variations  in  saturation  should  be  limited  to 
ranges  to  be  anticipated  in-situ  in  the  particular  type  of 
soil  being  studied.   Grossly  heterogeneous  materials  should 
be  avoided,  at  least  initially  and  perhaps  entirely. 

Variables  in  loading  to  be  included  in  an  experi- 
mental program  should  be  magnitude,  rate,  ratio  of  principal 
stresses,  orientation  of  principal  stresses  with  respect 
to  orientation  of  directional  fabric  elements,  and  pore 
pressure.   Both  normal  and  shear  loading,  in  relatively 
simple  configurations,  should  be  used.   Dynamic  testing, 
if  undertaken  at  all,  should  come  in  the  latter  stages  of 
an  experimental  testing  program. 

Another  objective  of  experimental  work  should  be 
to  help  lay  the  groundwork  for  in-situ  tests  that  could  be 
later  adapted  for  routine  use  by  the  field  practitioner. 
It  seems  advisable  to  begin  the  experimental  program 
very  simply,  with  a  small  number  of  pilot  experiments  designed 
to  indicate  which  categories  of  experiments  are  likely  to  be 
most  productive  and  should  be  explored  in  the  later  phases 
of  experimentation. 

Samples  of  well-graded  and  poorly  graded  sands  (SW 
and  SP) ,  saturated  and  dry,  tested  in  compression  and  in 
shear,  should  have  their  fabrics  recorded  before  and  after 
deformation.   Loads  should  be  carried  incrementally  to 
progressively  higher  values,  for  sets  of  like  specimens.   The 
fabrics  may  be  recorded  on  photographs  or  acetate  peels.   The 
analysis  of  the  fabrics  may  then  consist  of  conventional 
microcounting  and  micromeasuring  techniques,  opto-electronic 
counting  and  digital  processing,  or  optical  diffraction 
analysis.   Associated  with  whichever  of  these  techniques  may 
be  selected  should  always  be  old-fashioned  "eyeball"  evalua- 
tion overall  of  the  fabric. 
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Samples  of  inorganic  silts  and  clays  of  low  and  high 
liquid  limit  and  from  slight  to  high  plasticity  (ML,  CL,  MH, 
CH)  tested  in  compression  uniaxially  and  triaxially  and  in 
shear,  should  have  their  fabrics  recorded  before  and  after 
deformation  as  in  the  treatment  of  the  sands  above. 
Samples  of  organic  silts  and  clays  and  highly  organic 
soils  (OL,  OH,  P  )  should  also  be  treated  this  way.   As  the 
grain  size  decreases  and  the  organic  content  increases  there 
are  likely  to  be  practical  problems  in  using  acetate  peels 
to  record  the  fabrics  and  there  may  also  be  problems  in  macro 
and  microscale  photography.   Here  it  might  be  necessary  to 
turn  to  scanning  electron  microscopy  to  obtain  useful  fabric 
data. 

In  the  early  stages  of  experimentation,  samples 
without  obvious  stratification,  large  discontinuities  and 
other  conspicuous  heterogeneities  should  be  selected.   In 
addition  to  studying  effects  of  loading  on  fabric,  early 
experimentation  should  also  provide  for  evaluation,  adapta- 
tion, and  development  of  a  variety  of  experimental  techniques. 
For  example,  strain-birefringence  effects  in  quartz  grains, 
lattice-strain  determined  by  x-ray  diffraction,  and 
stereological  techniques  such  as  grain  intercept-area  ratios, 
should  be  considered  and  perhaps  attempted.   Partially 
indurated  soils,  in  which  natural  chemical  cements  or  ice 
serve  as  the  indurating  medium,  should  also  be  studied  in 
an  intermediate  phase  of  experimentation.   The  percentage 
induration  and  its  spatial  distribution  would  be  of  interest, 
particularly  in  frozen  soils. 

We  have  investigated  the  possibility  of  applying 
optical  diffraction  analysis  to  soil  fabrics  of  different  types 
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As  inputs,  14  photographs  from  Brewer  (1964),  reduced  four 
times,  have  been  used  to  generate  diffraction  patterns. 
The  results  are  presented  in  Figures  3-32  through  3-45. 
The  same  technique  can  be  applied  to  single  examples  at 
progressively  higher  loads  or  to  comparable  specimens 
loaded  in  different  modes  (uniaxial,  triaxial,  direct  shear). 

Details  of  optical  diffraction  analysis  and  inter- 
pretation are  covered  elsewhere  in  references  such  as  Pincus 
and  Dobrin  (1966),  Pincus  (1966,  1969(a),  1973),  and  Pincus, 
Power,  and  Woodzick  (1973) .   Results  of  spatial  filtering 
and  holographic  subtraction  as  applied  to  soil  fabrics  are 
not  presented  here,  but  will  likely  prove  to  be  powerful 
tools  in  subsequent  work. 

In  studying  inputs  and  their  transforms  or 
diffraction  patterns,  the  evaluation  of  scaling  factors 
is  necessary  if  the  transform  is  to  be  used  as  a  measuring 
device .   Typically,  a  small  piece  of  rectangular  coordinate 
paper  is  photographed  with  the  input,  and  is  then  used  to 
generate  a  calibration  transform.   Lines  intersecting  in  a 
rectangular  grid,  like  the  vertical  and  horizontal  lines 
in  the  input  of  Figure  3-17,  generate  series  of  dots  (along 

horizontal  and  vertical  lines  as  in  the  Figure  3-17  transform) 
also  distributed  in  a  rectangular  network.   The  spacing 
between  the  dots  varies  linearly  with  the  spatial  frequency 
of  the  input  lines.   Figure  3-31  is  the  calibration 
transform  for  Figures  3-32  to  3-45;  calibration  transforms 
also  appear  in  Figure  3-27 (b)  ,  3-49 (b)  and  3-50 (b) ; 
profiles  of  calibration  transforms  are  shown  in 
Figures  3-51(a)  and  3-52  (upper  left).  Table  3.3  summarizes 
scaling  relationships  for  inputs  and  outputs  in  this  work. 
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a)   Sphericity  and 
roundness  chart 
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b)  Assemblages  of 
particles  such 
as  those  in  a). 


c)   Transforms  of  b). 
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Figure  3-26.  Idealized  particulate  shapes  and  their  transforms 
(Larson,  1973). 
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Figure  3-29 


Some  reference  inputs  (line  drawings)  and 
their  respective  transforms  (Pincus,  1973). 
(Inputs  from  Jung,  1969). 
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Figure  3-31.  Calibration  Transform 

The  central  bright  cross  is  generated  by  the 
rectangular  aperture  used  to  illuminate  the  input  trans- 
parency. 

The  slightly  rotated  cross  composed  most  conspicuously 
of  bright,  circular  dots  is  the  calibration  transform, 
generated  by  the  millimeter  coordinate  paper  visible  in  the 
lower  right  corner  of  the  inputs  that  follow.   This  milli- 
meter paper  provided  the  scale  for  the  photographs;  the 
scale  for  the  specimens,  per  se,  can  be  obtained  from  the 
scale  bar  in  the  bottom  center  of  the  picture. 

The  upper  margin  of  the  transform  shows  the  interval 
corresponding  to  a  spatial  frequency  of  1  line/mm  in  the 
photograph  (lmm  spacing).   Twice  this  interval  represents 
2  lines /mm  (0.5mm  spacing)  in  the  photograph,  and  so  forth. 
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Figures  3-26  to  3-30  demonstrate  some  results  from 
earlier  studies  that  indicate  some  possibilities  for  analysis 
of  soil  fabrics.   Particle  shapes  (Figure  3-26)  and  shapes, 
sizes,  and  arrangements  of  grains  (Figures  3-27  to  3-29)  are 
all  amenable  to  this  type  of  analysis.   Spatial  filtering  of 
a  "real"  input  is  illustrated  in  Figure  3-30.   Figure  3-29 
is  taken  from  a  reference  collection  of  inputs  and  their 
transforms,  used  for  quick  comparisons  in  fabric  studies  of 
rocks . 

Figures  3-32  to  3-35  are  a  preliminary  sampling, 
with  inputs  from  Brewer  (1964)  ,  to  show  relationships 
between  a  variety  of  soil  fabric  elements  and  types 
of  fabrics  on  the  one  hand,  and  their  transforms,  (diffraction 
patterns)  on  the  other.   Input  captions  are  quoted  from 
Brewer. 

Input  transparencies  prepared  from  the  input 
photographs  were  prepared  on  35  mm  film  at  a  reduction  of 
4X.   The  first  six  input  transparencies  were  negatives  and 
the  remaining  eight  were  positives.   The  prints  of  the 
transforms  of  the  first  four  inputs  were  given  twice  as 
much  exposure  as  their  10  counterparts  in  order  to  bring 
out  useful  details . 

3.7   PLANS  FOR  MAPPING  AND  ANALYZING  SOIL  FABRICS 

Plans  for  mapping  and  analyzing  soil  fabrics 
must  always  take  into  account,  implicitly  if  not  explicitly, 
that  the  ultimate  general  objectives  are  to  determine  the 
orientation  of  principal  stress  directions  and  the  magnitudes 
of  principal  stresses  associated  with  the  fabric,  and  the 
stress  and  strain  history  of  the  soil. 

Where  some  features  of  soil  fabric  that  appear  to 
bear  upon  any  of  these  objectives  can  be  perceived  quickly, 
the  type  and  layout  of  mapping  can  be  tentatively  specified 
to  optimize  utilization  of  these  features.   For  example, 
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the  type  and  layout  of  mapping  will  not  be  the  same  in  one 
area  with  a  single  dominant  joint  direction  as  in  another 
area  with  three  mutually  perpendicular  joint  directions.   Nor 
will  the  mapping  and  analysis  procedures  be  the  same  in  an 
area  in  which  the  most  conspicuous  feature  is  a  single  set  of 
joints  as  in  another  area  showing  a  single  set  of  joints  and 
squeezed-out  peds,  indicating  compression  parallel  to  the 
direction  of  tension  inferred  from  the  joint  orientations 
alone.  Work  on  cross-cutting  relations  might  establish  the 
relative  ages  of  two  episodes  of  stressing  represented  by 
the  joints  and  the  elongation  of  the  peds. 

The  selection  of  type  of  mapping  of  soil  fabrics  , 
and  the  type  of  analysis  associated  with  such  mapping, 
depends  on  (a)  the  purpose  of  the  fabric  analysis  (and  the 
investigation  of  which  it  is  a  necessary  but  not  unique 
component)  (b)  the  type  of  soil  material  being  studied, 
(c)  the  form  and  size  of  the  body  composed  of  the  soil 
material  (such  as  a  mappable  unit  on  an  engineering  soil 
map),  (d)  the  spatial  distribution  of  anisotropies  in  the 
soil  body,  (e)  the  directional  distribution  of  anisotropies 
in  the  soil  body,  (f)  the  type,  form  and  requirements  of 
the  engineering  structure  to  be  designed  and  constructed, 
(g)  the  range  of  possible  orientation  relations  between 
the  engineering  structure  and  soil  anisotropies,  (h)  the 
range  of  possible  size  relations  among  the  engineering 
structure,  boundaries  of  soil  body,  and  soil  heterogeneities, 
(i)  the  scale  or  scales  at  which  fabric  is  to  be  mapped, 
and  (k)  the  resources  available  for  such  mapping.   All  of  the 
foregoing  factors  are  mutually  interdependent;  none  can  be 
considered  effectively  in  isolation. 

Mapping  techniques  may  range  from  conventional, 
infrared,  radar,  and  multispectral  photography,  down  to 
scanning  electron  microscopy  and  x-ray  diffraction.   It  is 
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desirable  that  soil  fabrics  be  studied  at  more  than  one 
scale  in  any  soil,  say  at  the  scales  of  grains  and  of  peds , 
laminations  and  large  joints,  grains  and  laminations,  and 
so  forth. 

At  the  scale  of  photographs  of  large  surface 
exposures,  the  recorded  differential  responses  of  different 
soils  to  incident  radiation  may  yield  only  scalar  fabric 
data;  although  vector  fabric  data  on  soils  may  be  recovered 
at  this  scale,  they  are  less  likely  to  be  as  available  as 
in  studies  of  comparable  rock  surfaces.  Multi-spectrai  photo- 
graphy of  large  surface  exposures  should  yield  additional 
useful  information  in  those  studies,  as  it  has  in  so  many 
remote  sensing  investigations. 

At  the  scale  of  small  surface  exposures,  say  those 
visible  in  test  pits,  shafts,  and  tunnels,  some  fabric 
elements  such  as  layering,  jointing  and  oriented  pebbles, 
can  be  recorded  quite  satisfactorily  by  careful  sketching; 
such  fabrics  can  be  recorded  more  precisely  by  good-quality 
photography  of  the  type  routinely  done  in  field  work  by 
earth  scientists  and  engineers  or  by  macrophotography. 
Figure  3-46  illustrates  some  of  the  types  of  features  that 
should  be  recorded. 

Cuxur  photography  is  not  universally  essential; 
in  fact  it  is  usually  easier  to  establish  standard  sets  of 
mutually  comparable  black-and-white  photographs  than  to 
attempt  this  with  color.   Photography  in  35  mm  and  120  formats 
is  probably  the  most  practical  in  this  work.   Voids  in  the 
exposed  surface  can  be  highlighted  by  impregnation  with 
florescent  material  and  illumination  with  ultraviolet 
radiation,  and  then  photographed  on  black-and-white  film.   A 
variety  of  impregnating  materials  for  both  dry  and  wet  soiis 
is  available  for  stabilization  as  well  as  for  florescence 
although  this  might  be  more  difficult  in-situ  than  in  the 
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Figure  3-46.  Surfaces  of  discontinuity  or  heterogeneity 

defined  by  preferred  orientation  and  preferred 
location  of  fabric  elements  (Turner 
and  Weiss,  1963) . 

a)  Preferred  orientation  alone 

b)  Preferred  location  alone 

c)  Combination  of  a)  and  b) 

d)  Combination  of  a)  and  b)  defines 
two  surfaces  oblique  to  each  other. 
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laboratory.   Recording  of  fabrics  in-situ  with  acetate  peels 
might  be  very  effective;  whether  or  not  prior  stabilization 
by  impregnation  will  commonly  be  necessary  is  not  known  at 
this  time . 

The  scale  of  small  surface  exposures  grades  into  that 
of  samples,  such  as  the  pedographer 's  soil  monolith  (,67m  x  ,15m 
x.lOm);  (26in.  x  5.9in.  x  3.9  in.)  and  so-called  undisturbed  core 
samples,  and  the  techniques  for  mapping  the  fabric  are  much  like 
those  to  accomplish  microscopic  work  in  the  laboratory,  where 
conditions  for  recording  the  fabric  can  be  more  carefully 
controlled  and  the  use  of  high-resolution  photographic  equip- 
ment is  more  feasibile.   Inherent  in  working  with  samples  is 
the  hazard  of  destroying  fabric  elements,  particularly  of 
the  vectorial  type,  in  sampling  and  subsequent  processing  of 
samples . 

Samples  used  for  microscopic  and  x-ray  analyses  may 
yield  valuable  information  at  the  intragranular  scale.   Some 
features  studied  at  this  scale  (fractures,  dislocations,  etc.) 
can  be  handled  much  the  way  they  are  in  larger  samples,  but 
some  of  the  techniques  used  (optical  and  x-ray  diffraction) 
yield  aggregated  data  directly.   The  preparation  of  thin 
sections  of  soils  can  be  accomplished  routinely,  but  it  is  not 
known  to  what  extent  the  preparation  of  thin  sections  will 
diminish  the  ability  to  resolve  slight  differences  in 
vectorial  fabrics. 

All  plans  for  mapping  fabrics  exposed  in  surfaces, 
other  than  at  aerial  scales,  should  provide  the  opportunity 
to  map  fabrics  on  three  mutually  perpendicular  surfaces 
bounding  presumably  homogeneous  materials.   This  might  not 
be  practical  to  achieve  in  test  pits,  shafts,  and  tunnels, 
but  it  should  be  relatively  easy  to  achieve  down  to  the  scale 
of  the  thin  section.   As  the  specimens  become  smaller  and 
approach  a  size  only  a  few  times  larger  than  the  elements  of 
a  subfabric,  the  specimens  become,  effectively,  less 
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homogeneous;  this  could  present  some  problems  in  optical 
and  scanning  electron  microscopy.   Conventions  for 
collecting  and  marking  oriented  rock  specimens  are  shown 
in  Figures  3-47  and  3-48. 

If  samples  can  be  withdrawn  from  boreholes  without 
damaging  the  soil  fabric,  some  consideration  should  be  given 
to  recording  the  fabric  through  360  degrees  using  acetate 
peels;  this  could  be  done  on  the  surface  of  the  core,  or 
better,  on  the  surface  inside  the  hole. 

Fabrics  should  also  be  studied  from  the  bottom  of  a 
section  to  the  top,  particularly  in  sections  where  some 
consolidation  has  taken  place.   This  could  demonstrate,  in  a 
practical  way,  the  sensitivity  of  fabric  analysis  to 
important  changes  in  material  properties . 

The  data  analysis  can  follow  at  least  some  of  the 
techniques  already  applied  successfully  in  conventional  studies 
of  rock  fabrics,  e.g.,  deformation  or  preferred  orientation, 
if  any,  of  elongated  grains  and  of  fractures.   In  addition, 
some  of  the  newer  techniques  used  in  studying  rock  fabrics, 
such  as  optical  diffraction  analysis,  can  be  used  to 
characterize  preferred  orientations,  elongations,  and 
concentrations  of  spatial  frequencies  (Figure  3-49  and 
Table  3.4;  Figure  3-50).   With  optical  diffraction  analysis, 
it  is  quite  a  simple  matter  to  compare  fabrics  over  the  whole 

range  of  scales  mentioned  earlier.   Much  of  the  photography 
mentioned  can  directly  provide  the  input  transparencies 
used  in  optical  diffraction  analysis. 

Quantitative  methods  are  available  for  characterizing 
the  information  in  the  transforms,  and  to  compare  them  with 
each  other o   For  example,  as  shown  in  Figure  3-51,  transforms 
of  progressively  deformed  material  can  be  compared;  changes  in 
the  width  of  crack  openings  and  the  initiation  of  fracturing 
at  higher  load  levels  leave  their  marks  on  these  profiles. 
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(a) 


Figure  3-47.  Methods  of  marking  oriented  hand 

specimens  (Turner  and  Weiss,  1963) 
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Figure  3-48»  Orientation  and  marking  of  hand  specimen 

and  thin  sections  (Turner  and  Weiss,  1963) 
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TABLE  3.4 

ELONGATIONS  OF  RECTANGULARLY 
STACKED  OVALS  (Pincus,  1969a) 


Input  ratio 

Tr 

ansform  ratio 

1.0 

0.99 

1.1 

1.08 

1.2 

1.23 

1.3 

1.30 

1.4 

1.40 

1.5 

1.49 

2.0 

1.92 

2.5 

2.50 

3.0 

3.05 
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a)   Elongation  =  1 


b)   Elongation  =  3 


c)   Elongation  =  3 
Figure  3-50a.   Rhombic  Stacking  of  ovals:   Inputs, 
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a)   Elongation  =  1 


b)   Elongation  =  3 


c)   Elongation  =  3 


d)   Calibration  Transform 


Figure  3-50b.   Rhombic  Stacking  of  ovals:   Transforms. 
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A  particularly  useful  region  of  the  profiles  is  between 
the  two  second-order  diffraction  spikes,  and  within  that  range, 
the  region  between  the  first  and  second-order  spikes  is  most 
useful.   This  corresponds  to  spacings  of  0.40-0.20  mm  (or 
spatial  frequencies  of  2.5-5  lines /mm)  in  the  inputs.   Note 
that  as  the  specimen  goes  from  no-load  to  the  first-loaded 
state  (b)  to  (c) ,  the  profile  becomes  smoother  and  the  profile 
level  in  the  region  of  the  first  order  spike  diminishes.   Going 
from  the  first-  to  the  second-loaded  state  (c)  to  (d),  the 
profile  level  remains  about  the  same,  but  the  curve  becomes 
rougher.   These  changes  are  consistent  with  the  closing  up 
of  cracks  under  initial  loading  and  the  onset  of  new  fractures 
with  increasing  loading.   Profiles  (b) ,  (c) ,  and (d)  are  oriented 
parallel  to  the  axis  of  the  cylinder,  which  is  also  the 
direction  of  compression. 

To  proceed  further  along  the  quantitative  route,  we 
have  devised  procedures  (Pincus ,  1973),  which  must  be  clearly 
understood  to  be  subject  to  further  refinement.   Figure  3-52 
represents  the  graphical  aspects  of  this  approach. 

First  note  the  similarity  between  Figure  3-52 (a)  and 
Figure  3-51,  and  recall  the  discussion  of  Figure  3-51.   In 
Figure  3-52  a) ,  we  have  smoothed  the  profile  in  the  region 
of  the  first  and  second-order  diffraction  spikes.   The 
a-curve  is  the  straight  line  connecting  the  ordinates  of  the 
smoothed  curve  at  the  positions  of  the  first  and  second- 
order  diffraction  spikes,  for  the  transform  at  load  a. 
The  b-curve  serves  the  same  purpose  for  the  transform  at 
load  b. 

Figure  3-52  (b)  shows  the  change,  from  an  a-curve  to 
a  b-curve  where  the  area  beneath  each  from  zero-frequency 
to  the  x-intercept  is  a  constant;  this  conservation  of  area 
(actually  total  spatial  frequency  in  this  interval)  is 
not  required  for  what  follows,  but  it  is  useful  analytically 
to  determine  whether  this  area  is  indeed  conserved,  say  in  a 
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deformation  series.   (The  area  would  not  be  conserved,  for 
example,  where  twinning  or  cleavage  develops  in  calcite 
grains  as  observed  in  our  experiments  with  Tennessee  marble; 
this  adds  a  high  frequency  component  to  spatial  frequencies 
already  present.) 

We  then  proceed  to  construct  some  indexes  based  on 
the  profile  in  the  interval  (I)  between  the  first  and  second 
order  diffraction  spikes. 

The  relative  change  in  area  from  the  a-curve  to  the 
b-curve  over  the  interval  I  is : 


AA 


m   1/2  (  a±  +     a2)/(  I  )  -  1/2  (  b±  +     b2)/(  I  ) 


or 


1/2   (fll  +  a2  )  /  (  I  ) 
AA   =   1  -  (F1  +  b^)/^  +  a2) 


where    a±     =  1/2  (a^   +  aj;  ,  b±  =     1/2   (b1  +  b[) ,  etc 

Reading  errors  for  small  values  of  b  or  a  can  be 
minimized  by  extending  the  b  and  a  curves,  reading  the 
intercepts  on  the  x  and  y  axes,  and  calculating  the  b  and 
a  values  at  the  first  and  second  order  spikes. 

When      b1>>   b2  and  a  »  J 

AA  »  1  -  (Fp/O^) 

When   AA  >  0,  there  is  a  decrease  in  area  over  I 
AA  =  0,  there  is  no  change  in  area  over  I 
AA  <  0,  there  is  an  increase  in  area  over  I 

In  Figure  3-52b) ,  AT  =  +  0.46.   In  this  example, 
areas  are  conserved.   The  shift  in  spatial  frequencies  is 
toward  higher  values  (+) . 

The  relative  change  in  slope  from  the  a-curve  to 
the  b-curve  is 
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A 


S  =  [(at  -  a2)/l]-[(b1  -  F2)/l] 
(ax  -  a2)/l 


or 


AS  =   1  -  (bL-  b2)/(a1  -  a2) 
when    b-.»b2  and  a,>>a2, 

AS  «   1  -  (EV)/(a.,),  which  is  also  the  approximate 

value  AA  under  the  same  conditions.   We  introduce  the  nota- 
tion 

AM  =  1  -  (F^/Ca^ 

to  identify  this  expression. 

In  Figure  3-52 (b)  ,  AS  =  0.84  and  AM  =  +0.52. 

When    AS^  >  0,  there  is  a  decrease  in  slope  over  I 
AS  =  0,  there  is  no  change  in  slope  over  I 
AS  *=  0,  there  is  an  increase  in  slope  over  I 

A  positive  value  of  AM  means  a  decrease  in  the 
first-order  coordinate;  a  negative  value  of  AM  means  an 
increase  in  the  first-order  coordinate. 

Finally,  Figure  3-52  (c)  presents  the  basis  for 
an  index  based  on  change  in  roughness  or  heterogeneity  of 
transform  profiles.   Here,  we  work  in  the  interval  between 
the  DC  spike  (zero  frequency)  and  the  first-order  spike. 
The  index  is  defined  as 

M     =      <Fa   '   V/(7a>    or     5E  "   !   -    (VFa> 

where   r   =  number  of  spikes  between  DC  (zero  frequency) 
a 

and  the  first-order  spike  with  height -"  some 
constant,  specified  fraction  of  the  height  of 
the  first-order  spike,  for  load  a. 
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r,    =  same  as  for  r   ,  for  load  b. 

b  a  ' 

r  =  mean  value  of  r   for  both  sides  of  the  profile, 
a  a 

r,  =  mean  value  of  r,  for  both  sides  of  the  profile. 


In  Figure  3-52 (c),  the  constant  fractional  value  of  the 
first-order  spike  is  taken  at  0.1,  and  AR  =  +0.75. 

In  our  work  in  rock  mechanics,  the  constant  fractional 
value  actually  used  is  10   , because  several  scale  changes  are 
necessary  to  chart  spike  amplitudes  for  the  very  intense 
calibration  transform  profile  and  the  transform  profiles  of  the 
rock  fabrics  studied.   All  of  the  transforms  of  rock  fabrics 
are  run  with  identical  instrumental  and  optical  settings. 

When    AR  >  0,  there  is  a  transition  to  a  smoother  profile. 
AR  =  0,  there  is  no  change  in  roughness. 
AR  <  0,  there  is  a  transition  to  a  rougher  profile. 

Changes  in  magnitude  but  not  sign  of  any  of  these 
indexes  are  easily  interpreted.   For  example,  suppose  that 
in  a  deformation  series  we  calculate  a  AR  for  the  first 
loading  point  with  respect  to  zero  load  and  find  it  to  be 
positive.   Then  suppose  that  we  calculate  the  AR  for  the 
second  loading  point  with  respect  to  zero  load  and  find  it 
also  to  be  positive  but  smaller  than  the  first  AR.   This 
means  that  the  transform  profile  becomes  smoother  from  zero 
load  to  first  load  and  less  smooth  from  the  first  load  to 
the  second.   All  of  the  foregoing  operations  can  be  digitized, 
thereby  increasing  the  volume  of  the  data  that  can  be 
handled.   Experiments  like  the  one  shown  in  Figure  3-53  can 
be  treated  in  this  way0   (Note  that  even  without  profiling, 
the  transforms  in  Figure  3-53 (b)  clearly  show  changes  in 
spatial  frequency  and,  at  some  levels  of  loading,  evidence 
for  preferred  orientation) • 
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Figure  3-53a.  Negative  photographs  of  series  of  uniaxially 
loaded  specimens  of  Tennessee  marble.   Strain 
increases  from  1-11  (Wawersik,  1968),  (Pincus, 
1969a). 
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Finally,  Figure  3-54  shows  that  the  specimen  surface 
need  not  be  planar  to  achieve  useful  results.   In  this  case, 
we  are  looking  at  the  surface  of  a  cylinder.   Cylindrical 
soil  samples  could  be  processed  in  the  same  way.   To  avoid 
propagation  of  projection  errors  in  calculating  indexes  such 
as  AA  or  AS,  it  is  advisable  to  generate  the  transform 

from  that  part  of  the  image  of  the  cylinder  in  which  the 
photographed  surface  departs  no  more  than  moderately  from 
perpendicularity  with  the  line  of  sight. 

3.8   FIELD  SAMPLING 

Any  sampling  scheme  should  be  designed  to  optimize 
the  information/cost  ratio  with  respect  to  the  collection  of 
samples  and  the  collection  of  data  from  in-situ  surfaces. 

The  same  principles  that  apply  to  location, 
recovery,  and  logging  procedures  in  collecting  samples  as 
a  part  of  an  engineering  soil  exploration  program  also  apply 
here,  in  general.   The  considerations  governing  location  and 
spacing  of  say,  informatory,  main,  and  supplementary  borings, 
and  of  test  pits,  shafts,  and  tunnels,  are  appropriate  here. 
These  considerations  can  be  thought  of  as  being  generally 
necessary,  but  they  will  usually  not  be  sufficient. 

In  sampling  for  fabric  studies,  orientation  considera- 
tions are  quite  important.   Thus,  completely  oriented  samples 
should  be  collected,  and  accordingly  it  might  be  necessary  to 
get  as  many  as  three  mutually  perpendicular  oriented  cores  at 
a  site  where  conventional  exploration  would  require  only  one 
unoriented  sample.   It  might  be  advisable  where  a  test  pit 
is  dug  to  collect  data  from  in-situ  surfaces,  to  put 
additional  effort  into  maintaining  planar  faces  and  in  trying 
to  cut  these  to  be  as  nearly  mutually  perpendicular  as  is 
practical.  Where  fabric  is  clearly  directional,  the  long 
axis  of  a  test  pit  or  the  axis  of  a  test  tunnel  might  better 
be  oriented  perpendicular  to  the  dominant  megascopic  fabric 
than  in  some  random  or  "convenient"  direction. 
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The  collection  of  the  rectangular  soil  monoliths 
used  by  pedographers  might  turn  out  to  be  much  more  efficient 
for  fabric  studies  than  are  the  "undisturbed"  core  borings 
of  engineering.   The  soil  fabrics  are  less  likely  to  be 
disturbed  by  monolith  sampling  and  the  mutually  perpendicular 
planar  faces  of  both  the  monolith  and  the  matching  in-situ 
faces  are  both  amenable  to  replication  and  interpretation. 
The  size,  shape,  and  orientation  of  the  monoliths  are  also 
more  readily  varied  than  are  the  corresponding  attributes  of 
conventional  undisturbed  soil  samples.   Typical  monolith 
dimensions,  1  to  2  ft  x  6  in.  x  4  in.  ( .30  to  .60  m  x  . 15  m  x 
.10m)  can  be  varied  to  say,  cubes  of  the  same  volume,  that 
is,  8.3  in.  (. 21  m)  or  6.6  in.  (.  17  m)  on  each  side,  where 
nothing  is  known  about  orientation  of  fabric  or  mechanical 
anisotropy.   Elongated  monoliths  could  be  cut  with  their 
long  axes  perpendicular  to  the  dominant  megascopic  fabric  or 
to  prominent  planar  heterogeneities,  or  in  some  direction 
(such  as  vertical)  that  might  have  some  special  significance 
for  soil  mechanics  property  testing  or  for  purely  design  or 
construction  requirements. 

In  determining  the  size  of  any  surface  studied, 
whether  this  be  in  samples  or  on  in-situ  faces,  it  appears 
wise  to  adopt  a  requirement  that  is  an  adaptation  of  ASTM 
practices  regarding  sample  sizes  of  soil  and  rock  specimens 
prepared  for  mechanical  testing.   That  is,  the  minimum 
dimension  or  span  of  any  surface  studied  should  be  at  least 
10  times  the  maximum  dimension  of  the  largest  fabric  element 
in  any  subfabric  of  interest  in  the  study. 

Further,  the  size  of  surface  to  be  studied  as  a  unit 
sample  of  fabric  should  not  be  extended  to  include  materials 
that  will  render  the  sample  heterogeneous,  unless  the  purpose 
of  extension  is  to  contrast  fabrics  (Figures  3-4  and  3-5). 
This  puts  a  maximum  limit  on  the  permitted  surface  dimensions. 
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Spacing  or  density  and  pattern  of  samples  and  of 
in-situ  surfaces  studied  will  depend  on  purposes  of  the 
investigation  and  characteristics  of  the  soil,  particularly 
on  its  heterogeneity  and  anisotropy.   Both  the  physical 
collection  of  samples  and  the  collection  of  fabric  data  from 
in-situ  surfaces  must  be  regarded  as  the  sampling  of  fabric 
data,  and  accordingly  must  be  coordinated. 

A  soil  with  planar  banding  (Figure  3-6),  for  example, 
should  have  smaller  spacings  perpendicular  to  the  banding 
than  parallel  to  the  banding,  if  banding  is  deemed  to  be  a 
legitimate  object  of  concern.   A  soil  with  curvilinear  banding 
may  require  a  curvilinear  sampling  grid. 

The  sampling  scheme  should  be  designed  three- 
dimensionally  to  sample  from  the  entire  volume  of  materials 
likely  to  influence  the  design  and  performance  of  the 
engineering  structure.   Ordinarily,  sampling  density  should 
vary  directly  with  heterogeneity  gradient,  but  many  exceptions 
to  this  guideline,  chiefly  for  practical  reasons,  will 
arise.   The  sampling  scheme  should  also  provide  for 
sequential  alteration  of  the  sampling  grid  to  optimize 
the  information  return  without  vitiating  the  basic  sampling 
design  and  without  increasing  the  cost  of  the  sampling 
operation;  the  most  important  ingredient  here  is  commonsense 
field  judgement. 

Collection  of  soil  samples  can  serve  the  dual 
purposes  of  fabric  analysis  and  mechanical  property  testing 
in  the  laboratory.   Collection  of  fabric  data  from  in-situ 
surfaces  has  the  advantage  of  providing  data  from  material 
less  likely  to  be  disturbed;  also,  much  data  can  be 
coordinated  with  data  from  in-situ  property  tests  and  with 
data  on  residual  stress. 
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3.9   RECOMMENDATIONS 

Expectations  for  determining  principal  stress 
directions  and  magnitudes  and  stress  and  strain  history 
have  been  discussed  in  3.2.   At  this  time  it  appears  that 
progress  in  refining  determination  of  principal  stress 
directions  is  likely  to  be  greater  than  determination  of 
absolute  stress  magnitudes.   Determination  of  relative 
magnitudes  of  principal  stresses  involves  about  the  same 
level  of  difficulty  and  is  in  fact  part  of  the  same  problem 
as  determining  principal  orientations . 

It  is  very  likely  that  some  situations  will  not  yield 
any  success  with  regard  to  any  of  the  three  major  objectives 
(directions,  magnitudes,  history).   It  should  be 
feasible,  with  additional  study,  to  spell  out  quite  explicitly 
the  conditions  that  will  permit  successful  determination  of 
any  of  the  three,  to  specified  accuracies. 

Some  experimental  work  that  might  superficially 
appear  to  be  quite  academic  could  have  very  important 
practical  payoffs.   For  example,  fabric  studies  in  con- 
junction with  in-situ  vane  shear  tests  would  be  very 
helpful  in  evaluating  quantitatively  the  strain  history 
along  a  variety  of  shear  surfaces  and  zones. 

In  contemplating  the  task  components  outlined  in 
Table  3.5,  each  listing  can  be  shown  to  be  actually  or 
potentially  useful.   It  is  most  unlikely  that  every  technique 
investigated  will  turn  out  to  be  practical  in  terms  of 
cost-effectiveness  and  some  techniques,  if  not  discarded, 
will  have  to  be  modified  to  be  practical.   Those  techniques 
that  do  prove  practical  should  be  classified,  in  the 
final  stages  of  the  work,  according  to  relative  reliability, 
cost  in  dollars  and  time,  and  genetic  relationship  to  standard 
soil  mechanics  tests  and  properties.   We  anticipate  that 
we  can  shed  light  on  the  basis  in  micro fabrics ,  for  example, 
for  the  values  of  Atterberg  limits  and  the  like. 
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Table  3.5  (Contfd.) 

MAN  -  YEAR  SUMMARY  RECOMMENDED 
FABRIC  ANALYSIS  STUDIES 


P.I.        A.I.        R.A.        Total 


Yr.  1 

1/2 

1/4 

1 

1  3/4 

Yr.  2 

1/2 

1/4 

1 

1  3/4 

YrG  3 

1/2 

1/4 

1 

1  3/4 

Yr.  4 

1/2 

1/4 

1 

1  3/4 

Yr.  5 

1/2 

1/2 

1 

Totals   2  1/2         1         4  1/2        8  * 


Additional  hourly  help  (technical  and  clerical  support) , 

an  average  1/4  -  3/8  man-year  per  year,  will  also  be  required, 
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Tasks,  objectives,  milestones,  and  estimated  net 
results  are  outlined  in  Table  3.6.   The  first  two-year 
unit  (Task  I)  will  concentrate  on  experimental  work  in  the 
laboratory,  with  some  minor  initial  field  work.   Task  II, 
the  second  two-year  unit,  will  overlap  the  first  unit  by 
one  year,  and  will  move  the  experimental  work  into  the 
field.   Task  III,  the  first  year  of  which  will  overlap  the 
second  year  of  Task  II,  will  last  two  to  three  years;  in 
addition  to  moving  the  work  directly  into  the  realm  of 
practical  operations,  Task  III  will  also  concern  itself 
with  increasing  resolution  and  decreasing  data  acquisition 
costs  in  terms  of  time  and  money.   Task  III  should  be 
coordinated  closely  with  operations  in  the  nonfabric 
components  of  the  program,  although  some  such  coordination 
should  prevail  at  a  meaningful  level  throughout  the 
project's  life. 

Of  course,  sufficient  flexibility  should  be  built 
into  the  program  to  allow  for  shifting  emphasis  and  for 
accommodation  to  exploit  breakthroughs  in  both  fabric  and 
nonfabric  components  of  the  program  and  to  utilize  advances 
developed  outside  the  program.   Year  5  is,  as  is  to  be 
expected,  the  most  uncertain.   This  could  range  from  0  to 
2  man-years,  depending  on  what  happens  in  years  1-4. 
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4.      ELECTROMAGNETIC  RESONANCE 

4.1  INTRODUCTION 

As  used  in  this  study,  electromagnetic  resonance 
(EMR)  is  defined  as  any  technique  that  exploits  the  absorp- 
tion of  electromagnetic  energy  by  resonating  centers.  A 
variety  of  resonant  centers  are  available  in  soil;  all 
absorb  energy  at  characteristic  frequencies.   Examples  of 
such  centers  include  polar  molecules  such  as  water,  nuclear 
magnetic  moment,  nuclear  quadrupole,  and  electron  spin. 
The  frequency  at  which  a  particular  center  resonates  depends 
on  its  type  and  environment.   The  environment  of  a  resonant 
center  may  include  crystal  lattice  distortion, the  distance 
separating  resonant  centers,  the  strains  imposed  on 
molecular  bonds  and  other  considerations  that  reflect 
gross  mechanical  stress  or  grain  orientation.  Electromagnetic 
resonance  is  separated  into  the  three  major  branches:   Nuclear 
Magnetic  Resonance  (NMR) ,  Electroparamagnetic  Resonance 
(EPR  or  ESR)  and  Nuclear  Quadrupole  Resonance  (NQR) . 

4.2  NUCLEAR  MAGNETIC  RESONANCE 

Nuclear  magnetic  resonance  is  for  the  investigation 
of  atoms  having  a  nonzero  nuclear  magnetic  moment „   The 
technique  measures  the  rate  at  which  the  magnetic  moment 
precesses  about  an  external  magnetic  field.   Approximately 
175  different  isotopes  of  the  elements  have  the  necessary 
nonzero  magnetic  moment  and,  in  theory,  are  observable  by 
NMR  spectroscopy..   Many  of  the  elements  that  are  found  in 
soil  samples  (Table  4.1)  theoretically  provide  suitable 
species  for  NMR  study.  There  are,  however,  many  factors 
that  must  be  considered  in  deciding  the  suitability  of  an 
isotope  for  NMR  investigation. 

The  important  parameters  used  in  determining  if  an 
NMR  study  of  a  particular  isotope  is  feasible  are  natural 
abundance  of  the  isotope  and  relative  sensitivity.   If 
the  values  for  these  two' parameters  (Table  4.1)  are 
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Table  4.1 

NMR-ACTIVE  ISOTOPES' 
FOUND  IN  SOIL  SAMPLES 


Isotope 

NMRb 
Frequency 
(MHz) 

Natural 
Abundance 

(%) 

Relative 
Sensitivity 

Magnetic0 
Moment 

Nuclear 
Spin 

h 

42.5759 

99.9844 

1.0000 

2.79268 

1/2 

13c 

10.705 

1.108 

1.59  x   10"2 

0.70220 

1/2 

14N 

3.076 

99.635 

1.01  x   10 '3 

0.40358 

1 

17o 

5.772 

3.7  x   10'2 

2.91  x   10"3 

-1.8930 

5/2 

27A1 

11.094 

100. 

0.206 

3.6385 

5/2 

29Si 

8.458 

4.70 

7.84  x   10"3 

-0.55477 

1/2 

25Mg 

2.606 

10.05 

2.68  x   10"3 

-0.85471 

5/2 

43Ca 

2.865 

0.13 

6.40  x   10"2 

-1.3153 

7/2 

aAfter  NMR  Table  by  Varian  Associates  5th  Edition  (1965) 
At  a  field  strength  of  10,000  Gauss 
cIn  units  of  Bohr  magnetons  (Bohr  magneton  =  eh/4irMc) 
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multiplied  any  isotope  with  a  combined  value  of  less  than 
0.1  are  considered  very  difficult  to  detect  with  current 
instrumentation,  excluding  fast  Fourier  transform  and 
computer  averaging  techniques.   With  computer  averaging 
and  FFT,  isotopes  with  a  combined  value  of  less  than  10 
become  undetectable. 

A  second  consideration  is  the  sample  and  sampling 
techniques.   Commonly  NMR  is  done  in  the  liquid  phase  to 
improve  the  homogeniety  of  the  magnetic  field  and  also  to 
decrease  the  magnetic  interaction  between  isotopes  which 
are  near  each  other.   IITRI  is  currently  engaged  in  perfecting 
a  high  resolution  NMR  technique   for  studying  solid  samples 
by  NMR.   This  is  an  attempt  to  extend  the  work  of  E.  R. 
Andrews  (1971),  who  has  done  extensive  work  on  the  technique. 
The  successful  completion  of  this  work  may  make  it  very 
attractive  for  the  study  of  species  found  in  soil  and  the 
interaction  of  these  species  under  stress.   At  present,  the 
broad  lines  one  observes  in  the  solid  state  limit  the 
amount  of  information  that  can  be  obtained  from  the  sample. 

It  would  be  desirable  to  be  able  to  investigate 
stresses  and  the  preferential  orientation  into  planes  of  the 
crystallites.   The  resulting  effect  on  the  nuclei  is  to 

(1)  distort  the  environment  in  which  the  nuclei 
occurs  and 

(2)  to  move  the  nuclei  relative  to  one  another. 

The  first  will  have  little  effect  on  an  NMR  spectrum 
characterized  by  broad  lines.   Much  'higher  resolution  will 

be  necessary  to  observe  the  result  of  small  perturbations  in 

the  field  of  the  nucleus  by  distortion  of  the  crystal  lattice. 

The  necessary  preliminary  work  on  such  a  high  resolution 

technique  is  being  done  at  IITRI . 

The  second  change  is  also  not  usually  observed  in 

1      19 
most  isotopes  other  than  H  and   F.   Cook  (1972)  reports 

that  the  line  width  "second  moment"  is  often  proportional 

to  r    the  inverse  cube  of  the  distance  separating  the 
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resonating  protons.   Since  this  distance  changes  under 
mechanical  stress,  NMR  may  be  used,  at  least  in  principal, 
to  measure  stress. 

The  largest  drawback  to  the  use  of  NMR,  especially 
for  in-situ  studies,  is  the  need  for  a  large  stable  magnetic 
field.   The  sensitivity  of  an  NMR  spectrometer  is  directly 
related  to  the  strength  of  the  magnetic  field  (i.e.,  the 
flux  measured  in  Gauss).   This  is  dictated  by  the  theory 
of  Boltzmann  distributions  and  the  conditions  for  observing 
NMR  resonance  transitions.   Suggestions  have  been  made  that 
the  earth's  magnetic  field  could  be  used,  however,  with  a 
field  strength  of  less  than  1  Gauss  as  compared  to  that  of 
the  average  NMR  spectrometer  (14,000  Gauss) ,  it 
is  easy  to  see  that  the  sensitivity  would  be  exceedingly 
low  and  that  detection  would  require  a  more  complex  elec- 
tronics package  than  is  now  feasible. 

4.3  NUCLEAR  QUADRUPOLE  RESONANCE 

The  nuclei  of  certain  atoms  behave  as  if  their  charge 
were  distributed  over  a  spheroid  instead  of  a  sphere.   This 
charge  distribution,  which  behaves  like  an  electric  quadru- 
pole,  will  experience  a  torque  in  an  electric  field  gradient. 
Such  gradients  are  produced  by  the  atoms'  electron  cloud  and 
by  charged  particles  near  by  in  the  crystal  lattice.   These 
provide  a  natural  field  about  which  the  quadrupole  precesses 
at  a  characteristic  Larmor  frequency.   Because  the  electric 
field  gradient  is  sensitive  to  the  distance  of  the  charges 
in  adjacent  atoms,  NQR  should  be  sensitive  to  volume  changes, 
hence  the  mechanical  stress  in  the  crystal.   Nuclear  quadru- 
pole resonance  is  similar  to  NMR,  except  that  no  external 
magnetic  field  is  required. 

Cook  reports  that  the  resonant  frequency  in  crystals 
usually  increases  with  pressure  and  notes  that  changes 
reported  in  the  literature  are  on  the  order  of  1  Hz/psi. 
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The  bulk  of  this  work  has  been  done  on  optically  clear 
gem-quality  semiprecious  stones  or  relatively  rare  minerals. 
Some  success  has  been  found  with  powdered  samples  if 
sophisticated  techniques  are  employed.   It  appears  that 
NQR  is  a  promising  tool  provided  techniques  applicable  to 
soils  can  be  developed. 

Therefore,  NMR  may  be  appropriate  for  future  use 
in  studying  stresses  in  soil  for  model  systems.   Continued 
efforts  in  the  field  of  high  resolution  solid  NMR  is 
imperative  and  may  possibly  yield  results  which  will  allow 
for  the  development  of  an  in-situ  method. 

4.4   ELECTRON  PARAMAGNETIC  RESONANCE 

Of  the  three  main  branches  of  magnetic  resonance, 
EPR  is  unique  in  that  it  detects  or  observes  the  electrons 
of  an  atom,  and  unlike  NMR  or  NQR,  can  directly  investigate 
the  effect  produced  by  changes  in  the  environment  of 
electrons,  the  actual  participants  in  molecular  bonding. 
This  property  will  be  useful  in  evaluating  the  strains 
imposed  on  selected  molecular  bonds.   It  is  expected  that 
empirical  relationships  can  be  found  between  the  soil 
effective  stresses  and  the  response  of  those  unpaired 
electrons  associated  with  some  molecular  bonds.   The  soil 
effective  stresses  contribute  to  forces  on  individual  soil 
grains,  resulting  in  distortion  in  molecular  bonds.  This 
distortion  of  resonant  centers  provides  a  potential 
mechanism  by  which  EPR  may  be  employed  to  measure  soil 
stresses  if  appropriate  centers  can  be  found  or  created. 
There  are  reproducible  differences  between  tetrahedral 
SiO/  in  quartz  tubing  and  in  kaolin  clay.   These  differences 
are  associated  with  the  differing  environments  in  which  the 
paramagnetic  species  is  situated 0  The  line  width  is  a 
measure  of  the  degree  of  crystallinitye  The  quartz,  being 
a  highly  crystalline  form  of  silica,  has  the  narrowest  line 
widths „  The  silica  in  kaolin  clay  is  probably  distorted 
from  pure  tetrahedral  by  being  included  in  a  crystallite 
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with  alumina „  The  presence  of  the  alumina  also  causes  the 
crystal  to  be  more  anisotropic,  an  effect  which  produces  a 
broader  EPR  resonance „ 

The  differences  between  various  forms  of  soil- 
related  compounds  such  as  calcium  carbonate  are  even  more 
pronouncedo   Many  similarly  related  forms  of  calcium 
carbonate  exist c   X-ray  powder  patterns  cannot  distinguish 
between  them,  however,  the  EPR  spectrum  of  Mn   ,  which 
substitutes  for  the  Ca   ions,  can  easily  identify  the 
various  forms.   In  addition,  EPR  is  not  limited  to  one 

phase.   Solid  samples  are  as  easily  investigated  as  liquids 

-12 
and  gases.   Very  small  quantities  (10   mole)  or  larger 

samples  (0.1  percent  in  20  grams  of  material)  can  be 

investigated. 

The  limitation  to  the  use  of  EPR  in  in-situ  experi- 
ments is,  as  with  NMR,  the  need  for  a  high-strength  magnetic 
field.   However,  most  transition  metal  ions  produce  spectra 
in  the  absence  of  a  magnetic  field  and  many  examples  of 
"zero-field  EPR"  are  reported  in  the  literature. 

Two  indirect  methods  also  exist  which  lend  them- 
selves to  model  studies  of  stresses  in  soils  and  to 
possible  use  in  in-situ  measurements.   The  first  is  the 
use  of  ionizing  radiation  to  create  paramagnetic  centers. 
This  method  could  be  applied  to  the  development  of  an 
in-situ  technique  utilizing  a  y-source  in  the  field  to 
create  paramagnetic  centers.  A  core  sample  would  be 
removed  and,  since  the  molecules  were  irradiated  under  a 
stressed  condition,  their  EPR  spectrum  might  be  indicative 
of  the  stresses  under  which  the  soil  was  subjected  and 
would  be  different  than  for  soil  similarly  irradiated  in 
the  absence  of  stress. 

The  second  technique  is  more  applicable  to 

I  | 
montmorillonite  soils „  A  paramagnetic  ion  such  as  Cu 

can  be  exchanged  with  all  exchangeable  ions  in  the  soil 

sample o   The  Cu   ion  is  then  investigated  by  EPR  and  the 
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difference  between  stressed  and  unstressed  soils  would  be 
indicated  by  differences  in  the  EPR  spectrum  of  the  Cu 
ions. 

Of  the  three  tools ,  EPR  affords  the  best  sensitivity 
to  the  kinds  of  differences  expected  for  stressed  and 
unstressed  soils 0   Experimental  studies  have  shown  that 
the  potential  for  developing  a  successful  method  of  deter- 
mining stresses  in  model  systems  is  higho   The  subsequent 
utilization  of  the  technique  for  in-situ  measurements 
depends  on  the  success  of  the  research  on  model  systems. 

4o5  EXPERIMENTAL  RESULTS 

In  an  attempt  to  establish  that  electron  paramagnetic 
resonance  (EPR)  spectroscopy  can  be  used  to  investigate 
resonances  arising  from  defects  in  soil  samples,  some 
preliminary  studies  were  conducted  on  ^-irradiated  kaolin 
clay0   The  EPR  of  kaolin  clay  (98  percent  pure)  was 
investigated  at  room  temperature  and  at  approximately  77 
degrees  K  in  its  unirradiated  state „   No  EPR  resonances  were 
found,  indicating  that  no  naturally  occurring  paramagnetic 
centers  associated  with  the  kaolin  crystallites  were  present „ 
The  lack  of  EPR  resonances  in  the  unirradiated  state  also 
indicated  the  lack  of  paramagnetic  impurities  existing 
substitutional^  in  the  kaolin  crystallites.   Since  most 
soil  deposits  are  not  98  percent  pure,  it  would  be  expected 
that  actual  soil  samples  would  contain  many  paramagnetic 

I  I         I  I         I  I         I  I  J L 

impurities  such  as  Mn   ,  Fe   ,  Co   ,  Ni   ,  and  Cu   either 

in  interstitial  or  substitutional  positions  in  the  crystal 

lattice  of  the  various  clays  and  inorganic  matter  (such  as 

calcium  carbonate)  from  which  the  soil  was  composed. 

fin 
The  kaolin  clay  was  Y-irradiated  by  two   Co  sources 

for  2-1/2  hours  at  a  flux  of  7500  rads/hour  for  each  source. 

The  total  dose  was  37,500  rads .   Ionizing  radiation,  such 

as  that  produced  by   Co  when  done  at  liquid  nitrogen 

temperatures  (77  degrees'  K)  and  in  a  vacuum  (50  p.  of  Hg 
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pressure),  produces  paramagnetic  ions  which  can  be  observed 
be  EPR  spectroscopy.   An  example  of  the  spectrum  obtained 
for  kaolin  clay  is  shown  in  Figure  4.1.   The  single  asym- 
metric resonance  at  g  =  2.0064  can  be  identified  with  a 
species  occurring  in  the  silica  portion  of  the  kaolin  lattice 
A  similar  EPR  resonance  in  y-irradiated  quartz  glass  has 
been  observed  in  previous  experiments  at  IITRI.   The  g-values 
and  line  widths  vary  with  the  kind  of  quartz  and  degree  of 
crystallinity.   Both  the  quartz  and  kaolin  clay  contain 
silicon  and  oxygen  atoms  in  a  tetrahedral  arrangement,  thus 
the  similarity  of  the  EPR  spectra  is  not  surprising. 

4.6  RECOMMENDATIONS 

Electron  paramagnetic  resonance  holds  potential  for 
investigating  the  environments  of  specific  entities  (ions, 
electron  bonds,  nuclei)  within  the  soil.   Little  is  known 
of  the  EPR  response  of  soil  except  that  it  is  influenced 
by  stress  state .   Considerable  basic  research  is  required 
to  quantify  soil  EPR  under  stress.   Baseline  soils  data, 
which  are  not  now  available,  must  be  acquired  by  laboratory 
experimental  work.   This  should  be  a  relatively  low- level 
long-range  program.   A  high  level  of  effort  is  not  justified 
because  of  the  lack  of  suitable  electronics.   Current 
hardware  is  bulky  and  unsuited  for  use  downhole  but  this 
drawback  will  be  overcome  as  electronics  packages  become 
smaller  and  EPR  techniques  more  sophisticated.   The  general 
plan  should  be  to  first  identify  paramagnetic  centers  in 
several  soil  types ,  such  as  has  been  done  for  silica  and 
the  ^-irradiated  kaolin.   Once  identified,  the  response 
of  these  centers  to  such  macroscopic  conditions  as  stress, 
orientation,  and  moisture  can  be  evaluated.   It  is  antici- 
pated that  approximately  one  man-year  of  effort  will  be 
required  for  the  initial  phase  and  that  subsequent  work 
will  depend  upon  the  results  obtained. 
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5.   RECOMMENDED  IN- SITU  STRESS  GAGE  DEVELOPMENT 

There  exists  immediate  need  for  hardware  capable  of 
measuring  the  full  stress  tensor  in  soil.   There  currently 
are  no  instruments  that  are  capable  of  making  rapid  determina- 
tions of  more  than  a  single  component  of  the  in-situ  stress 
tensor.   The  only  instruments  known  to  attempt  measurements 
of  multiple  stress  components  are  self-boring  stress  probes 
currently  under  development  in  England  and  France.   These 
devices  employ  high  modulus  interface  stress  gages  in  place 
of  pressuremeter  bladders.   They  are  otherwise  similar  to 
the  self-boring  pressuremeters  previously  discussed  in 
Section  2.4.2.   The  use  of  high  modulus  gages  results  in 
stress  relaxation  times  of  hours,  as  was  previously  shown 
in  Figure  2-39.   A  production  rate  of  one  or  two  tests  per 
day  may  be  acceptable  in  a  research  tool,  but  is  inappropriate 
for  a  production  testing  device.   Under  some  conditions  a 
self -boring  pressuremeter  is  a  rapid  stress  measuring  tool,  but 
it  will  be  shown  that  it  can  not  be  used  if  the  lateral 
in-situ  stresses  are  not  equal.   What  is  required  is  a 
measurement  system  capable  of  obtaining  the  in-situ  stress 
tensor  rapidly  for  a  wide  range  of  soil  types  and  stress 
conditions.   Speed  of  testing  means  that  significant  pore- 
pressure  changes  should  be  eliminated,  both  during  placement 
and  operation  of  the  instrument.   A  rapid  drained  test  results 
from  operating  at  vanishingly  small  strains,  using  low 
modulus  transducers  to  avoid  lengthly  stress  relaxation 
delays.   The  operating  parameters  during  insertion  must  be 
simple  in  theory  and  operation.   Sophisticated  development 
is  required  to  identify,  remove  or  automate  those  system 
elements  that  are  operator  sensitive,  while  maintaining  an 
acceptable  total  system  cost. 

5 . 1  CRITERIA  AND  CAPABILITIES 

Criteria  for  such  an  in-situ  stress  gage  are  as 
follows : 
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o  Capability  for  measuring  the  full  stress 
tensor  a,  ^  a„  f   cr^ 

o  No  inherent  limitations  to  operating  depth 

o  Applicability  to  a  wide  range  of  soil  types 

o  Measurement  accuracy  and  precision 

o  Relatively  inexpensive  in  terms  of  both 
time  and  costs 

o  Highly  skilled  operators  not  required 

These  criteria  and  the  state-of-the-art  suggest  that  a 
method  employing  direct  measurement  in-situ  should  be 
developed.   Regardless  of  the  care  with  which  an  "undisturbed" 
specimen  is  taken,  the  in-situ  stresses  aire  released  when  it 
is  removed  from  its  location,  and  stress  measurement  must 
then  depend  on  secondary  effects,  as  in  the  methods  of 
Skempton  and  of  Poulos  and  Davis.   The  disturbance  involved 
in  placing  a  probe  in-situ  should  be  reviewed  so  that  the 
measurements  can  be  made  directly  on  the  in-situ  stress  field. 

In  addition  to  the  development  of  new  hardware,  the 
precision  of  the  anisotropic  vane  shear  should  be  evaluated. 
The  simplicity  and  low  cost  of  this  technique  are  very 
attractive,  so  that  the  small  effort  required  to  determine 
its  operating  limitations  will  be  well  rewarded. 

5.2  ANALYSIS  OF  INSERTION  DISTURBANCE 

5.2.1  ANALYSIS  METHODS  AND  ASSUMPTIONS 

To  support  the  feasibility  of  the  recommended  develop- 
ment effort,  preliminary  analyses  were  conducted  to  evaluate 
the  disturbance  associated  with  probe  insertion  and  operation. 
While  low  disturbance  devices  have  been  developed  and  tested 
by  other  investigators,  (Baguelin  and  Jezequel,  Kenney,  and 
Wroth  and  Hughes) ,  none  was  specifically  designed  for  use 
in  a  full  triaxial  stress  field,  and  several  were  intended  to 
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measure  soil  stiffness  as  well  as  in-situ  stress.   These 
devices  indicate  that  probes  can  be  inserted  with  in- 
significant disturbance  in  fine  grained  soils  under  K  stress 
conditions . 

Six  quantities  must  be  measured  to  fully  determine  the 
full  stress  tensor.   Because  of  the  difficulty  in  making  a 
measurement  of  in-situ  shear  stress --recognizing  that  shear 
failure  occurs  on  the  transducer  surface  during  insertion  — the 
recommended  hardware  must  operate  by  measuring  the  normal 
stress  in  different  directions.   Specifically,  the  stress 
ellipse  in  the  plane  normal  to  the  axis  of  the  probe  was 
selected  for  measurement.   The  full  stress  tensor  may  be 
determined  by  three  measurements  made  with  the  probe  axis  in 
three  noncoplaner  orientations. 

The  two  dimensional  state  of  stress  in  a  plane  is 
given  by  the  stress  tensor 
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which  requires  the  measurement  of  three  independent  quantities 
This  tensor  contains  in-plane  principal  stresses  a        and  a 
regardless  of  the  orientation  of  the  plane.   In  general,  the 
in-plane  two  dimensional  principal  stresses  a   ,a   are  not 
equal  to  any  of  the  three-dimensional  principal  stresses 
aaa  unless  the  plane  happens  to  be  normal  to  a  principal 
direction.   The  full  stress  tensor  may  be  obtained  from 
measuring  o^    in  three  directions  and  solving  the 
simultaneous  equations  given  by  the  tensor  transform 

o    =  a .   a .   a .  .  . 
mn     im  jn   ij 

For  many  cases,  the  vertical  may  be  taken  as  a 

principal  direction  and  the  vertical  stress  as  a    =  vz 

zz 
For  these  conditions,  which  are  basic  to  the  design  of  all 

existing  low-disturbance  probes,  the  full  stress  tensor  is 
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determined  by  determination  of  the  two  dimensional  stress 

tensor  a  0  in  the  horizontal  plane. 

aft  r 

Our  analysis  is  based  on  the  assumption  that  the 

hardware  will  be  self-boring,  with  the  features  typical  of 

this  class  of  devices  (Figure  5-1) . 
These  include: 

o   smooth  outer  surface  with  no  projections 

o   soil  removal  through  the  interior  by 
slurry  fed  through  a  soil  chopper 

o   in-situ  stress  ahead  of  the  probe  maintained 
by  friction  on  the  interior  probe  walls  or 
by  slurry  pressure 

Much  of  the  analysis  to  follow  was  conducted  using 
the  finite  element  method.   This  is  a  departure  from  previous 
investigations  of  the  expanding  cavity  problem  (e.g.,  Vesic, 
1972;  Gibson  and  Anderson,  1961)  in  which  closed  form  analytic 
solutions  were  obtained.   The  FEM  was  used  for  the 
following  reasons : 

o  noncircular  geometry  may  be  investigated, 

o  general  loading  conditions  are  easily 
applied  , 

o   soil  constitutive  laws  can  be  changed,  and 

o   stress  and  strains  are  available  for  the 
entire  mesh,  whereas  closed- form  solutions 
are  usually  solved  for  the  response  at 
some  selected  point. 

An  elastic-plastic  finite  element  program  was 
employed,  but  the  plasticity  option  was  suppressed  for  these 
preliminary  analyses.   This  allowed  a  considerably  greater 
number  of  runs  and  is  justified  if  only  small  stresses  and 
strains  are  considered.   Because  the  intention  is  to  analyze 
a  low  disturbance  device,  large  stress  changes  are  to  be 
avoided.   Should  the  disturbance  become  significant  with  respect 
to  the  in-situ  stresses,  design  modifications  or  special 
operational  procedures  would  be  modeled  in  an  effort  to 
reduce  disturbance.   The  object  is  not  to  predict  large 
changes  in  stress  accurately,  but  to  avoid  them. 
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Section 


Interior  Geometry  and 
Soil  Chopper  Not  Yet 
Designed;  Design  Choices 
Must  be  Based  on  the 
Type  of  Instrumentation 
Employed  in  the  Section 
Above . 


In-Situ  Axial  Stress  a 


Note:   In  the  general  case,  the  radial  in-situ 

stress  will  vary  with  6,  and  in-situ  shear 

stresses  may  be  present  in  o  ,  o   and  an 

J  r  z   r      e 

Figure  5-1.   Components  Typical  of  Self  Boring  Devices 
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The  finite  element  mesh  employed  for  much  of  the 

analysis  is  shown  in  Figure  5-2.   A  comparison  of  radial 

deflections  made  with  the  closed  form  solution  for  the 

5     2 
case  a  =100  psi  (6.9xl05  N/m2)  ,  aKss?1  =   50  Psi  (3.4x10  N/m  ) 

shown  close  agreement  between  the  two  solutions. 
Contours  of  maximum  shear  stress  in  the  vicinity  of  the  hole 
were  obtained  from  the  FEM  solution  as  shown  in  Figure  5-4. 
Displays  such  as  this  are  helpful  in  visualizing  the 
stress  field,  b$t  are  not  easily  obtained  from  closed- form 
solutions . 

5.2.2   VERTICAL  STRESS  BELOW  PROBE 

The  soil  immediately  in  advance  of  the  probe  must  not 
experience  a  change  in  stress.   A  pressure  equal  to  the  in- 
situ  normal  stress  in  the  axial  direction  must  be  supplied 
by  the  axial  thrust  of  the  probe.   Depending  upon  the  design 
details,  this  may  be  supplied  either  by  interface  shear 
resistance  acting  on  the  interior  surface  of  the  probe,  or 
by  the  pressure  of  slurry  acting  on  the  top  of  an  intact  soil 
column  within  the  probe  mouth.   By  varying  these  parameters, 
the  soil  may  be  made  to  displace  either  inward  into  the  probe 
or  outward  away  from  the  probe.   The  three  groups  of  in- 
vestigators cited  have  shown  that  it  is  possible 
to  operate  a  low  disturbance  probe  in  such  a  fashion  that 
zero  radial  movement  occurs  in  advance  of  the  probe.   With 
this  assurance,  detailed  design  of  the  probe  interior  and 
operation  procedure  was  not  investigated.   This  analysis 
should  be  conducted  during  the  development  phase  so  that  an 
optimum  geometry  may  be  produced  that  will  be  effective  over 
a  full  range  of  soil  types.   It  is  anticipated  that  a  trans- 
ducer will  be  required  to  measure  the  axial  thrust  directly  so 
that  it  may  be  continuously  controlled.   This  may  take  the 
form  of  a  slurry  pressure  gage  or  a  load  cell  to  monitor 
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Part   A 

K " 

Figure  5-2. 


Finite  Element  Idealization  of  a  Part 
(Part  A)  of  a  Cylindrical  Cavity  in 
Any  Infinite  Soil  Mass 
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Shear  Stress  Contours 
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Figure  5-4, 


Maximum  Shear  Stress  Contours  Around  a 
Cylinderical  Hole  Subjected  to  Uniform 
Internal  Pressure 
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axial  thrust  on  the  probe  interior.   Detailed  FEM  analysis 
may  suggest  measurement  locations  on  the  probe  exterior  to 
monitor  anomalous  stress  concentrations  resulting  from  improper 
operating  parameters.   This  would  provide  a  direct  control  on 
the  disturbance  effects. 

5.2.3   DRIVING  SHEAR  STRESS 

During  insertion,  interface  shear  failure  conditions 
exist  on  the  outer  surface  of  the  probe.   Basic  research  on 
soil/interface  shear  resistance  (Leonards,  1965;  Potyondy,  1961) 
is  reviewed  by  Huck  et  al. ,  (1973).   The  results  are  summarized  in 
Table  5-1.   Using  Potyondy' s  data  we  can  express  the  interface 
shear  resistance  tt  as 

Tj  =   cfc  +  a  tan  (0  f^) 

Table  5-1  shows  that  the  angle  of  interface  shear 
resistance  6  can  be  as  high  as  0.8  0  for  polished  steel.   The 
data  of  Leonards  indicate  that  the  interface  shear  resistance 
may  be  reduced  by  a  graphite  coating,  but  a  teflon  coating 
is  ineffective  for  granular  materials. 

The  effect  of  the  interface  shear  stress  induced 

on  the  exterior  of  the  probe  is  shown  in  Figure  5-5  for  a 

normally  consolidated  soil  having  K   -  =   0.5  and  interface 

shear  resistance  y  =   0.5.   This  is  a  conservative  example 

in  that  it  is  equivalent  to  the  effect  of  y  =   0.25  in  an 

overconsolidated  soil  with  K   =  1  .   As  shown  in  Figure 

o  ° 

5-5 (a)  the  lateral  normal  stress  a   is  unchanged,  but  a 

shear  stress  x    =   a  /2  is  superimposed.   The  principal 

stress  directions  rotate  22-1/2  degrees  and  the  minor  principal 

stress  is  reduced  from  a   to  0.71  a  .   In  a  real  soil  dis- 

x  x 

playing  dilatational  behavior,  the  lateral  stress  a  will 
change  in  response  to  volume  changes  accompanying  shear. 
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Figure  5-5A.   Rotation  of  Principal  Stresses  Resulting 

From  Driving  Shear  Stress  Induced  on  Probe 
Surface 
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Figure  5-5B.   Distribution  of  Induced  Wall  Shear  Stress 
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The  magnitude  of  this  change  depends  upon  the  soil 
behavior  and  the  volume  of  soil  affected.   The  soil 
dilatational  behavior  may  cause  increase  or  decrease  in 
volume,  depending  on  void  ratio  and  other  properties.   The 
amount  of  soil  affected  is  known,  since  the  shear  stress 
decays  with  the  inverse  radius,  as  shown  in  Figure  5-5 (b) 
It  is  desirable  that  the  influence  of  the  induced  wall 
shears  be  reduced  or  eliminated.   This  may  be  accomplished 
by  materials  selection,  reducing  y  =   tan  6.   The  results 
of  reducing  y  are  shown  in  Figure  5-6.   For  a  value  of  y  = 
0.2,  the  minor  principal  stress  is  reduced  by  4%  with  an 
accompanying  rotation  of  11  degrees,  if  K  =0.5.   For 
K  =1,  [i  must  be  reduced  to  0.1  to  achieve  equal  benefits. 

5.2.4  LATERAL  PERTUBATIONS  WHILE  DRIVING 

It  is  expected  that  the  insertion  process  will  be 
accompanied  by  a  certain  amount  of  vibration  and  accidental 
lateral  displacements.   To  evaluate  this  effect,  a  finite 
element  model  was  assembled  to  determine  the  effect  of  a 
0.001-inch  (2.5x10"  m)  lateral  displacement  of  a  rigid  circular 
inclusion,  as  shown  in  Figure  5-7.   In  this  run,  the  node 
at  0 =   0   ,  at  the  top  of  the  circular  inclusion,  was  per- 
mitted to  move  freely  in  the  y -direction  because 
slippage  might  be  expected  at  this  location.   All  other  nodes 
on  the  inclusion  boundary  were  deflected  an  equal  amount,  A,  in 
the  y-direction.   The  slippage  allowed  the  single  node  in- 
duced stress  concentrations  in  the  adjacent  nodes.   In  an 
actual  soil  these  would  be  distributed.   The  stresses  on  the 
boundary  are  shown  in  Figure  5-7(b).For  the  conditions  in- 
dicated, the  radial  stress  change  produced  by  a  0.001  in. 

(2.5xl0~5m)  deflection  is  on  the  order  of  1.75  psi  (12xl03  N/m2) 
The  lateral  force  required  to  produce  this  displacement  is 

112  lb  (488N)  on  an  18-inch  long  by  3-inch  diam.  (,46m  x  ,076m) 
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Plate   Dimensions- 
=   o.33  Modeled  17"  x   34" 

_5        (.43  x    .86  m) 
=   0.001   inch    (2.54x10     m)    Hole  Radius    - 

Modeled   at    1    1/2" 
(3.8  x   10 -2m) 


-\*r 


19  Nodes  on 
Hole  Boundary 


rh     rh     rrr 
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Note:   Top  node  on  hole  boundary  allowed  to 
slip.   All  other  hole  boundary  nodes 
displaced  0.001  inch  (2.54  x  10_5m) 
in  x  direction. 


Figure  5-7A.   Rigid  Body  Motion  of  Hole  in  Semi-Infinite 
Plate 
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Resulting  From  Lateral  Deflection  of  0.001 
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probe.   The  behavior  on  rebound  depends  upon  the  soil  be- 
havior.  If  the  soil  unloading  modulus  is  twice  the  loading 
modulus,  the  stresses  resulting  from  a  displacement-relaxation 
cycle  would  be  equal  to  but  opposite  in  sign  to  those  in 
Figure  5 -7  (b)assuming  the  probe  is  forced  back  to  its  initial 
position. 

This  analysis  is  subject  to  the  assumption  that  the 
perturbation  stresses  are  small  relative  to  the  existing 
in-situ  stresses.   If  this  is  true,  then  the  soil  will  behave 
elastically,  and  the  stresses  will  be  everywhere  compressive 
with  no  debonding  from  the  probe  surface.   In  general,  it  may 
be  said  that  lateral  perturbing  forces  on  the  order  of  100  lb 

(4501$  may  be  acceptable,  yet  loads  on  the  order  of  1000  lb 

5    2 
(4500N)  will  cause  15  psi  (10  N/m  )  changes  in  stress,  which 

will  probably  not  be  acceptable  except  at  depths  where  the 

in-situ  stresses  are  high.   Sustained  vibrations  may  cause 

problems,  depending  upon  the  cyclic  loading  behavior  of  the 

particular  soil. 

5.2.5   RELEASE  OF  IN-SITU  SHEAR  STRESS 

A  very  low  interface  shear  resistance  is  desirable 

to  avoid  significant  driving  wall  shears  on  the  outer 

surface  of  the  probe.   If  a  completely  frictionless 

cylindrical  probe  is  inserted  in  a  triaxial  stress  field, 

the  in-situ  lateral  shear  stress  x   will  be  released  at 

xy 

the  probe  surface.   The  case  a   =  2a   =  a   =  100  psi 

5    2  x     y    z 

(6.90  x  10   N/m  )  was  evaluated  using  the  quarter  plane  mesh 

shown  in  Figure  5-2.   The  existing  in-situ  normal  and  shear 

stresses  on  the  probe  surface  are  shown  in  Figure  5-8 (a). 

Upon  full  release  of  the  in-situ  shear  stress  x   ,  the  in- 

xy 

situ  normal  stresses  at  the  probe  surface  change  as  indicated 
in  Figure  5-8  (b).   The  major  principal  stress  increases  by 
15%  and  the  minor  principal  stress  decreased  by  30%. 
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4xl05  N/m2) 


Figure  5-8A.   Normal  and  Shear  Stresses  on  a  Circular 
Locus  in  a  Biaxial  Stress  Field 
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Figure  5-8B.   Normal  Stress  Changes  Upon  Full  Release 
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It  is  likely  that  full  release  of  x   will  not  occur, 
but  the  shear  distribution  after  partial  release  is  not 
obvious.   If  the  partially  released  shears  are  assumed  to  have 
a  maximum  value  of  pa  p5  ),  but  to  be  unchanged  in  distribution, 
a  linear  interpolation  may  be  made  to  approximate  the  effect 
of  changes  in  the  interface  shear  resistance  ratio  y.   This 
is  shown  in  Figure  5-9,  together  with  the  effect  of  y  on  the 
driving  shear  stresses.    By  the  use  of  this  figure,  a 
qualitative  choice  for  an  optimum  value  of  p,  may  be  made 
of  0.2.   The  trends  acting  in  the  decision  are: 

o  For  larger  values  of  K  ,  a  lower  y 

is  desirable,  and  y   =   0  improves 

the  problem  of  driving  shears  for  any 

stress  ratio, 
o  For  equal  lateral  stress,  y  has  no 

effect  on  the  release  of  in-situ 

lateral  stress  t   which  is  zero  for 

xy 

a        =      o    .   If  the  lateral  stress  ratio 
x      y 

o   /o       =     K  is  increased  above  1,  the 
y  x 

effect  of  shear  release  is  negligible 
for  y>-  rfK. 

A  range  of  reasonable  choices  is  provided  by  the  maximum 
values  of  the  stress  ratios  K  and  a   /a    .      Additional 
analysis  is  required  in  the  developmental  program  to  in- 
corporate the  effects  of  dilatational  and  nonlinear  soil 
response.   This  is  basically  a  detailed  design  and  materials 
selection  problem  requiring  parallel  analytic  and 
experimental  analysis. 

The  entire  problem  of  selecting  y  may  be  side-stepped 
if  a  technique  is  found  to  provide  y-0  in  the  axial  direction 
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Figure  5-9. 


Effect  of  Wall  Friction  on  In-Situ  Stresses 

For  Assumed  Case  a      =  o      =   100  psi,  (6.9  x  10  N/nr) 

cr  =  50  psi,  (3.4  x  lO^N/m2). 
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and  |j=>0.5  in  the  circumferential  direction.   Possible 
techniques  by  which  this  might  be  accomplished  are  driving 
through  a  lubricated  membrane  (Figure  5-10)  or  departure 
from  circular  geometry  as  in  Figure  5-11.   The  inverting 
lubricated  membrane  would  require  a  mechanical  clamp  to 
release,  cut  and  hold  the  membrane,  and  only  a  limited 
number  of  measurements  can  be  made  in  any  given  hole  until 
the  stored  membrane  is  expended.   Both  schemes  require 
experimental  evaluation  to  determine  their  relative  advantages. 

5.3   STRESS  MEASUREMENT  CONCEPTS 

The  mechanism  by  which  biaxial  stresses  may  be 
measured  in  the  plane  normal  to  the  axis  of  the  instrument 
is  considered  in  the  following  paragraphs.   This  will  be 
accomplished  by  analysis  of  increasingly  elaborate  techniques, 
until,  at  the  end  of  the  discussion,  the  advantages  of  designing 
increasingly  sophisticated  devices  will  be  clearly  established. 
The  major  concern  is  how  accurately  a  biaxial  lateral  stress 
tensor  can  be  evaluated  without  knowledge  of  the  detailed 
constitutative  relations  of  the  soil. 

5.3.1  MENARD  PRESSUREMETER 

Early  pressuremeter  designs  required  first  that  a  bore- 
hole be  drilled  into  which  a  pressuremeter  could  be  dropped 
and  inflated.   The  normal  stresses  on  the  borehole  wall  are 
first  reduced  to  zero,  then  increased  hydrostatically.   The 
resulting  radial  deflections  are  plotted  in  Figure  5-12 
as  a  function  of  internal  pressure.   For  this  particular 
stress  ratio  and  Poisson's  ratio,  the  initial  deflections  are 
everywhere  inward  when  the  borehole  is  opened  and  the  boundary 
stress  drop  to  zero.   Pressurizing  the  pressuremeter  to  50  psi 
(3.4  x  105  N/m2)  and  100  psi  (6.9  x  105  N/m2)  results  in  the 
curves  shown.   The  linear  elastic  solution  employed  indicates 
large  disturbances  before  the  pressuremeter  is  lowered  into 
the  borehole.   The  use  of  this  device  requires  detailed 
knowledge  of  the  constitutive  relations  of  the  soil. 
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Rolled  Membrane  Storage 


Extruded  Membrane 


Figure  5-10.   Driving  Through  Extruded  Lubricated 
Membrane  Prior  to  Test 


cross  section 


Corrugation  Depth,  Height 
and  Interface  Shear  Resistance 
Optimized  to  Minimize  Driving 
Shear  Stress  and  to  Maintain 
In-Situ  Lateral  Shear  Stress 


Figure  5-11.   Non-Circular  Geometry  to  Provide 

Anisotropic  Interface  Shear  Resistance 
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Figure  5-12. 


Radial  Deflections  of  1.5  Inch  (3.81  x  10-2m) 
Radius  Borehole  For  Various  Uniform  Internal 

Pressures. 
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This  device  was  designed  as  an  analogue  to  foundation 
strength  and  deformation  characteristics  (Menard, 
1957;  Smith  and  Smith,  1968;  Baguelin  and  Jezequel,  1973) 
and  is  not  desirable  for  in-situ  stress  measurements. 

5.3.2   SELF-BORING  PRESSUREMETER 

The  disturbance  associated  with  opening  a  borehole  is 
eliminated  by  the  use  of  a  self-boring  probe.   The  stress 
field  about  such  a  device  may  be  somewhat  disturbed.   We 
shall  consider  the  response  of  a  self  drilling  pressuremeter 
to  the  stress  in  which  it  is  immersed,  regardless  of  disturb- 
ance.  The  basic  assumptions  are  that  the  in-situ  lateral 
shear  stress  is  not  released  and  that  inward  deflection  as 
large  as  0.001  in.  (2. 54x10"  m)  may  occur.   These  assumptions  are 
not  the  same  as  either  the  French  or  the  English  self-boring 
pressuremeters ,  but  do  approximate  the  response  expected 
of  this  general  class  of  devices. 

As  pressure  is  bled  into  the  device,  no  deflections 
occur  until  the  internal  pressure  exceeds  the  soil  normal 
stress  acting  on  the  membrane  at  some  point.   As  the  internal 
pressure  increases,  membrane  deflection  is  initiated  in  the 
direction  of  the  smallest  normal  stress  and  an  expanding  zone 
of  hydrostatic  pressure  develops  until  the  entire  membrane  is 
pressurized.   The  radial  deflections  of  points  on  the  probe 
boundary  are  shown  as  functions  of  the  internal  pressure  in 
Figure  5-13.  The  elastic  curves  corresponding  to  internal 

pressures  of  50  psi  (345xl03  N/m2) ,  76  psi  (525xl03  N/m2)  and 

3    2 
100  psi  (690x10  N/m  )  are  shown  in  Figure  5-14.   The  stress 

3    2 
conditions  for  PT  =  50  psi  (345x10  N/m  )  are  identical  to 

the  existing  in-situ  stresses  and  the  deflections  are  every- 
where zero0  An  approximation  in  the  analysis  results  from  the 
superposition  techniques  employed  in  this  analysis.   The 
deflections  shown  in  Figures  5-13  and  5-14  result  from  the 
wedge-shaped  stress  diagram  shown  cross-hatched  in  Figure  5-13 . 
This  is  the  difference  in  the  hydrostatic  internal  pressure 
under  the  membrane  and  the  existing  in-situ  stresses,  which 
were  supported  by  the  rigid  backing  of  the  membrane  prior  to 
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pressurization.   At  internal  pressures  between  the  values  of 
the  maximum  and  minimum  in-situ  normal  stresses,  the  normal 
stresses  are  slightly  decreased  outside  the  pressurized 
zone  as  the  soil  deflects  away  from  the  probe  body. 
Normal  stresses  cannot  be  reduced  to  a  value  less  than 
the  fluid  pressure  under  the  membrane,  which 

will  flow  into  any  region  having  normal  stresses  less  than  PT . 

3    2 
The  deflection  curves  for  PT   =   50  psi  (345x10  N/m  )  and 

3    2 
100  psi  (690x10  N/m  )  are  not  affected  by  this  approximation, 

since  the  membrane  is  not  pressurized  at  all  in  the  first 

case  and  is  fully  pressurized  in  the  second  case. 

This  type  of  device  ideally  indicates  the  normal 
stress  in  any  given  direction  by  the  level  of  the  internal 
pressure  at  the  time  the  membrane  initially  deflects  away 
from  the  probe  body.   The  minor  lateral  principal  stress  is 
easily  measured,  since  membrane  deflection  in  this  direction 
(point  A    Figure  5-13)  begins  when  the  internal  pressure  is  • 
equal  to  the  normal  stress.   The  major  lateral  principal 
stress  is  difficult  to  measure,  since  the  membrane  at  point 
B  may  begin  deflection  when  the  internal  pressure  is  either 
less  than  or  greater  than  the  normal  stress.   The  deflections 
at  this  point  are  affected  by  Poisson's  ratio  effects- 
This  uncertainity ,  or  the  converse  need  to  measure 
Poisson's  ratio,  may  be  eliminated  by  pressurizing  onlv 
small  zones  under  the  membrane. 

5.3.3   ZONE  PRESSURIZATION 

If  a  self  drilling  probe  is  placed  without  disturbance 
in  the  ground,  the  membrane  may  then  be  pressurized  over  a 
small, randomly  selected  zone,  say  30  degrees  or  less  of  the 
circumference  of  the  probe.   If  the  pressurizing  fluid  is  not 
allowed  to  escape  from  the  selected  area,  only  second- order 
effects  will  occur  outside  the  zone  pressurized.   The  mechanical 
details  may  be  arranged  in  a  number  of  fashions.   For  our 
later  discussion,  we  will  assume  that  this  is  accomplished 
in  the  manner  shown  in  Figure  5-15.   In  this  configuration, 
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Figure  5-15, 


Cut-Away  View  of  Hydraulic  Zone  Pressuriza- 
tion  Concept  for  30°  Pressurized  Zone 
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fluid  pressure  is  introduced  into  the  center  groove  and 
drained  by  the  edge  grooves ,  so  that  a  region  of  locally 
increased  pressure  is  established  in  the  vicinity  of  the 
center  groove.   If  the  fluid  pressure  is  increased  to  a  value 
slightly  above  the  in-situ  normal  stress  at  the  location  of 
the  center  groove,  the  membrane  will  lift  slightly,  allowing 

flow  to  occur  beneath  the  membrane  to  the  relief  grooves  at  the 

3    2 
side.   For  a  pressure  10  psi  (69x10  N/m  )  greater  than  the 

normal  stress  in  the  selected  direction  and  for  soil  properties 

E  =  104  (69xl06  N/m2)  and  v  =  0.33,  a  deflection  of  0.005  in. 

-4 
(1.3x10  m)  will  occur  on  the  pressurized  area.   Because  the 

area  outside  the  pressurized  zone  does  not  deflect,  the 

existing  normal  stresses  are  slightly  reduced.   The  stress 

distribution  around  an  arbitrarily  located  30  degree  zone 

3    2 
pressurized  to  20  psi  (138x10  N/m  )  above  the  existing 

in-situ  normal  stresses  is  indicated  in  Figure  5-16.   The 

solid  line  represents  the  existing  in-situ  stresses  around 

the  probe,  the  dashed  line  is  the  distribution  with  pressure  in 

the  zone,  and  the  dotted  line  shown  the  stresses  after  the 

zone  is  unpressurized  if  the  unloading  modulus  is  twice  the 

loading  modulus.   It  can  be  seen  that  the  normal  stresses 

more  than  20  degrees  away  from  the  pressurized  zone  are  not 

changed.   This  procedure  may  be  employed  to  measure  the  normal 

stress  at  up  to  seven  locations  around  the  probe  (Figure  5-17) . 

The  seven  zones  are  numbered  sequentially  in  order  of 

pressurization  in  Figure  5-17 (a)  ,  the  normal  stresses  at 

these  seven  locations  define  the  numbered  points  on  the  stress 

diagram,  Figure  5-17 (b).   From  this  diagram,  the  magnitude 

and  orientation  of  the  in-situ  stresses  are  easily  found. 

Note  that  measuring  the  average  normal  stress  over  a  30 -degree 

arc  tends  to  decrease  the  radius  of  the  indicated  Mohr's  circle 

without  shifting  its  center.   This  is  a  geometric  error  that 

is  easily  accounted  for  by  knowing  the  width  of  the  pressurized 

zones.   It  appears  that  if  the  membrane  deflections  can  be 

measured  accurately,  then  the  lateral  stress  tensor  may  be 

defined  with  precision  by  the  zone  pressurization  technique. 
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The  accuracy  of  membrane  deflection  techniques 

can  be  judged  by  considering  past  experience  with  borehole 

pressuremeters.   The  maximum  error  in  the  device  of  Baguelin 

and  Jezequel  (1973)  is  given  as  3.3x10   in.  (8. 5x10"  m) 0 

The  Wroth  and  Hughes  (1972)  device  employs  feeler  gages  to 

measure  radial  deflections  on  the  .039m  (1.5  in.)  diameter 

device.   The  accuracy  of  the  experimental  data  after 

_3 
differentiation  is  said  to  be  within  0.5%  strain,  or  3.8x10  in. 

(9.5x10  m)  on  the  radius.   The  raw  data  appears  to  measure 

-3  -5 

raw  deflections  to  within  2x10   in.   (5x10  m) „   Bronstein, 

et  al.  (1973)  suggest  that  radial  deformation  should  be 

-4        -3        -5 
measured  with  an  accuracy  of  4x10   to  6x10   in.  (10  m  to 

lo5xl0  m).   It  can  be  seen  that  quite  precise  measurements 

of  radial  deflection  are  possible.   In  the  transducer  design 

mentioned  earlier  in  this  section,  a  transducer  similar  to  a 

hydraulic  earth  pressure  cell  was  suggested.   The  flowrate 

in  such  a  device  increases  with  the  square  of  the  membrane 

deflection.   The  response  of  such  a  system  is  shown  in 

Figure  5-18.   It  is  assumed  that  small  pressure  levels  beneath 

the  membrane  cause  slight  penetration  of  the  membrane  into 

the  spaces  between  the  soil  grains 0   If  this  is  taken  as 

arbitrarily  equivalent  to  a  uniform  separation  of  0.001  in. 

-5  3  3 

(2.5x10  m)  ,  a  base  flow  rate  of  1.5  in.  /min  (25  cm  /min) 

is  established.   When  the  internal  pressure  increases  above 

the  normal  stress,  gross  membrane  deflections  take  place  in 

the  pressurized  zone.   The  flowrate  curves  that  occur  at 

three  different  locations  on  the  probe  depart  radically  from 

the  base  flowrate  when  the  normal  stress  is  exceeded.   At 

3    2 
10  psi  (69x10  N/m  )  overpressure,  the  flowrate  has  more  than 

doubled  even  though  the  effective  deflection  has  increased 

only  50%  over  the  assumed  base  rate.   The  pressure-flowrate 

curve  easily  indentifies  the  point  at  which  gross  membrane 

deflections  begin.   It  should  be  noted  that  this  transducer 

is  not  necessarily  the  best  that  can  be  conceived  but  it  is 

included  to  indicate  that  workable  designs  can  be  evolved  to 

measure  the  necessary  deflections  accurately. 
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5.4  PROBE  REQUIREMENTS 

An  in-situ  gage  for  use  in  a  biaxial  lateral  stress 
field  must  be  designed  to  perform  accurately  and  rapidly. 
The  recommended  design  includes  several  critical  elements: 

o  Detailed  design  of  the  interior  of  the 
probe  together  with  appropriate  in- 
strumentation and  operational  techniques 
will  insure  that  the  axial  in-situ  stress 
is  not  disturbed  by  the  approach  of  the 
device . 

o  Materials  should  be  selected  to  produce 
an  interface  shear  resistance  on  the  outer 
probe  surface  of  0.2  times  the  normal 
stress  or  less.   If  experimental  work 
indicates  that  grooving  or  other  minor 
departure  from  circular  geometry  will 
support  the  in-situ  lateral  shear  stresses, 
then  even  lower  friction  coefficients 
should  be  attempted. 

o  A  technique  (e.g.,  zone  pressurization) 
that  will  measure  in-situ  stresses 
without  a  knowledge  of  the  soil  moduli 
should  be  used.   If  pressuremeter 
techniques  are  to  be  used  instead,  major 
corrections  must  be  made  to  account  for 
the  effect  of  Poisson's  ratio  in  biaxial 
lateral  stress  fields. 

o  The  use  of  high-modulus  stress  transducers 
should  be  avoided  because  of" the  long  times 
needed  for  stresses  to  relax.   A  technique 
that  operates  by  stress  transfer  at 
vanishingly  small  strains  will  offer  the 
advantage  that  measurements  can  be  made 
very  rapidly  without  troublesome  pore 
pressure  effects. 

A  tentative  view  of  the  device  emerges  from  these 
requirements.   It  is  a  heavy-wall  self -boring  device  with  the 
instrumented  section  located  near  the  bottom.   The  outer 
surface  is  smooth,  of  a  material  with  lower  frictional  re- 
sistance than  polished  steel.   The  outer  surface  may  have 
shallow  grooves  to  increase  the  frictional  resistance  in  the 
circumferential  direction,  as  shown  in  Figure  5-11.   Alternative 
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treatments  of  the  interface  shear  resistance  problem  may  in- 
volve driving  through  an  inverting  membrane,  as  in  Figure 
5-10,  or  extruding  a  thin  film  of  lubricant  in  certain 
soils.   This  problem  can  be  best  evaluated  by  preliminary 
experimental  work.   The  pressure  transducers  will  be  small 
zones  spaced  around  the  circumference  of  the  device  that  are 
individually  pressurized.   Borehole  pressuremeters  universally 
use  rubber  membranes ,  but  a  better  choice  for  the  recommended 
stress  gage  may  well  be  thin  metallic  membranes.   This  choice 
is  permitted  because  of  the  small  maximum  deflections  measured, 
and  would  have  the  advantage  of  low  wall  friction.   This  again 
is  a  problem  for  design  detail  after  preliminary  laboratory 
investigation.   It  may  be  desirable  that  the  device  be  useable 
as  both  a  stress  gage  and  a  borehole  pressuremeter .   This  would 
be  possible  if  zone  pressurization  and  a  rubber  membrane  are 
used  in  that  the  stress  circle  could  first  be  traced  by 
pressurizing  individual  zones,  then  the  gross  soil  stiffness 
measured  by  pressurization  of  the  entire  membrane  in  pressure- 
meter  mode. 

5.5   PROPOSED  DEVELOPMENT  PLAN  FOR  IN-SITU  STRESS  GAGE 

The  detailed  desien,  fabrication  and_evaluation 
efforts  should  be  undertaken  with  the  goal  of  de- 
livering an  operational  in-situ  stress  gage  within  a  period  of 
18  months.   The  major  emphasis  of  this  effort  plus  reporting 
is  to  deliver  proven  hardware  for  operation  under  general  soil 
and  full  stress  conditions.   The  recommended  program  plan 
is  shown  in  Figure  5-19  . 

Phase  IA-Analysis 

Phase  IA  is  an  intensive  theoretical  analysis  of  the 
proposed  in-situ  stress  gage.   The  broad  design  will  have 
been  established  and  it  is  necessary  to  evaluate  the  performance 
of  the  device.   The  following  areas  should  be  treated: 

o   effect  of  soil  tVDe,  using  realistic 
constitutive  relations, 
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o    effect  of  lateral  stress  ratios, 

o    permissible  degree  of  soil  disturbance, 

o    instrumentation  requirements, 

o    optimum  gage  geometry  -  length,  diameter,  nose 
and  excavator  geometry,  transducer  sensing 
area  and  location, 

o    insertion  operational  parameters,  and 

o    test  and  data  reduction  procedures 

Phase  IB-Mechanics 

In  Phase  IB,  the  mechanical  subsystems  will  be  evaluated 
Models  of  the  various  subsystems  will  be  designed,  fabricated 
and  tested  in  the  laboratory.   The  following  areas  will  require 
treatment : 

o    reduction  of  shear  stress  along  the  outer 
soil-instrument  interface, 

o    real  time  measurement  and  control  of  the 

vertical  loading  imparted  by  the  interior  of 
the  instrument  to  the  underlying  soil  during 
insertion, 

o    elimination  of  other  disturbance  sources  such 
as  lateral  thrust,  vibration  or  others, 

o     driving  and  instrumentation  subsystems 

coupling  the  device  to  the  surface  hardware 
and 

o     automatic  control  of  testing  functions. 

Much  of  Phase  IB  should  consist  of  experimental  feasibility 
and  optimization  efforts  conducted  in  the  laboratory 

Phase  IC- Instrumentation 

In  Phase  IC,  candidate  transducer  designs  will  be 
evaluated  and  compared  with  the  requirements  developed  in 
Phase  IA.   The  following  technical  and  operational  areas 
should  be  treated: 

o     effective  transducer  modulus,  range,  accuracy, 
and  precision, 

o     data  retrieval, 

o     simplicity,  durability  and  cost  of  mechanical, 
transducer  and  recording  systems, 
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o    applicability  to  full  circle  or  zone 
pressurization,  and 

o    detailed  design  of  competitive  candidate 
transducer. 

Phase  IIA-Design  and  Fabrication 

The  Phase  IIA  effort  is  to  tie  together  the  three 
corroboratory  Phase  I  efforts.   During  the  design  and  fabrica- 
tion of  the  in-situ  stress  gage,  the  minor  hardware  required 
to  evaluate  the  anisotropic  vane  shear  technique  will  also  be 
fabricated. 

Phase  IIB-Laboratory  Evaluation  and  Modifications 

The  Phase  IIB  effort  is  a  major  laboratory  hardware 
evaluation  of  the  in-situ  stress  gage.   It  is  expected  that 
design  improvements  may  be  revealed  that  will  require  cycling 
through  a  modif ication-reevaluation  loop.   Testing  will 
include  the  insertion,  measuring  and  data  recording  functions 
under  realistic  laboratory  modeling.  These  tests  will  evaluate 
at  least  the  following  parameters: 

o  variations  in  the  magnitudes  and  the  ratios 
of  the  three  principal  stresses, 

°  ODeration  at  arbitrary  angles  to  the 
principal  directions, 

o  soil  type,  including  a  range  of  cohesive  and 
cohesionless  soils, 

o  degree  of  saturation  in  fine-grained  soils. 

These  tests  will  reveal  the  accuracy  of  the  device,  optimum 
operating  parameters,  conditions  to  be  avoided  and  useful 
software  for  data  reduction  and  presentation. 

While  these  evaluations  of  the  in-situ  stress  gage 
proceed,  a  parallel  but  much  smaller  evaluation  of  the 
anisotropic  vane  shear  technique  will  proceed.   This  will 
include: 

o   cohesive  and  cohesionless  soils, 

o   effect  of  testing  rate,  and 

o   the  possibility  of  using  data  at  large 
rotation. 
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The  anisotropic  vane  shear  evaluation  is  not  to  constitute 
a  large  portion  of  the  testing  program. 

Phase  III  Field  Testing 

Phase  III  constitutes  the  final  evaluation  of  the 
in-situ  stress  gage.   The  field  work  will  include  comparative 
testing  of  the  in-situ  stress  gage  and  at  least  one  existing 
technique.   The  following  areas  will  be  studied: 

o   accuracy  of  in-situ  stress  measurement. 

o   technical  operational  parameters. 

o   economic  factors  associated  with  the  use  of  the 
new  instrument. 

Reporting 

The  theoretical  and  experimental  efforts  will  be  fullv. 
documented  in  a  technical  report.   In  addition,  operational 
requirements  and  limitations  in  the  application  of  the  in-situ 
stress  gage  will  be  delineated  to  provide  the  practicing 
geotechnical  engineer  with  all  the  information  needed  to 
determine  when  and  how  to  use  the  instrument.   An  operator's 
manual  will  be  carefully  prepared  and  appended  to  the 
report.   This  will  include  step-by-step  operating  procedures, 
recommended  data  nandling  schemes,  or  nomographs  and  example 
problem  solutions. 

The  level  of  effort  recommended  for  each  phase  is 
shown  in  Table  5-2.   The  total  effort  is  expected  to 
amount  to  28  engineer-months  and  20  technician-months. 

In  addition  to  the  tentative  level  of  effort  as  shown, 
secretarial  time,  materials  costs  including  computer  time,  and 
rental  charges  on  standard  field  equipment,  should  be  added. 
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Table  5.2 

RECOMMENDED  LEVEL  OF  EFFORT  FOR  IN- SITU  GAGE 

DEVELOPMENT 


Phase 


Engineer 
(Months) 


Technician 
(Months) 


IA  Stress  Analysis 

IB  Mechanical  Design  Modeling 

IC  Instrumentation  Design 


4 
5 
3 


IIA  Design  &  Fabrication 
IIB  Laboratory  Testing 
III  Field  Evaluation 

Reporting 


3 
7 
4 

2 


3 
4 

4 


TOTAL 


28 


14 
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6 .   SUMMARY 

6.1  OVERVIEW 

A  number  of  techniques  are  currently  available  for 
use  in  the  partial  determination  of  in-situ  stress  fields. 
Existing  methods  measure  a  single  component  of  stress.  If 
it  may  be  confidently  assumed  that  a  given  site  has  a 
history  of  K  loading,  measurement  of  the  lateral  stresses 
and  the  assumption  that  the  vertical  stress  is  given  by 

a  =  vz  define  the  full  stress  tensor.   If  the  site  is  not 

z    ' 

level  or  has  experienced  events  to  modify  the  lateral 
stresses,  more  extensive  testing  is  required.   To  provide 
the  capability  required  by  the  latter  conditions,  a  three- 
level  recommendation  is  made: 

a)  Immediate  development  of  an  in-situ  stress 
gage  capable  of  determining  the  two-dimen- 
sional stress  tensor  in  the  plane  normal  to 
its  axis. 

b)  Medium-range  development  of  fabric  analysis 
as  an  effective  soil  mechanics  tool. 

c)  Long-range  development  of  methodology  using 
electroparamagnetic  resonance  as  a  soil  stress 
determination  tool. 

6.2  CURRENTLY  EXISTING  TECHNIQUES 

Current  techniques  are  classified  into  (a)  laboratory 
K  tests  (b)  laboratory  stress  tests  and  (c)  in-situ  stress 
tests.   The  first  classification  includes  a  number  of 
laboratory  test  methods  and  devices  by  which  K  loading 
conditions  may  be  imposed  on  a  sample  of  soil.   Use  of 
these  methods  requires  that  the  stress  history  at  the  site 
be  known  and  that  the  soil  respond  to  laboratory  testing 
in  a  fashion  similar  to  the  geological  behavior  in-situ 
(Section  2.3) . 

Two  examples  of  test  methods  by  which  the  actual 
in-situ  stress  is  estimated  by  laboratory  testing  were  found. 
These  are  the  Skempton  method,  which  employs  the  effect  of 
stress  on  the  pore  pressure  parameter  A,  and  the  lateral 
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consolidation  test  reported  by  Poulos  and  Davis.   These 
techniques,  both  of  which  are  applicable  to  fine  grained 
soils  only,  are  reviewed  in  Section  2.3.5. 

Five  in-situ  test  methods  were  discussed  in 
Section  2.4,  and  compared  in  Section  2.4.6.   The  most 
applicable  geotechnical  conditions  and  the  quantity  measured 
by  each  method  are  as  follows: 

Table  6.1 
CHARACTERISTICS  OF  VARIOUS  IN-SITU  TECHNIQUES 


test 


sonic  velocity 
(2.4.1) 


borehole 
pressuremeter 
(self  drilling) 
(2.4.2) 

borehole 

stressprobe 

(2.4.3) 

hydraulic 

fracturing 

(2.4.4) 

anisotropic 
vane  shear 
(2.4.5) 


most  applicable 
geotechnical 

conditions 

granular  or 
cohesive  soils 


best  in  cohesive 
soil 


frictional 
or  cohesive 


fine-grained 
KoSl 


best  in 

frictional 

soil 


quanity  measured 


stress  directions  by 
multiple  use;  estimate 
mean  stress  by  single 
use. 

mean  lateral  stress 
also  soil  stiffness. 


normal  principal 
stress  regardless  of 
orientation. 

minimum  principal 
stress  regardless  of 
orientation. 

stress  ratio. 


The  choice  of  technique  must  be  made  on  the  basic  circumstances 
associated  both  with  the  site  and  the  experience  and 
facilities  available  to  the  engineer.   We  cannot  include  all 
the  various  capabilities  and  limitations  of  each  technique 
and  the  table  is  not  a  substitute  for  sound  engineering 
j  udgment . 

6.3   RECOMMENDED  DEVELOPMENT 

In  Sections  3,  4  and  5  recommendations  were  made  for 
additional  development  of  fabric  analysis,  electroparamagnetic 
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resonance  (EPR) ,  and  the  in-situ  stress  gage.   The  use  of  EPR 
as  a  soil  stress  gage  is  a  long-range  recommendation.   The 
concept  shows  considerable  promise,  but  very  little  of  the 
basic  data  on  soil  response  is  known.   Because  of  the 
uncertainties  involved  in  converting  a  research  tool  to  a 
production  testing  device,  a  detailed  estimate  for  such  a 
program  cannot  be  reliable. 

The  development  of  fabric  analysis  as  a  soil  stress 
gage  is  a  medium- range  recommendation.   This  technique  may 
confidently  be  expected  to  resolve  the  principal  stress 
directions  and  the  relative  magnitudes  of  the  principal 
stresses  but  it  is  less  certain  that  the  absolute  magnitudes 
of  the  principal  stress  can  be  determined  with  any  degree  of 
accuracy.   Fabric  analysis  can  currently  be  used  with  some 
success  for  some  simple  cases  not  involving  rotation  of 
principal  stress  directions  or  other  complications. 

The  development  of  an  in-situ  stress  gage  is  the 
strongest  recommendation.   This  unit  will  be  capable  of 
resolving  the  two-dimensional  stress  tensor  in  the  plane 
normal  to  its  axis  in  a  single  operation.   This  is  not 
possible  with  existing  hardware.   For  conditions  under  which 
the  device  axis  can  be  oriented  along  the  direction  of  a 
known  principal  stress,  as  in  K  conditions,  the  full 
three-dimensional  stress  tensor  can  be  obtained.   In  the  case 
of  no  prior  knowledge  concerning  the  principal  stress 
magnitudes  or  directions,  three  orientations  of  the  axis 
in  oblique  planes  are  required  to  fully  resolve  the  full 
stress  field  in  a  region.   This  will  require  six  measure- 
ments using  existing  hardware.   In  addition  to  the  greater 
quantity  of  data  provided  by  the  in-situ  stress  gage,  the 
stress  data  are  expected  to  be  of  a  high  order  of  accuracy. 
Disturbances  in  a  full  triaxial  stress  field  are  given  more 
extensive  treatment,  and  the  proposed  measurement  techniques 
effectively  reduce  the  effect  of  soil  stiffness  by  the  use 
of  stress  transfer  methods.   Because  of  the  low  strain 
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levels  involved,  the  testing  time  in  saturated  soils  will 
be  reduced,  as  little  if  any  pore  pressure  will  require 
dissipation  during  the  test.   The  result  is  a  much  reduced 
testing  time  as  compared  with  pressuremeters ,  which  operate 

at  large  strain  levels,  or  stress  probes,  which  require 
lengthy  relaxation  times . 

The  recommended  system,  by  providing  economical 
and  accurate  in-situ  stress  data  will  permit  the  more 
confident  use  of  powerful  computational  tools  such  as  the 
finite  element  method. 
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APPENDIX  A 


SOIL  STRESSES 


A.l  STRESS  TENSOR 

Stress  is  defined  as  the  quotient  of  force  divided 

by  area,  having  the  units  pounds  per  square  inch  (psi) 

2 
or  Newtons  per  square  meter  (N/m  ) .   Geotechnical  engineers 

frequently  also  use  the  units  tons  per  square  foot  (TSF)  and 

the  dimensionally  incorrect  units  kilograms  per  square  centi- 

2  2 

meter  (Kg/cm  ).   This  report  will  deal  in  psi  and  N/m  . 

Stress  notation  employs  the  symbols  a  or  t  to  indi- 
cate individual  stress  components.   Subscripts  are  used  to 
indicate  the  direction  in  which  each  particular  component 
acts .   Either  of  two  notation  systems  may  be  employed  to  note 
the  normal  stress  acting  perpendicular  to  the  plane  of 
interest  and  shear  stresses  that  act  parallel  to  the  plane. 
In  the  second  (matrix)  notation,  repeated  subscripts  (xx) 
indicate  normal  stress  and  different  subscripts  (xy)  indicate 
shear  stress.   The  difference  in  normal  and  shear  stresses 
is  shown  in  Figure  A.l. 
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X 


deformed  shape 


original  shape 
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Normal  Stress  JL  to  y-z  Plane 


deformed  shape 
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I       f      T 


xy 


yx 


Shear  Stress    II     to    x-y  Plane 


Figure  A-l. 


Normal  and  Shear  Stresses  in  the 
x-y  Plane 
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The  concept  of  stress  is  typically  introduced  by  the 

consideration  of  forces  acting  on  an  imaginary  cube  cut  from 

a  solid  body  having  an  arbitrary  shape  and  external  loads 

as  shown  in  Figure  A -2.   In  the  absense  of  gravitational  or 

other  body  forces,  a  system  of  normal  stresses  o   .    g   .    a   . 
J  '  J  x'   y*   z* 

and  shear  stresses  t   ,  t   ,  and  t   are  seen  to  act  on  the 

xy*   xz'       yz 

cube.   The  magnitude  of  the  normal  and  shear  stresses  acting 
on  a  surface  depends  on  the  orientation  of  that  surface. 
Three  mutually  perpendicular  planes  exist  on  which  the 
shear  stresses  vanish  and  the  normal  stresses  attain  a 
maximum,  or  "principal"  value.   The  two  dimensional  case 
in  the  x-y  plane  is  represented  in  Figure  A -3 . 


If  matrix  calculus  is  employed,  it  is  convenient  to 
consider  stress  to  be  a  second  order  symmetric  tensor 
represented  by 


a .  . 


xx 


yx 
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xy 


yy 


zy 


xz 
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a 


zz 


In  this  notation,  the  normal  stresses  are  represented  by  the 
double  indices  a. 


xx 


,  a   ,  and  a   corresponding  to  the 


previously  developed  notation  o   .    o    «    and  a  .   The  shear 

f  J  f  X»       y>  Z 

stresses  are  represented  by  ctt___  =  a_ 
and  a  ^   =   a 


xy'   xz 


Txz> 


yx     "xy      xy'  "xz     zx 

t   .   The  matrix  is  symmetrical 
zt       zy     yz  J 

across  the  major  diagonal.   Finally,  the  eigenvalue,  or 
principal  stresses  become 
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Solid  of  Arbitrary- 
Shape  and  Loading 
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coordinate 
directions 


Figure  A-2.   Development  of  the  Stress  Concept 
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(b)   Principal  Stresses  in  the  x-y  Plane 
Figure  A-3.   Principal  Stresses  in  a  Plane 
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in  as  much  as  the  shear  stresses  vanish  in  the  principal 
stress  directions. 

A  convenient  representation  of  the  shear  tensor  in 
two  dimensions  is  the  Mohr's  circle  indicated  in  Figure  A-4. 
On  a  plot  of  shear  vs  normal  stress,  the  stress  tensor  is 
represented  by  a  circle  having  a  center  at  (a-,  +  Go)/2,0| 
and  radius  of  I  (a,  -  ao)/2  I  .   The  stresses  at  any  angle  a 
are  represented  by  the  coordinates  of  points  on  the  circle  at 
angle  2a  as  indicated  in  the  figure. 

Finally,  the  normal  stresses  about  a  point  may  be 
represented  by  an  ellipse  in  normal  stress  space  as  shown 
in  Figure  A-5.   The  normal  stress  components  at  an  inclina- 
tion a  and  a  +  90°  are  represented  by  the  G-  and  o~ 
coordinates  of  the  stress  ellipse,, 


263 


/^^                               \,(a' 

t) 

t^ 

CO 

Gl    "   °3        /]\ 

CO 
0) 

2        /              \ 

*-l 

4-> 
CO 

al  +  G3          /             1 

U 
CO 

cu 

|G3 

2       v/*a 

1CT1 

in 

a 
Normal  Stress    a 

Figure  A-4. 


Mohr's  Representation  of  the  Stress 
Tensor  in  Two  Dimensions 
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Figure  A-5. 


Stress  Ellipse  showing  Normal  Stresses 
Around  a  Point 
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A- 2.   EFFECTIVE  STRESS  CONCEPT 

Stress   tensors  in  geotechnical  applications  are 
frequently  represented  as  the  sum  of  the  effective  stress 
a1  and  the  pore  pressure  u  that  results  from  the  pressure 
of  fluid  within  the  soil  pores.   The  effective  stress  is 
associated  with  the  grain  to  grain  contact  forces.   The 
concept  of  total  stress  a   is  expressed  by 

a     -     a1  +  u 

The  effective  stress  concept  is  convenient  because  the 
Mohr  Rupture  theory  frequently  used  to  express  soil  shear 
strength  is  given  by 

S   =   c  +  af  tan0 

=   c  +  (a  -  u)  tan0 

in  which 

S  =  shear  strength 

c   =   cohesion 

0  =  angle  of  internal  friction 

Thus  the  effective  stress  concept  is  used  to  evaluate  the 
behavior  of  the  soil,  while  the  total  stress  determines  the 
forces  that  a  structure  must  resist. 

A. 3   STRESS  TRANSFORMATION  EQUATIONS 

To  deal  with  the  stresses  at  an  arbitrary  orientation, 
the  stress  transformation  equations  may  be  employed. 

In  matrix  notation 


a    =  a .   a .    a .  . 

mn       lm   jn    ij 
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where  a   and  a. .  are  the  stress  tensor 
mn      ij 

referred  to  two  coordinate  systems  having 
different  orientations  with  direction 
cosines  a^ 

Because  each  coordinate  system  is  composed  of  three 
orthogonal  axes,  the  direction  cosines  are  constrained  by 


a  Q  a    =   6Q  Kronecker  delta 

aB   aY       By 


a  a      a  D   =    6  Kronecker  delta 

aB   yB       aY 


The  usual  index  notation  and  range  conventions  apply. 

If  algebraic  notation  is  employed,  the  stress 
transformation  equations  for  coordinate  systems  xyz  and 
x'y'z'  became: 


2  2  2 

GX»X»    =   QXX  C°S  (X'>X)  +  CTyy  COS  (X  '  >  y)  +  QZZ  C°S  (X'>Z) 

+  2t   cos(x*,x)  cos(x',y)  +  2t   cos(x',y)  cos(x»,z) 

Ay  y  Z 


+  2t   cos(xf,z)  cos(xf,x) 


zx 


Tx'yf   =  °xx  cos(x'>x)  cos(y',x)  +  ayy  cos(x',y)  cos(y',y) 
+  Qzz  cos(x'>z)  cos(y',z) 

+  Txy   [cos(x',x)  cos(y',y)  +  cos(x',y)  cos(y',x)j 

+  Ty2   [cos(x',y)  cos(y»,z)  +  cos(x',z)  cos(y',y)l 

+  Tzx   [cos(x*,z)  cos(y»,x)  +  cos(xf,x)  cos(y!,z)] 

where  the  direction  cosines  are  given  by  cos (A, B) . 
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Similar  forms  are  used  for  the  remaining  four  stress 
components  in  x'y'z'  coordinate  system.   The  matrix  and 
algebraic  equations  are,  of  course,  mathematically  identical. 

For  the  case  in  which  the  stresses  are  measured 

on  an  arbitrary  plane  x'y1  in  a  a  a   a   stress  field,  the 

J    r  j  x   y   z  ' 

transformation  equations  give  a  •  • ,  a  •  • ,  t  •  •  in  the 

n       &    x 'x' *  y  y    x'y1 

plane,  and  a  ,  ,,  t  ,  ,,  T  f  f  acting  normal  to  the  plane, 
z  z     x  z     y  z 

It  is  convenient  to  formulate  the  equations  relating 
the  normal  stresses,  a  ,  in  arbitrary  directions  a-i  ,  a2j  ...  a. 
to  the  principal  stresses  a,  and  a«,  and  their  orientation 
B.   This  may  be  treated  by  reference  to  the  stress  circle 

a,  +  a«     a,  -  a0 

a   =   -i ^  +  -i ^  cos  2p 

n         o  o 


an      a   +  Aa  cos2S 
n       m 

where  a  and  Aa  are  the  mean  and  deviator  stresses, 
m 

If  measurements  of  a  are  made  at  angles  (a-.  4-  8),  (ao  +  8) 
...  (at  +  B) 

cr  i   =  a  +  Aa  cos   2 (a,  +  B) 

n  1      m  1   ■  ' 

CTn2  =   crm  +  Aa  cos   2(a2  +  3) 

a  .   =  a  +  Aa  cos   2 (a.  +  p) 
ni      m  -L 

In  this  system,  the  measured  a  and  a.  are  known,  the  a  , 
J  '  n      l  '      m' 

Aa  and  angle  B  are  found  by  simultaneous  solution  of 

three  equations.   In  the  body  of  this  report  it  is  suggested 

that  normal  stresses  be  measured  in  seven  directions.   This 

would  permit  a  least  squares  fit  to  be  made  on  the  values 

of  a  ,  Aa  and  B  since  up  to  35  sets  of  simultaneous  equations 

could  be  written  if  such  a  thing  were  desired. 
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A. 4   DISCONTINOUS  IN- PLANE  STRESSES  IN  A 
LAYERED  SEMI -INFINITE  BODY 


Let  Figure  A. 6  represent  a  layered  semi -in finite 
continuous  elastic  body,  isotropic  in  each  layer,  with 
layer  elastic  properties  E/-\,  v('\   anc*  layer  densities 
Y/.s  where  the  subscript  (i)  identifies  the  layer  in  the 
figure.   These  properties  and  densities  may  vary  from  layer 
to  layer  but  are  constant  within  each  layer.   Let  the  x, y, z 
axes  shown  in  the  figure  represent  a  cartesian  coordinate 
system  with  its  origin  at  the  top  surface  of  the  first  layer 
and  the  z-axis  in  the  direction  normal  to  the  layers,  with 

z.,  i  =  0,  1,  2,  ,  representing  the  boundary  locations 

between  the  layers.   Furthermore,  let  the  coordinate  system 
origin  be  located  at  an  arbitrary  point  on  the  top  surface 
and  the  x,y-axes  have  an  arbitrary  orientation  on  this 
surface. 

The  equations  governing  the  static  behavior  of  this 

(1) 
body  are  as  follows.   y   The  equilibrium  Equations  (la,b,c) 

where  X,Y,Z  represent  body  forces  per  unit  volume.   The 

surface  traction  boundary  conditions,  Equations  (2a,b,c), 

where  X,  7,  Z  represent  the  surface  tractions  and  l,m,n  are 

the  direction  cosines  of  the  external  normal  to  the  surface 

of  the  body  at  the  point  under  consideration.   Equations 

(3a,b,c)  and  (4a,b,c)  which  represent  the  stress-strain 

relations  and  Equations  (5a,b,c),  (6a,b,c)  which  are  the 

conditions  of  compatibility  of  strains  with  the  existence  of 

continuous  displacements. 


269 


Figure  A-6.   Layered  Semi-Infinite  Continuous  Elastic  Body 
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It  will  be  shown  below  that  an  in-plane  state  of 
stress,  discontinuous  from  layer  to  layer,  can  exist  in  the 
body  provided  the  following  sufficient  conditions  are 
satisfied: 

a)  The  in-plane  normal  strains  are  uniform  for  all 

layers.   That  is,  e  .  e  ,  are  constant,  independent 

x   y 

of  the  in-plane  orientation  of  the  x, y-axes  and 
independent  of  the  x,y,z  coordinates.   For  Fig.  1 

ex  =  ey  =  €xf  =  V  =  C;  Z  "  °        (?) 

where  c  is  a  constant. 

b)  Only  one  body  force  is  present :  that  due  to  gravity. 
•  Let  that  be  in  the  z-direction.   Hence, 

X  =  0;  Y  =  0;  Z  -  Y/.jj  z>0        (8) 

for  every  layer. 

c)  Only  one  uniform  surface  traction  is  present  at  z  =  0; 
that  normal  to  the  surface.   Hence, 


X  -  0;  Y  =  0;  Z  =  Z  (9) 

o 


where  Z   is  a  constant, 
o 

The  stress  state  resulting  from  the  existence  of  these 

conditions  is  derived  as  follows.   The  subtraction  of 

Equation  (3b)  from  (3a)  and  the  substitution  of  (7)  into  the 

result  gives 

ax  =  ay  =  a;  z  ~   °  (10) 

where  a  expresses  the  commonality  of  a  and  a  . 
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Since  a  does  not  vary  with  the  orientation  of  the 
z 

x,y-axes  it  follows  from  (3a, b),  (7)  and  (10)  that  a  is  a 
principal  stress.   It  follows  from  orthogonality  that  az  is 
likewise  a   principal  stress.   Hence, 

T   =  T   =  T   =0;   z  ^  0        (11) 

xy    xz    yz     '  ' 

The  substitution  of  (8),  (10)  and  (11)  into 
Equations  (la,b)  shows  that 

a  -  a  (z)  z  2=  0        (12) 

where  the  function  a  (z)  may  be  zero,  a  constant  or  a 
suitable  function  of  z  only  plus  a  constant.   The  substitu- 
tion of  (7)  and  (12)  into  (3a)  gives 

az  -  az  (z)  (13) 

i.e.,  a  is  a  function  of  z  and/or  a  constant  only. 
From  Equations  (1c),  (8)  and  (11)  we  obtain 

—  "  "  Y(1)         »!.!-  «-  «t        (14) 

i  =  1,2,3  -  — 

Equation  (13)  is  a  step  function.   The  function  is  continuous 
and  single  valued  within  the  open  interval  z.  .  <  z  -*=  z.. 
At  the  boundary  points  of  a  layer,  z.  •■  and  z.,  it  is  mathe- 
matically undefined.   It  must  therefore  be  defined  from 
additional  considerations  beyond  the  mathematical  ones,  such 
as  the  physical  requirements  of  continuity  of  displacements 
and  satisfaction  of  equilibrium.   A  sufficient  definition 
meeting  these  requirements  is 
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1Z5 

a  z 


m     -  Y 


(1) 


z  =  z 


(15a) 


Sz 


(i) 


z  =  z. 

1 


(15b) 


i  =  1,2,3,— 


By  combining  (15)  with  (14)  we  have  now  a  single  valued 
function  for  z  —  0  satisfying  all  the  required  conditions. 


az 


-  Y 


(1) 


z0-  z  -  zx 


(16a) 


da 


— =  -  Yd); 


z .  ,  -«=  z  «<  z . 
1-1     —   1 


i  =  2,3,4,  — 


(16b) 


From  integration  of  Equations  (16)  and  taking  into 
consideration  (9)  and  (13)  we  get 


a„  =  z   -  Y/i \  (z  -  z   )  ; 

Z      O      (1)   v       0  '   ' 


z  -«=  z  «<  z, 
o  -    -   1 


(17a) 


a  = 
z 


i-1 


z  -  v  Y/«\(z.  -  z.  - ) 
.2-,  (j)  J  J-l^ 
J=l 


-  Y/  •  \  (z-z.  -.)  ; 
T(i)v    l-l7' 


z.  -«:z^z.X17b) 
l-l        i  r 


i  =  2,3,4,  — J 


The  substitution  of  (17),  (12),  (10)  and  (7)  into  (3a)  gives 


1-v 


z  *  z  2  z,    (18a) 
o        I 


(1) 
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(l)       (l)  2 

1  "  V(D 


'i-1' 


z  -  zi    (18b) 


i  =  2,3,4, 


Equation  (18b)  representing  the  in-plane  normal  stresses  is 
a  single  valued  sectionally  continuous  function.   The  dis- 
continuities across  the  boundary  points  z.  come  about  because 
of  the  discontinuous  elastic  properties  E,.k,  V/.\  from 
layer  to  layer. 

It  can  readily  be  shown  by  substitution  of  the 
stress  functions  (18),  (17),  (11),  and  (9)  into  (2)  that 
the  boundary  condition  equations  are  satisfied.   An  example 
of  a  stress  field  based  on  the  above  formulae  is  shown  in 
Figure  A-7 . 

The  conditions  of  compatibility  (5)  and  (6)  are 
likewise  satisfied  as  can  be  seen  from  the  following.   The 
substitution  of  (17)  and  (18)  into  (3c)  gives 

2 


1-v 


+  2v 


(1) 


2vmc 


(i) 


l-v 


(1) 


1-V/  .  N  +  2v / . N 
(i)      (i) 


l-v 


(i) 


£+  mmm  £*  ——     "   «< 

O  1 


(i) 


l-v 


(i) 


19b' 


i-l 


z  ;f  z 


i  =  2,3,4, 


where    the   expression   for  ct      in    (19a)    is   given  by   (17a) 
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and  for  (19b)  by  (17b).   From  (4)  and  (11)  we  have 

=  0  ;     z  2!  0    (20) 


T 

= 

0; 

y 

= 

0; 

y 

xy 

zx 

7 

yz 

It  then  follows  readily  by  differentiation  of  Equations  (7), 
(19)  and  (20)  and  substitution  into  Equations  (5)  and  (6) 
that  the  compatibility  conditions  are  satisfied.   Hence  it 
has  been  shown  than  an  in-plane  state  of  stress,  discontin- 
uous from  layer  to  layer,  can  exist. 

The  satisfaction  of  the  compatibility  conditions 
insures  the  existance  of  a  continuous  displacement  field. 
That  this  is  the  case  here,  can  be  readily  shown  as  follows. 
Let  u,v,w  be  the  displacement  in  the  x,y,z  direction 
respectively  of  a  point  with  position  x,y,z  in  Figure  1. 
Then  : 


ow 


du  dv 

te~*s'  fy'*1"  Tz  =  *' 


(21) 


du       dv  _  du    ,   dw  dv    ,    dw  _ 

dy  "*~  dx  ==  7-~";  dz    [    Tc  ~  Tl1'         Tz  ~  T~y  ~  7vr  (22) 


2\dy      dz)-"*>         2\TZ~li)-^'        2\di~d~Tj)^^ 


(23) 


where  o>   o>   co  are  the  rotations  of  the  direction  of 
x'   y'   z 

principal  strains  when  moving  from  point  to  point  in  the 
body.   The  absence  of  shear  strains,  Equation  (20),  precludes 
such  rotation.   Hence 


o>x  -  0,  cuy  =  0;  a>z  =  0;       z  ±  0      (24) 

It  follows  by  addition  and  subtraction  between  Equations  (22), 
(23)  and  inclusion  of  (24)  that 
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hy 


-     0;  2H     =      0; 


2H 

£z 


£z 


21.     =     0 


(25) 


221     =      0 


c3w 
hy 


22.     =     0; 


Consequently   the  displacements   are   functions   of  single 
variables    only  as    follows 


u  =  u(x) 
v  =  v(y) 
as      =     w(z) 


(26a, b,c) 


Since  the  body  of  Figure  1  is  semi-infinite  we  can  select 
an  aribtrary  point  on  the  body  as  the  point  of  zero  dis- 
placement and  measure  the  displacements  of  all  other  points 
relative  to  it.   Let  that  point  be  x=y=z=0.   Then 


u  (0)  -  0 
v  (0)  =  0 
w  (0)   =  0 


(27a, b,c) 


In  view  of  Equations  (26)  and  (27)  we  can  use  (21),  (19), 
and  (7)  to  obtain  u,v,w  by  integration 
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y 
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y=0 
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V 


dy  = 


dv 

dy 


dy   = 


dy; 


z  ^  0 


y=0 


y=0 


(29) 


w  = 


^dz 


J 


dw 
dz 


dz  = 


e  dz 
z 


z  i  0    (30) 


z=0 


z=0 


z=0 


Hence 


u  =  ex 


z  *  0     (31) 


v  =  cy 


z  ^  0    (32) 
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ri-1  f 
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L  j=0  J 


Gzdz 


z . 
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ezdz 
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Z  •   1  «— «   Z  — »   Z  • 

1-1       1 


(33) 


i  =  1,2,3, 


As  can  be  seen  from  (31),  (32),  (33)  u,  v,w  are  continuous. 
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APPENDIX  B 
FABRIC  ANALYSIS  GLOSSARY 


ARRANGEMENT.   Arrangement  includes  distribution  and 
orientation. 

basic : arrangement  within  a  group  of  like  individuals. 

related: arrangement  between  groups  of  unlike  individ- 
uals  . 

referred : arrangement  of  individuals  with  regard  to 
specific  reference  feature. 

Arrangement  of  peds  at  each  level  includes  accomodation 
(degree  to  which  adjacent  faces  are  molds  of  each  other) , 
packing  (regular  to  random) ,  and  inclination  (orienta- 
tion with  respect  to  vertical  or  horizontal) .   (Brewer) 

BASIC  STRUCTURE.   The  structure  of  the  s -matrix,  that  is, 
the  size,  shape,  and  arrangement  of  simple  grains 
(plasma  and  skeleton  grains)  and  voids  in  primary  peds 
or  apedal  soil  material,  excluding  pedological  features 
other  than  plasma  separations.   (Brewer) 

CLOD.   (not  a  ped)  Caused  by  disturbance,  such  as  plowing  or 
digging,  that  molds  the  soil  to  a  transient  mass  that 
slakes  with  repeated  wetting  and  drying.   (Brewer) 

COLLECTIVITY  (of  fabric).   The  number  of  the  arrangements 
(collectives)  into  which  fall  the  fabric -elements 
distinguished.   (Sander) 

COMPONENTAL  MOVEMENT  (related  to  a  mechanical  deformation) . 
Any  movement  of  components  (in  rocks)  in  consequence  of 
which  a  rock  retains  its  continuity  but  acquires  a  new 
shape  in  the  period  of  time  considered.   (Sander) 

CONCRETION.   (not  a  ped)  Caused  by  local  concentrations  of 
compounds  that  irreversibly  cement  the  soil  grains 
together.   (Brewer) 

CRACK.   A  partial  or  incomplete  fracture.   (AGI  Glossary) 

CRYSTALLARIA .   Single  crystals   or  arrangements  of  crystals 
of  relatively  pure  fractions  of  the  plasma  that  do  not 
enclose  the  s -matrix  of  the  soil  material  but  form 
cohesive  masses;  their  morphology  (especially  shape  and 
internal  fabric)  is  consistent  with  their  formation 
and  present  occurrence  in  original  voids  in  the  enclos- 
ing soil  material.   Examples  include  crystal  tubes, 
crystal  chambers,  crystal  sheets,  and  intercalary 
crystals.   (Brewer) 
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CRYSTALLINITY.  (a)  The  degree  to  which  a  rock  (especially 
an  igneous  rock)  is  crystalline,  whether  holocrystal- 
line ,  glassy,  or  crypto-crystalline .  (b)  The  degree  to 
which  the  crystalline  character  of  an  igneous  rock  is 
developed  (e .g .,  macrocrystalline,  microcrystalline,  or 
cryptocrystalline)  or  is  apparent  (e  .g .,  phaneritic  or 
aphanitic) .   (AGI  Glossary) 

CUTAN.   A  modification  of  the  texture,  structure,  or  fabric 
at  natural  surfaces  in  soil  materials  due  to  concen- 
tration of  particular  soil  constituents  or  in-situ 
modification  of  the  plasma.  Gutans  can  be  composed  of 
any  of  the  component  substances  of  the  soil  material. 
(Brewer) 

DEFORMATION  LAMELLA.   A  type  of  slipband  that  is  produced 
especially  in  quartz  by  active  slip  within  a  mineral 
grain  during  tectonic  deformation;  also,  a  similar 
feature  produced  by  shock.   (AGI  Glossary) 

DOMAIN  (fabric) .   A  volume  defined  by  boundaries  such  as 
structural  or  compositional  discontinuities  within 
which  the  rock  fabric  is  uniform.   (AGI  Glossary) 

DYNAMIC  ANALYSIS  (fabric) .   Dynamic  inferences  concern  the 
nature  of  the  stresses  in  the  rocks  at  the  time  of 
deformations. 

Experimentally  deformed  specimens  are  compared  with 
their  naturally  deformed  counterparts  to  determine 
statistically  the  orientation  and  relative  magnitude 
of  the  principal  stresses  in  rocks  at  the  time  of 
deformation. 

Important  deformation  processes   include  fracturing, 
faulting,  gliding  flow,  rotation,  and  recrystallization. 
(Friedman) 

ELEMENTARY  STRUCTURE.   A  simplified  level  of  primary 

structure;  it  is  an  integration  of  a  charac t er is t ic 
size,  shape,  and  arrangement  of  specific  pedological 
features  and  the  basic  structure,  or  structure  of  the 
s -matrix.   (Brewer) 

FABRIC .   The  spatial  data  for  the  interior  of  the  domain 
under  consideration  describe  its  fabric.   (Sander) 

The  physical  constitution  of  a  soil  material  as 
expressed  by  the  spatial  arrangement  of  the  solid 
particles  and  associated  voids.   (Brewer) 

Fabric  is  the  element  of  structure  (size,  shape  and 
arrangement  of  solid  particles  and  voids)  which  deals 
with  arrangements.   (Brewer) 
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Texture-structure  and  degree  of  crystallization 
(crystallinity)  of  solid  particles.   (Brewer) 

After  Turner-petro  fabrics  is  the  study  of  all  struc- 
tural and  textural  features  of  a  rock  as  manifested 
in  every  recognizable  rock  element  from  crystal 
lattices  to  large-scale  features  requiring  field 
investigation.   (Friedman) 

A  fabric,  like  a  crystal  structure,  is  a  three- 
dimensional  array  of  discontinuities  or  heterogene- 
ities in  structure.   (Paterson  and  Weiss) 

(struct,  petrol.)  The  sum  of  all  the  structural  and 
textural  features  of  a  rock.   The  term  incorporates 
the  notion  of  function  or  behavior  (correlative  physical 
properties)  as  well  as  of  form  (arrangement  of  structural 
and  textural  components) .   (AGI  Glossary) 

FABRIC  ELEMENTS .   Actual  or  theoretical  like  parts  of  the 

fabric,  resembling  each  other  in  a  particular  respect. 
(Sander) 

Deformation  examples  include  crystallographic  parameters, 
deformation  lamellae,  kink  bands,  lineations,  fold  axes, 
fractures,  faults,  and  other  planes  of  mechanical  dis- 
continuity.  (Friedman) 

A  surface  or  line  of  structural  discontinuity.   It  may 
be  a  crystallographic  feature  such  as  a  lattice  plane, 
or  noncrystallographic,  such  as  foliation.   A  spatial 
array  of  a  particular  fabric  element  is  called  a 
sub  fabric  .   (AGI  Glossary) 

FAULT.   A  general  term  for  a  dislocation  in  a  crystal. 

A  surface  or  zone  of  rock  fracture  along  which  there 
has  been  displacement,  from  a  few  centimeters  to  a 
few  kilometers  in  scale.   (AGI  Glossary) 

FRACTURE.   (mineral)  The  breaking  of  a  mineral  other  than 
along  planes  of  cleavage.   A  mineral  can  be  described 
in  part  by  its  characteristic  fracture,  e.g.  uneven, 
fibrous,  conchoidal. 

(struct,  geol.)  A  general  term  for  any  break  in  a  rock, 
whether  or  not  it  causes  displacement,  due  to  mechanical 
failure  by  stress.   Fracture  includes  cracks,  joints, 
and  faults 

(exper.  struct,  geol.)  Deformation  due  to  a  momentary 

loss  of  cohesion  or  loss  of  resistance  to  differential 

stress  and  a  release  of  stored  elastic  energy. 
(AGI  Glossary) 
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FRAGMENT.   (not  a  ped)  Caused  by  rupture  of  the  soil  mass 
across  natural  surfaces  of  weakness.   (Brewer) 

GENEITY  (state  of) .   Characterization  of  the  position 
between  homogeneity  and  heterogeneity.   (Sander) 

GLAEBULE.   A  three-dimensional  unit  within  the  s -matrix  of 
the  soil  material,  and  usually  approximately  prolate 
to  equant  in  shape;  its  morphology  (especially  size, 
shape,  and/or  internal  fabric)  is  incompatible  with 
its  present  occurrence  being  within  a  single  void  in 
the  present  soil  material.   It  is  recognized  as  a  unit 
either  because  of  a  greater  concentration  of  some 
constituent  and /or  a  difference  in  fabric  compared  with 
the  enclosing  soil  material,  or  because  it  has  a 
distinct  boundary  with  the  enclosing  soil  material. 
(Includes  nodules,  concretions,  septaria,  and  related 
forms  but  not  the  spherulites  and  other  regular  crystal 
growths  whose  internal  fabric  suggest  formation  and 
present  occurrence  in  a  single  large  void) .   (Brewer) 

HETEROTACTIC  FABRICS.   Fabrics  in  which  subfabrics  do  not 
conform  to  a  common  symmetry.   (AGI  Glossary) 

HOMOGENEOUS  (mass) .   Exhibiting  essentially  the  same  physical 
properties  at  every  point  throughout  the  mass . 
(ASTM-D653-67) 

HOMOTACTIC  FABRICS.   Fabrics  in  which  all  subfabrics  are 
conformable  in  symmetry.   (AGI  Glossary) 

ISOTROPIC  (mass) .   Having  the  property  (or  properties)  in 
all  directions.   (ASTM-D653-67) 

JOINT.   A  surface  of  actual  or  potential  fracture  or  parting 
in  rock,  without  displacement;  the  surface  is  usually 
plane  and  often  occurs  with  parallel  joints  to  form 
part  of  a  joint  set.   (AGI  Glossary) 

KINEMATIC  ANALYSIS  (fabrics) .   Kinematic  inferences  concern 
the  displacements  that  have  transformed  the  initial 
fabric  into  the  observed  fabric .   It  is  assumed  that 
the  symmetry  of  the  rock  fabric  reflects  the  symmetry 
of  the  movement  responsible  for  the  evolution  of  that 
fabric .   (Friedman) 

MESOSCOPIC.   A  term  introduced  by  Dennis  (1967)  to  describe 
a  tectonic  feature  large  enough  to  be  observed  without 
the  aid  of  a  microscope  yet  small  enough  that  it  can 
still  be  observed  directly  in  its  entirety. 
(AGI  Glossary) 
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MICROSCOPIC.   Said  of  an  object  or  phenomenon  or  of  its 
characteristics  that  cannot  be  observed  without  the 
aid  of  a  microscope.   (AGI  Glossary) 

PED.   An  individual  natural  soil  aggregate  consisting  of 
a  cluster  of  primary  particles,  and  separated  from 
adjoining  peds  by  surfaces  of  weakness  which  are 
recognizable  as  natural  voids  or  by  the  occurrence  of 
cutans .   (Brewer) 

PEDALITY.   The  physical  constitution  of  a  soil  material  as 
expressed  by  the  size,  shape,  and  arrangement  of  peds. 
(Brewer) 

PEDOLOGICAL  FEATURES.   Recognizable  units  within  a  soil 

material  which  are  distinguishable  from  the  enclosing 
material  for  any  reason,  such  as  origin  (deposition 
as  an  entity) ,  differences  in  concentration  of  some 
fraction  of  the  plasma,  or  differences  in  arrangement 
of  the  constituents  (fabric)  (includes  features 
inherited  from  the  parent  rock  or  formed  by  the  pro  - 
cesses  of  deposition  of  transported  parent  material) . 
(Brewer) 

PEDOTUBULE.   A  pedological  feature  consisting  of  soil 
material  (skeleton  grains  or  skeleton  grains  plus 
plasma,  as  distinct  from  concentrations  of  fractions 
of  the  plasma)  and  having  a  tubular  external  form, 
either  single  tubes  or  branching  systems  of  tubes; 
its  external  boundaries  are  relatively  sharp.   Tubular 
form,  in  this  context,  means  that  the  feature  as  a 
unit,  or  its  impression  in  the  enclosing  soil  material, 
has  a  relatively  uniform  cross -sectional  size  and 
shape,  most  commonly  circular  or  elliptical;  that  is, 
the  impression  of  the  pedotubule  conforms  to  the 
definition  of  channels.   (Brewer) 

PLASMA.   That  part  of  soil  material  which  is  capable  of 

being  or  has  been  moved,  reorganized,  and /or  concen- 
trated by  the  processes  of  soil  formation.   It  includes 
all  the  material,  mineral  or  organic,  of  colloidal  size 
and  relatively  soluble  material  which  is  not  bound  up 
in  the  skeleton  grains.   (Brewer) 

PLASMIC  STRUCTURE.   The  structure  of  the  plasma  of  the 

s -matrix,  that  is,  the  size,  shape,  and  arrangement 
of  the  plasma  grains  and  associated  simple  packing 
voids .   (Brewer) 

PRIMARY  PEDS.   The  simplest  peds  occurring  in  a  soil  mate- 
rial; they  cannot  be  divided  into  smaller  peds,  but 
they  may  be  packed  together  to  form  compound  peds  of 
a  higher  level  of  organization.   (Brewer) 
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PRIMARY  STRUCTURE.   The  structure  within  an  apedal  soil 
material  or  within  the  primary  peds  in  a  pedal  soil 
material;  it  is  an  integration  of  the  size,  shape, 
and  arrangement  of  all  the  pedological  features 
enclosed  in  the  s -matrix  and  the  basic  structure  or 
structure  of  the  s -matrix.   (Brewer) 

S -MATRIX.   (matrix  including  skeleton  grains)  The  material 

within  the  simplest  (primary)  peds,  or  composing  apedal 
soil  materials,  in  which  the  pedological  features 
occur;  it  consists  of  the  plasma,  skeleton  g rains,  and 
voids  that  do  not  occur  in  pedological  features  other 
than  plasma  separations.   (Brewer) 

SCALAR  FABRIC  DATA  (undirected) .   Examples  include  individual 
fabric  element  considered  in  isolation;  about  pairs  of 
fabric  elements,  e.g.  fabric  habit  of  kinds  of  grain, 
shape  of  boundaries  between  fabric  elements  reacting 
with  one  another,  pore  volume,  percentage  of  particular 
species  of  grain  in  a  given  volume,  and  homogeneity 
(volume  of  smallest,  statistically  homogeneous  domain) . 
(Sander) 

SECONDARY  STRUCTURE.   The  size,  shape,  and  arrangement  of 
the  primary  peds,  their  interpedal  voids,  and  asso- 
ciated interpedal  pedological  features  in  a  soil 
material.   (Brewer) 

SHAPE  (particle) .   The  spatial  or  geometric  form  of  the 

particles  in  a  sediment  or  rock;  a  fundamental  property 
of  a  particle  that  determines  the  relation  between  its 
volume  and  surface  area.   It  depends  upon  the  sphericity 
and  roundness  of  the  particle,  although  the  term  is 
frequently  applied  to  sphericity  as  distinguished  from 
roundness.   Syn. :  grain  shape .   (AGI  Glossary) 

SKELETON  GRAINS.      Individual  grains   which  are   relatively 
stable   and   not   readily  translocated,    concentrated   or 
reorganized  by  soil-forming  processes;    they  include 
mineral  grains   and  resistant   siliceous   and   organic 
bodies    larger  than  colloidal   size.       (Brewer) 

SLIPBAND.   One  of  the  parallel  lines  known  as  Luders  lines 
or  Hartmann  lines  in  the  crystalline  grains  of  a 
material  stressed  beyond  its  elastic  limit  that  are 
visible  only  under  a  microscope  and  are  produced  by 
intracrystalline  slip.  .(AGI  Glossary) 

STRUCTURE  (soil) .   The  physical  constitution  of  a  soil 

material  as  expressed  by  the  size,  shape,  and  arrange- 
ment of  the  solid  particles  and  voids.  (Brewer) 
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SUBCUTANIC  FEATURE.   A  pedological  feature  (recognized  by 
a  difference  in  texture,  structure,  or  fabric  as 
compared  with  the  enclosing  s -matrix)  that  has  a  con- 
sistent relationship  with  natural  surfaces  in  the  soil 
material,  but  does  not  occur  immediately  at  the 
surfaces.   (Brewer) 

SYMMETRY.   The  presence  of  any  identities  in  such  positions 
in  space  that  they  replace  each  other  as  a  result  of 
symmetry-operations.  (Identities  are  color,  hardness, 
shape,  etc.).   (Sander) 

TERTIARY  STRUCTURE.   The  size,  shape,  and  arrangement  of 
the  secondary  peds  of  a  soil  material  (compound  peds 
resulting  from  the  packing  of  primary  peds),  their 
interpedal  voids  and  associated  interpedal  pedological 
features .   (Brewer) 

TEXTURE  (petrology) .   The  general  physical  appearance  or 
character  of  a  rock,  including  the  geometric  aspects 
of,  and  the  mutual  relations  among,  the  component 
particles  or  crystals;  e.g.,  the  size,  shape,  and 
arrangement  of  the  constituent  elements  of  a  sedimentary 
rock,  or  the  crystallinity,  granularity,  and  fabric  of 
the  constituent  elements  of  an  igneous  rock.   The  term 
is  applied  to  the  smaller  (megascopic  or  microscopic) 
features  as  seen  on  a  smooth  surface  of  a  homogeneous 
rock  or  mineral  aggregate.   The  term  structure  is 
generally  used  for  the  larger  features  of  a  rock. 
"Texture   should  not  be  used  synonymously  with  "struc- 
ture" although  some  textural  features  such  as  foliation 
or  flow  texture  may  parallel  major  structural  features 
of  a  rock.   (AGI  Glossary) 

The  physical  constitution  of  a  soil  material  as 
expressed  by  its  structure,  and  by  the  degree  of 
crystallization  (crystalinity)  of  the  solid  particles; 
fabric  is  a  part  of  structure,  which  is  a  part  of 
texture  in  this  context.   (Brewer) 

TROPY  (state  of) .   Characterization  of  the  position  between 
isotropy  and  anisotropy.   (Sander) 

VECTORIAL  FABRIC  DATA  (directed) .   All  fabric  data  which 
show  directional  (mostly  symmetrical)  influences  . 
Examples  include  preferred  orientation  according  to 
crystal  shape  and  crystal  structure,  directionally 
dependent  paths  for  transport  of  material,  response  to 
stress,  etc.,  and  growth  fabrics.   (Sander) 
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VOID.  Interstice;  an  opening  or  space  between  one  thing 
and  another,  as  an  opening  in  a  rock  or  soil  that  is 
not  occupied  by  solid  matter.   (AGI  Glossary) 

Space  in  a  soil  mass  not  occupied  by  solid  mineral 
matter.   (ASTM-653) 

Packing  voids  Voids  due  to  random  packing. 

Vughs  Relatively  large  voids  (other  than  packing) , 
usually  irregular  and  not  normally  interconnected 
with  other  voids  of  comparable  size. 

Vesicles  Voids  that  differ  from  vughs  in  that  their 
walls  consist  of  smooth,  simple  curves. 

Channels  Voids  that  are  significantly  larger  than 
those  which  would  result  from  normal  packing  of  single 
grains,  and  have  a  generally  cylindrical  shape. 

Chambers  Voids  that  differ  from  vughs  in  that  their 
walls  are  regular  and  smoothed  (metavoids)  and  from  both 
vesicles  and  vughs  in  that  they  are  interconnected 
through  channels  and  usually  have  a  characteristic 
shape. 

Planes  Planar  voids.   (Brewer) 
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APPENDIX  C 
ELECTRON  MICROSCOPY  SPECIMEN  PREPARATION 

A  critical  factor  in  the  application  of  electron 
microscopy  is  specimen  preparation.   Because  the  electron 
microscopy  views  the  sample  surface  directly,  more  careful 
preparation  is  required  than  is  the  case  for  other  tools  that 
respond  to  conditions  deeper  within  the  surface.   Specimens 
must  either  be  dehydrated  or  sealed  because  most  microscopy 
requires  that  they  be  exposed  to  vacuum.   Because  most  soils 
shrink  upon  dehydration,  special  preparation  techniques  are 
required . 

Barden  and  Sides  (1971)  discuss  disturbance  associ- 
ated with  the  following  specimen  preparation  techniques : 

o  peeling 

o  impregnation 

o  partial  impregnation 

o  grinding 

o  microtoming 

PEELING 

The  peeling  technique  requires  that  water  be  removed 
and  the  dry  clay  fractured  to  expose  an  interior  surface, 
then  peeled  to  remove  fracture  debris.  Water  removal  may  be 
accomplished  by  air  drying,  freeze  drying,  or  critical-point 
drying.  Air  drying  is  fastest  and  easiest,  but  may  cause 
shrinkage.   Critical-point  drying  reduces  shrinkage  due  to 
the  lack  of  liquid-gas  meniscus  forces  at  the  critical  point 
of  water,  which  is  3200  psi  at  706  F.   Freeze  drying  and 
critical-point  drying  are  discussed  by  Sridharan,  Altschaeffl 
and  Diamond  (1971) . 
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Once  dry,  the  specimen  is  fractured  to  expose  a 
surface  along  the  desired  plane,  and  the  fracture  debris  is 
peeled  off  by  50  to  100  applications  of  adhesive  tape  to 
remove  surface  debris. 

IMPREGNATION 

A  number  of  products  have  been  suggested  to  impreg- 
nate microscope  specimens.   The  intent  is  to  replace  (at 
least  partially)  the  pore  water  with  a  substance  that  will 
cement  the  clay  grains  together.   Depending  on  the 
permeability  of  the  soil  and  the  viscosity  of  the  cement, 
a  long  time  may  be  required  for  impregnation. 

PARTIAL  IMPREGNATION 

The  partial  impregnation  procedure  attempted  by 
Barden  and  Sides  (1971)  involved  impregnating  only  the  first 
2  or  3  mm  of  soil  below  the  specimen  surface.   After  harden- 
ing, the  untreated  clay  was  washed  away.   However,  a  zone 
of  partial  impregnation  was  found  ahead  of  the  fully  impreg- 
nated zone,  making  the  technique  unreliable  in  its  present 
form. 

GRINDING,  MICROTOMING  AND  ETCHING 

Grinding  was  used  on  both  impregnated  and  peeled 
specimens,  but  found  to  cause  re-orientation  of  the  surface 
grains.   Microtoming  was  discarded  without  trial  because 
of  previously  reported  difficulties  and  the  presence  of  silt- 
sized  particles  in  most  natural  clays.   Etching  was  used 
with  some  success,  but  the  open  structure  of  highly 
flocculated  clays  may  be  damaged. 

CONCLUSION 

In  general,  the  peeling,  impregnation  with  gelatine, 
and  gelatine  impregnation  with  etching  were  found  satisfactory 
for  optical  work.   Corbowax  impregnation  allowed  surface  damage 
and  was  unsatisfactory  for  the  scanning  electron  microscope. 
It  was  adequate  for  x-ray  diffraction  which  penetrates  30  to 
50  microns.   A  review  of  Table  D.l  indicates  that  the  most 
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TABLE  D-l 

USE  OF  VARIOUS  SAMPLE  PREPARATION  TECHNIQUES 
FOR  INVESTIGATIONS  OF  SOIL  MICROSTRUCTURE 
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common  technique  used  is  impregnation.   Since  this  is  a  time 
consuming  process,  the  peeling  technique  is  attractive  when 
air  drying  shrinkage  is  not  prohibitive. 
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